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foreword 


\ 

The  Conjugate  Point  Symposium  was  held  at  Boulder,  Colorado,  during 
June  13-16,  1967.  The  symposium  concerned  upper  atmospheric  phenomena 
whose  manifestations  uniquely  identify  pairs  of  northern  and  southern 
hemisphere  locations  at  the  earth's  surface  which  are  connected  by  the 
earth's  field  lines.  The  papers  presented  at  the  meeting  were  divided 
among  four  sessions:  Luminosity.  Ionosphere  and  Particle  Precipitation, 

EM  and  HM  Guided  Waves,  and  Geomagnetic  Field  and  Magnetosphere. 

This  Technical  Memorandum  is  a  reproduction  of  the  papers  submitted 
by  the  speakers.  In  a  few  cases  only  the  abstract  is  available.  The  pages 
are  numbered  for  the  session  and  paper  order  with  some  changes  from  the 
original  order  to  improve  continuity.  The  reproduction  in  this  Memorandum 
does  not  constitute  formal  publication;  it  is  done  as  a  service  to  the 
attendees  and  the  symposium  sponsors.  Primary  responsibility  for  the 
technical  content  of  this  document  must  rest,  of  course,  with  the  individual 
authors  and  their  organizations.  Readers  should  refer  to  these  papers  as 
presentations  to  the  Conjugate  Point  Symposium. 
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Rapid  Geomagnetic  Field  Vnriat.i<  j  ::  r;b served  at 
Conjugate  locations 

* 

by 

W.  H.  Campbell 

Institute  for  Telecommunication  Sciences  and  Aeronomy 
Environmental  Science  Services  Administration 
Boulder,  Colorado 

1 .  Introduction 

This  paper  reviews  the  observational  evidence  for  the  earth  field 
line  dependence  of  rapid  geomagnetic  field  variations  in  the  frequency 
range  from  4.0  to  0.002  c/s.  In  the  literature  of  geomagneticians  the 
"micropulsations"  have  been  referred  to  using  the  letters  £  for  "pulsation" 
followed  by  t,  c,  £,  or  £  for  the  amplitude  trace  characteristics  of 
"transient",  "continuous",  "giant"  (Ann.  IGY,  1959)  or  "pearl"  (Sucksdorff, 
1936)  appearance.  In  1963  8111  international  committee  (Gen.  Assem.  IUGG, 
1963)  recommended  a  £i  and  £c  system  of  naming  for  *the  field  fluctuations, 
using  i^  to  indicate  an  irregular  pattern  and  c  a  regular  mainly  continuous 
appearance  the  amplitude  trace.  A  suffix  1  to  5  is  proportional  to  the 
pulsation  period  (fig.  l). 

This  later  naming  will  provide  an  order  for  presenting  the  topic  of 
this  paper.  The  £i  will  be  given  first,  then  the  £c  following  from  lower 
to  higher  frequencies.  The  review  will  emphasize  the  important  character¬ 
istics  of  the  micropulsations  which  have  been  observed  simultaneously  at 
comparable  earth  field  locations  in  opposite  hemispheres  (Campbell  and 
Matsushita,  1967).  Table  1  summarizes  the  station  locations  of  the  authors 
whose  work  is  discussed  herein.  A  review  of  the  general  properties  of 
geomagnetic  pulsations  has  been  given  by  Campbell  (1967b).  With  a  focus 
upon  the  phenomena,  there  will  be  no  attempt  to  evaluate  priority  of  first 
observations. 


2.  Pulsations  with  Irregular  Forms 


This  class  of  pulsations  has  frequency  components  spread  through  the 
entire  range  of  4.0  to  0.002  c/s.  They  are  known  to  be  closely  related 
in  occurrence  and  amplitude  to  auroral  zone,  disturbed  ionospheric 
contions.  It  is  usually  difficult  to  indicate  a  distinct  separation 
between  the  lower  £i  frequencies  and  types  which  typically  occur  simul¬ 
taneously  with  the  rapid  variations  such  as  the  magnetic  disturbances  of 
SSC  (storm  sudden  commencement),  positive  and  negative  bays,  SI  (sudden 
impulses),  etc. 

Figure  2  illustrates  for  College,  Alaska,  and  Macquarie  Island  some 
of  the  similarities  and  difference  observed  in  jri  1  amplitude  records  on 
December  28,  1961  (Campbell  et  al.,  1962).  The  top  two  records  have  a 
pass-band  for  oscillations  between  30  sec  and  5  sec  period,  the  lower 
two  records  for  5  sec  to  1  sec  period.  Note  the  missing  higher  frequencies 
at  College  at  the  first  burst  of  the  event  and  the  similar  onset  times 
for  both  major  activity  pe;.ks. 

Figure  3  illustrates  a  pi  1  event  at  the  Great  Whale-Byrd  conjugate 
stations  on  3  August  1963.  In  this  presentation,  the  pulsation  events 
appear  as  impulsive,  vertical  darkening  of  the  chart.  The  bursts  of 
activity  are  often  simultaneous  but  some  non-conjugate  behavior  is  evident 
even  though  these  stations  are  in  closer  conjugacy  than  College -Macquarie. 

A  station  100  km  away  from  an  E-region  line-current  location  would  sense 
only  about  a  30 %  signal  amplitude  reduction;  exact  conjugacy  of  sites  would 
not  seem  to  be  a  necessary  requirement  in  data  acquisition  for  such 
phenomena.  Difference  in  the  occurrence  of  the  vertical  bands  of  activity 
in  figure  3  should  be  ascribed  to  a  large  but  brief,  changed  location  in 
the  earth's  field  lines  connecting  the  two  sites  or  to  an  asymmetry  of 
the  auroral  particle  precipitation  between  the  two  hemispheres.  Often  the 
spectral  records  show  horizontal  regions  of  enhanced  or  faded  signal  level. 
Since  these  are  usually  band-like,  at  similar  frequencies  for  many  impul¬ 
sive  bursts  at  one  station  but  differing  at  conjugate  stations  (figure  37, 
resonant  transmission  or  selective  absorption  at  ionospheric  heights  is 
probably  responsible. 
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Mather  and  Wescott  (1962)  compared  earth  current  records  for  gi 
events  occurring  at  Cold  Bay  in  the  Aleutipji  Islands  and  Oamaru  in  New 
Zealand.  The  major  disturbances  on  the  two  station  records  occurred 
simultaneously,  had  similar  duration,  and  usually  showed  peak  for  peak 
correspondence.  Other  Alaskan  stations  generally  were  disturbed  at  the 
same  time  but  failed  to  be  correlated  in  the  detailed  variations.  Figure  4 
shows  the  similarity  of  amplitude  changes  at  the  two  sites  for  an  event 
on  April  2,  1962.  Mather  and  Wescott  (1962)  interpreted  the  relative 
loss  of  shorter  period  gi  flucutations  at  Cold  Bay  due  to  the  unique 
ground  conductivity  at  that  site. 

A  series  of  papers  concerning  micropulsation  studies  at  Great  Whale 
River  and  Byrd  Station  have  been  published  by  the  Canadian  research 
groups  at  Pacific  Naval  Laboratory  and  the  University  of  British  Columbia. 
Lokken  et  al.(l96l,  1962 ),  Yanagihara  (1963 ),  Wright  (1964),  Jacobs  (1964), 
Lokken  (1964)  and  Wright  and  Lokken  (1965)  presented  conjugate  data 
obtained  in  January  1961.  Jacobs  and  Wright  (1965)  added  to  the  analysis 
the  conjugate  events  for  the  year  from  February  1963  to  February  1964. 

In  their  comparison  of  Byrd  gi  events  with  those  at  Canadian  stations, 

Great  Whale  River  proved  to  be  best  correlated.  At  times,  though,  it 
seemed  that  a  Byrd  gi  burst  corresponded  best  with  Ft.  Churchill  data. 

The  field  line  location  was  thought  to  flucutate  enough  that  occasionally 
the  two  stations  were  not  at  all  conjugate.  Using  72  selected  events  with 
good  timing,  Jacobs  and  Wright  (1965)  found  two  groups  of  26  and  29  cases 
when  one  station  led  the  other  by  an  average  of  1.1  min  and  one  group  of 
17  cases  with  zero  time  difference. 

In  figure  5,  the  5-nJ.n  average  amplitudes  at  six  world  stations  are 
displayed  for  a  brief  period  of  gi  activity  in  November  1963.  The  distur¬ 
bance  starting  about  1800  UT  on  7  November  probably  reached  from  Sweden 
to  Alaska  across  Siberia.  Northern  Canada  seemed  to  catch  most  of  the 
activity  at  0200  on  November  8.  The  pulsations  starting  near  0520  UT 
were  spread  over  the  western  hemisphere  and  extended  to  low  latitude. 

The  conjugate  stations  of  Great  Whale  River-Byrd  exhibited  a  correlation 
coefficient  of  0.7  for  5-nJ.n  scalings  of  the  data  during  the  period  of 
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figure  5;  other  pairs  of  stations  gave  much  lower  coefficients.  In  one 
hemisphere,  great  differences  in  thi'  amplitude  changes  for  pi  typically 
occur  at  distances  of  1000  km  and  disappear  for  200  km  separation.  The 
conjugate  region  may  be  within  such  dimensions.  The  distribution  of 
activity  for  pi  events  in  figure  5  in  quite  similar  to  the  regions  of 
ionospheric  absorption  of  cosmic  noise  with  which  the  pi  are  clearly 
related  (CampbeU  and  Leinbach,  1961). 

Equatorial  locations  sometimes  have  been  called  "self  conjugate". 

In  a  study  of  the  long-period  irregular  fluctuations  which  are  most 
numerous  during  daytime  hours,  Onwumechilli  and  Ogbuehi  (1962)  reported 
that  about  9Cfjo  of  the  events  at  Ibadan,  Nigeria,  could  be  identified  at 
other  observatories  along  the  equator  and  within  the  auroral  zone. 

Wescott  and  Mather's  (1965c)  results  for  long  period  variations  at  San 
Juan,  Puerto,  Rico,  and  Trelew,  Argentina,  near  conjugate  locations  with 
L  of  about  l.U,  were  in  agreement  with  the  Nigerian  study.  The  close 
association  of  the  amplitude  of  the  fluctuations  with  the  strength  of 
the  equatorial  Sq  current  system  implies  that  the  field  variations  are 
actually  changes  in  Sq  communicated  to  the  equatorial  regions  from  the 
higher  latitudes. 

The  longer-period  irregular  variations  were  compared  for  College  and 
Little  America  by  Ahmed  and  Scott  (1955)-  Activity  at  conjugate  locations 
in  middle  and  auroral  latitudes  were  more  recently  reported  in  a  series 
of  papers  by  Wescott  (1961,  1962),  Wescott  and  Mather  (1963,  1965a,  b) 
and  Brown,  Barcus,  and  Parsons  (1965).  Between  Cold  Bay,  Alaska,  and 
Oamaru,  New  Zealand,  activity  correlation  was  about  the  same  or  slightly 
poorer  than  auroral  zone  pairs.  Macquarie  Island  records  compared  best 
with  Alaskan  stations  at  Healy  and  Kotzebue.  The  region  of  coherence  of 
signals  in  Alaska  was  elongated  along  the  line  of  constant  L  about  350  km 
E-W  and  120  km  N-S.  There  was  a  seasonal  change  in  the  relative  amplitude 
of  the  field  disturbance  variations  such  that  the  southern  station  was 
more  active  in  the  December  solstice  months  and  the  northern  station  more 
active  in  the  June  solstice  months.  It  appears  that  ionospheric  conduc¬ 
tivity  and  currents  play  an  important  role  in  the  conjugate  occurrence 
region  and  amplitudes  for  these  auroral  zone  field  variations. 
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Irregular  field  fluctuations  seen  on  magnetometers  at  Shepherd  Bay, 
Canada,  and  Scott  Base,  Antarctica,  both  within  the  auroral  zone,  exhibited 
their  greatest  similarity  during  the  nighttime  hours,  a  secondary  maximum 
of  moderate  correlation  about  midday,  and  usually  dissimilar  behavior 
between  the  two  maxima  (Wescott  and  Mather,  1963).  One  interpretation 
of  figure  6  is  that  during  the  midday  hours  Sq-like  currents  carry  the 
variations  in  from  the  auroral  zones  whereas  at  nighttime  either  the  high 
latitude  field  lines  are  connected  and/or  the  auroral  magnetic  distur¬ 
bance,  associated  with  particle  precipitation,  is  spread  over  a  broader* 
region. 

Magnetic  storm  sudden  commencements  (SSC)  and  magnetic  field,  sudden 
impulses  (Si)  are  typically  associated  with  pi  type  pulsations  in  the 
auroral  zone  (Campbell  and  Matsushita,  1962).  More  than  50$  of  SSC's 
are  associated  with  pulsations  of  2.5  to  10  minutes  period  (Sugiura  and 
Wilson,  1961).  For  that  reason  some  results  of  three  studies  of  SSC  and 
SI  phenomena  are  mentioned  here. 

The  SSC's  are  rapid  changes  in  the  geomagnetic  field  which  last  about 
1  to  6  minutes,  breaking  a  "quiet  period  of  recording,  and  occur  at  the 
beginning  of  a  magnetic  storm.  Sugiura  and  Wilson  (1961)  and  Nagata  et 
al.  (1966)  found  that  conjugately  located  auroral  stations  have  SSC's  with 
similar  appearance,  occurrence  frequency  and  polarization  vectors  showing 
elliptic  countci clockwise  rotation  (viewed  along  the  field  line)  from 
22  h  to  10  h  local  geomagnetic  time.  In  this  feature  the  SSC  behaves  like 
pc  5  pulsations.  The  relative  sizes  of  the  SSC's  at  Reykjavik,  Iceland, 
and  Syowa,  Antarctica,  seemed  to  indicate  slightly  larger  amplitudes  in 
the  winter  hemisphere.  Sato  (1962)  gave  excellent  examples  of  the  concur¬ 
rent  ionospheric  absorption  effects  associated  with  the  SSC  simultaneously 
occurring  in  Alaska  and  Antarctica. 

The  Si's  are  sudden,  world-wide  changes  in  the  magnetic  field,  much 
like  the  SSC's  but  not  followed  by  a  magnetic  storm.  The  general  conjugate 
features  were  the  same  as  that  reported  above  for  the  SSC. 
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Longer  period  continous  pulsations  of  the  £c  4-5  variety  are  often 
seen  within  the  Si.  2  spectra.  Figure  7  illustrates  that  type  of  conjugate 
pi  micropulsation  received  on  ITSA  recorders  at  Great  Whale  River  and  Byrd 
Station.  Three  aspects  of  this  figure  should  be  pointed  out.  There  is  a 
good  agreement  in  the  occurrence  time  for  most  of  these  £i  bursts.  Saito 
(1967)  used  an  interesting  photographic  technique  for  displaying  this 
coincident  spectral  pattern.  A  higher  frequency  portion  of  activity 
between  08  h  and  09  h  UT  is  missing  at  Great  Whale  River  whereas  a  higher 
frequency  portion  between  09  and  10  UT  is  absent  at  Byrd.  There  is  a  well 
developed  pc  4-5  in  the  early  hours  of  the  event  that  becomes  the  major 
micropulsation  activity  after  10  h  UT. 

There  are  a  few  reports  of  the  behavior  of  the  long  period,  damped 
oscillation  type  pulsations  (jgt)  at  conjugate  locales.  Some  aspects  of 
these  may  place  them  in  the  transitory  region  between  pi  pc  pulsation 
types  or  miniature  version  of  the  SI  mentioned  earlier.  Saito’ s  (1961) 
hodographs  of  the  pt  disturbance  vector  in  the  horizontal  plane  showed  a 
northern  hemisphere  convergence  into  a  region  lying  along  the  auroral  zone 
near  the  local  midnight  meridian  while  vectors  in  the  southern  hemisphere 
diverged  from  the  approximate  conjugate  point.  Jacobs  and  Wright  (1965) 
show  one  example  of  a  large  4  to  5  min  period  £t  .  observed  simultaneously 
at  Great  Whale  River  and  Byrd  (figure  8). 

3.  Continuous  Pulsations  with  Periods  of  One  to  Ten  Minutes 

The  long  period  continuous  pulsations,  pc  4-5  generally  have  large 
amplitudes  reaching,  at  times,  to  several  hundred  gammas.  They  occur 
principally  in  the  auroral  latitudes  in  periods  of  high  solar-terrestrial 
activity.  The  forms  are  quite  sinusoidal,  often  damped,  and  the  duration 
ranges  from  a  few  minutes  to  several  hours.  Many  events  have  harmonic 
structure.  The  spectra  in  figure  7  from  1015  to  1200  h  UT  is  typical  of 
the  detailed,  fine  structure  coincidence  found  at  conjugately  located  high 
latitude  stations. 

The  reports  by  Sugiura  (1961)  on  the  conjugate  College-Macquarie  pair 
of  stations,  by  Nagata  et  al.  (1963,  1966)  and  Kckubun  and  Nagata  (1965) 
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on  the  Reykj avik-Syowa  stations  pointed  to  the  correspondence  of  such 
pulsations  at  conjugate  locations  and  their  hydromagnetic -wave -like 
behavior.  The  signals  seemed  not  unlike  the  magnetic  impulses  of  SCMS 
(Wilson  and  Sugiura,  1961)  with  the  disturbance  vector  rotation  sense 
invariant  along  the  whole  line  of  force  in  space;  a  counterclockwise 
direction  (fig.  9)  in  the  horizontal  field  vector,  H,  during  the  morning 
(looking  into  the  northern  hemisphere)  and  reversed  in  the  afternoon  with 
a  transition  about  11  h  local  time.  The  main  axis  of  the  polarization  was 
almost  perpendicular  to  the  direction  of  the  main  geomagnetic  force. 
Correlation  coefficients  between  successive  maxima  and  minima  were  about 
0.8.  Amplitudes  were  close  to  the  same  value  at  the  two  stations  but 
amplitudes  of  the  morning  pulsations  were  appreciably  larger  than  the 
afternoon  variety.  Frequency  spectra  show  usually  one  but  often  two 
discrete  peaks.  The  pulsations  were  ascribed  to  hydromagnetic  waves, 
resonating  with  the  eigen-oscillation  of  the  geomagnetic  field  line  and 
excited  by  the  arrival  of  the  solar  wind. 

Ahmed  and  Scott  (1955)  found  gc  4-5  with  similar  occurrence  and  period 
at  College  and  Little  America.  We scott  and  Mather  (1965b)  mention  £c  4-5 
events  observed  at  the  Cold  Bay-Oamaru  station  pair.  Jacobs  and  Wright 
(1965)  and  Saito  (1967)  published  some  similar  £c  4-5  events  at  the  Great 
Whale-Byrd  station  pair. 

Recently  Wilson  (1966)  reexamined  some  of  the  College-Mac quarie  data 
and  concluded  that  ionospheric  current  systems  were  induced  by  the  circu¬ 
larly  polarized  standing  hydromagnetic  waves  on  the  field  lines.  He 
ascribed  differences  in  signal  appearance  at  near  conjugate  locations  to 
the  position  of  the  observing  station  relaJve  to  the  terminal  point  of 
the  oscillating  field  line.  His  model  seemed  to  leave  unexplained  the 
interrelationship  of  the  j>c  4-5  with  ionospheric  absorption  (Sato,  1964, 
1965),  auroral  radar  reflection  (Kaneda  et  al.,  1964)  and  precipitating 
auroral  electrons  (Anger  et  al.,  1963,  and  Barcus  and  Rosenberg,  1965). 

Sato's  (1964,  1965)  study  used  the  data  from  a  world  network  of 
stations.  He  found  most  pulsations  in  Antarctica  to  be  a  local  phenomenon 
with  occurrence  maxima  different  in  each  area  and  phases  not  coincident 
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between  stations.  Figure  10  illustrates  the  amplitude  and  phase  similarity 
of  the  pulsations  he  studied  at  the  conjugate  College -Macquarie  pair  of 
stations.  The  correspondence  with  concurrent  riometer  measured  ionospheric 
absorption  led  Sato  to  the  conclusion  that  most  high  latitude  pulsations 
cf  this  type  are  caused  by  currents  arising  from  the  high  latitude,  particle 
bombardment  of  the  ionosphere. 

4.  Pulsations  with  5  to  50  sec  Period 

Mather  et  al.  (1964)  studied  an  active  period  in  April  1962  during 
■which  pulsations  with  about  40  sec  period  occurred  at  Cold  Bay  and  Oamaru. 
The  signal’s  behavior  is  rather  like  that  expected  of  ££  4  type  pulsations. 
The  amplitude  and  the  periods  were  somewhat  shorter  during  the  daytime. 

There  was  some  tendency  to  slightly  longer  periods  at  the  time  of  auroral 
disturbance  onset.  Figure  11  is  an  example  of  their  study  on  the  direction 
of  rotation  of  the  disturbance  vector  for  the  two  stations.  Four  of  their 
ten  simultaneous  conjugate  disturbances  maintained  exactly  the  same  phase 
relationship  with  oppositely  rotating  vectors.  The  polarization  sense  was 
uniform  at  night,  counterclockwise  (looking  into  the  ground)  at  Cold  Bay 
and  clockwise  at  Oamsiru.  The  rotation  pattern  was  more  irregular  during 
the  day  but  generally  reversed  to  the  nighttime  behavior.  The  earth-sun 
line  division  of  the  sense  change  seemed  at  variance  with  the  results  cf 
Wilson  and  Sugiura  (1961)  and  Nagata  et  al.  (1963)  for  the  longer  period 
variations. 

-  ~  true  2-3  pulsations  are  typically  a  daytime  phenomena  with 

longer  periods  at  midday  and  amplitudes  proportional  to  the  Kp  index  of 
world  wide  geomagnetic  disturbance.  The  signals  usually  stay  under  1  gamma 
in  size  and  maximize  in  the  auroral  region.  Examples  of  simultaneous  high 
latitude  appearance  of  these  pulsations  have  been  published  by  Lokken  et  al. 
(1962)  for  Victoria  and  Byrd,  by  Campbell  (1963)  for  College  and  Macquarie 
(fig.  ll),  and  by  Jacobs  and  Wright  (1965)  for  Great  Whale  and  Byrd. 

At  low  latitudes  Kate  and  Okuda  (l95o)  reported  similar  30  to  50  sec 
pulsations  at  Japan  and  Ceylon.  Schlumberger  and  Kunetz  (1946)  reported 
a  possible  reversal  of  signal  polarization  between  France  and  Madagascar 
for  these  signals. 
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The  £c  2-3  are  the  most  poorly  studied  pc  group  of  pulsations.  The 
polarization  sense  and  major  direction  of  the  horizontal  disturbance  vector 
undergoes  a  regular  diurnal  variation  at  most  latitudes  (Terada,  1917; 
Kawamura  et  al.,  1961)  which  can  account  for  some  of  the  hemispherical 
differences  reported.  The  peak  for  peak  identification  at  conjugate 
stations  does  exclude  the  local  ionospheric  resonance  as  a  generative 
model . 


5.  Regular  Pulsations  in  the  1  sec  Period  Range 

-  The  pc  1  pulsations  occur  principally  in  the  band  from  0.15  to  3.5 
c/s.  They  usually  have  a  rising  frequency  structure  but  variations  have 
been  observed.  The  average  frequency  and  emission-recurrence  period  seem 
to  be  linearly  related.  Bursts  of  activity  typically  last  about  one-half 
hour;  the  general  disturbance  level  may  be  as  long  as  twelve  hours.  The 
averaged  diurnal  variation  in  mid-frequencies  is  similar  at  all  world 
locations:  minima  of  about  .4  or  .5  c/s  occur  from  l400  to  1800  hours 
local  time  and  maxima  of  about  .9  or  1.0  c/s  occur  from  0100  to  0500  hours 
local  time.  The  occurrences  are  usually  near  the  local  noon  hours  at 
auroral  latitude  stations  and  they  shift  into  the  morning  to  predawn  hours 
at  middle  and  low  latitudes. 

At  times  the  events  are  found  simultaneously  at  widely  spaced  stations 
over  a  hemisphere  of  the  earth  (Tepley  et  al.,  1965;  Kenney  and  Knaflich, 
1967).  The  usual  case,  however,  is  for  a  rather  local  appearance  of  the 
events.  This  is  evidenced  in  the  unique  local  occurrence  maxima  as  a 
function  of  latitude  and  in  the  local-time  diurnal  variation  of  mid¬ 
frequency  for  the  events.  Campbell  and  Stiltner  (1965)  and  Campbell  (1967a) 
demonstrated  the  conjugate  station  similarity  of  these  two  characteristics. 
They  placed  the  typical  range  for  coherent  signal  to  be  under  1000  km  in 
the  high  latitudes.  In  contrast,  Borsoukov  and  Ponsot  (1964)  found  phase 
coherence  at  1100  km.  Campbell  and  Stiltner  (1965)  only  found  significant 
correlation  coefficients  for  data  from  stations  conjugately  located. 

Figure  13  shews  the  cross-correlation  for  joint  sampling  periods  at  College- 
Macquarie,  3reat  Whale -Byrd,  Bale  St.  Paul  (Quebec) -Eights.  The  highest 
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correlation  was  obtained  from  this  last  pair.  These  stations  were  not 
only  nearest  to  computed  conjugacy  but  closer  to  each  other  in  local  time 
and  at  lower  L  values  than  the  other  two  pairs. 

Figure  l4  is  an  example  of  the  spectra  of  ]3c  1  everts  received  at 
conjugate  locations.  In  addition  to  the  last  two  references  a  number  of 
authors  have  published  examples  of  conjugate  £c  1  signals.  Dawson  (1965 ) 
discussed  occurrence  and  similar  average  frequency  changes  at  College  and 
Macquarie.  Lokken  et  al.  (l963)>  Yanagihara  (1963),  51111  Wright  and  Lokken 
(1965)  presented  examples  of  concurrent  events  at  the  Great  Whale-Byrd 
pair  of  stations.  The  data  received  at  Kerguelen  and  Sogra  or  Borok 
appeared  in  publications  of  Borsoukov  and Ponsot  (1964),  Laurent  et  al. 

(1964),  Troitskaya  et  al.  (1964),  Gendrin  and  Troitskaya  (1965),  Lacourly 
et  al.  ( 1965 ) »  Loran  et  al.  (1965),  Matveyeva  and  Troitskaya  (1965)  and 
Ger.drin  et  al.  (1966).  Schlicht  (1963)  reported  some  similarities  in  the 
records  taken  at  Chambon-la-Foret  and  Garchy,  France,  and  at 
Ksrnanus,  South  Africa.  The  pc  1  at  the  near  conjugate  locations  of  Kauai 
and  Tcngatapu  have  been  presented  by  Tepley  (1964,  1965,  1966). 

A  number  of  unique  features  have  come  to  light  in  the  conjugate  studies 
which  should  give  insight  into  the  magnetospheric  propagation  character¬ 
istics  and  separate  out  the  local  ionospheric  effects  on  the  signal  recorded. 
At  conjugate  stations  with  greatly  differing  geographic  longitude,  the 
average  occurrence  and  diurnal  changes  in  mid -frequency  are  best  matched 
at  the  two  sites  when  a  common  field-line  time  is  assumed  (Campbell  and 
Stiltner,  196?).  Some  of  the  gross  appearance  of  the  signal  is  immediately 
evident  in  viewing  such  examples  as  figure  1- .  There  is  often  a  general 
activity  recurrence  pattern  of  about  one-half  cr  one  hour.  The  general 
envelope  of  the  event  is  often  similar  at  the  two  stations,  for  example 
in  the  early  hours  of  figure  14.  These  features  and  the  horizontal 
"streaking"  cf  signal  enhancement  and  disappearance  must  not  be  dependent 
upon  the  transient  ionospheric  conditions  which  vary  independently  in 
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a  £c  1  events.  However,  as  the  figure  shows,  both  senses  are  observed. 
Theoretical  models  seem  to  predict  a  left-hand  polarized  (looking  into 
the  northern  hemisphere),  Alfven-mode  of  field  guided  propagation  in  the 
magnetosphere.  The  results  of  this  figure  must  modify  that  assumption 
or  introduce  increased  importance  to  propagation  processes  within  the 
ionosphere . 

The  interlacing  of  £c  1  fine  structure  has  been  demonstrated  by 
several  independent  researchers.  Figure  15  also  illustrates  this  alter¬ 
nation  between  the  two  hemispheres.  Lokken  et  al.  (1963 )  and  Yanagihara 
(1963),  examining  Great  Whale  River  and  Byrd  Station  records,  found  a  few 
remarkable  cases  when  the  amplitude  envelope  of  activity  alternated  in 
the  northern  and  southern  auroral  zones  with  the  maxima  of  activity  at 
one  station  occurring  midway  between  the  maxima  at  the  other.  Tepley's 
(1964)  studies  in  the  Pacific  area  indicated  that  emissions  observed 
simultaneously  in  the  same  hemisphere  were  simultaneous  within  an  accuracy 
of  ±  12  seconds  whereas  emissions  observed  at  Kauai  and  Tongatapu,  at  near 
conjugate  locations,  had  structural  elements  displaced  in  time  by  one-half 
the  fine  structure  periodicity.  He  also  observed  occasional  doubling  of 
the  fine  structure  repetition  frequency  at  middle  and  low  latitude  stations. 
Borscukov  and  Ponsot  (1964)  calculated  correlation  coefficients  and  auto¬ 
correlation  functions  for  simultaneous  £c  1  events.  They  found  no  time 
differences  greater  than  0.1  sec  between  Borok  and  Lovozero,  a  distance 
of  about  1100  km;  the  oscillations  were  in  phase.  The  authors  found  no 
phase  coherence  between  the  conjugate  regions  of  Borok  and  Kerguelen  for 
the  individual  oscillations.  Figure  16  shows  the  correlation  function 
for  various  time  shifts  between  the  envelopes  of  activity  received  alter¬ 
nately  in  the  two  hemispheres.  The  figure  indicates  that  the  signals 
appeared  at  Kerguelen  about  70  sec  after  its  appearance  at  Borok  whereas 
the  repetition  period  was  120  sec.  Gendrin  and  Troitskaya  (1965)  point 
out  that  this  discrepancy  may  be  due  to  some  heating  between  overlapping 
emissions.  Wentworth  et  al.  (1966)  reported  that  for  Pacific  stations  in 
the  northern  hemisphere  the  gc  1  seem  to  penetrate  the  ionosphere  at  high 
latitudes  and  propagate  equatorward  with  an  average  apparent  velocity  of 
about  900  km/sec.  They  interpreted  such  propagation  time  as  responsible 
for  the  difference  from  exact  l80°  phase  shift  measured  for  signals 
received  in  opposite  hemispheres. 
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Apparent  attenuation  of  the  ionosphere  upon  pc  1  signals  was  investi¬ 
gated  by  Campbell  (1967).  Studying  the  0.5  to  1.0  c/s  pc  1  emissions  at 
the  conjugately  located  stations  he  selected  Northern  summer  solstice  times 
when  the  ionosphere  was  dark  over  Eights,  Antarctica,  but  cycle  normally 
over  Bale  St.  Paul,  Canada.  Table  2  indicates  that  the  signal  level  was 
reduced  under  the  daytime  ionosphere.  At  the  higher  frequencies  and  with 
more  disturbed  (higher  Kp)  ionospheric  conditions  greater  attenuation  was 
found.  Slight  changes  in  seasonal  occurrence  and  signal  frequency  were 
shown  to  be  consistent  with  the  estimated  attenuation.  The  signal  ampli¬ 
tude  reduction  was  matched  to  theoretical  t=::pectations  by  the  assumption 
of  Alfven  to  modified-Alfven  mode  coupling  of  the  propagated  signal  in  the 
ionosphere. 

Several  features  of  the  pc  1  now  seem  quite  evident.  The  phenomenon 
is  of  magnetospheric  origin.  It  is  field  guided  into  high  latitude 
stations  without  appreciable  ionospheric  attenuation.  The  signals  may 
travel  to  lower  latitudes  in  other  ducted  modes.  There  are  still  numerous 
unexplained  characteristics  of  the  form  of  these  pc  1  events  which  cannot 
be  ascribed  to  local  anomalies. 
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♦Computations  for  November  1966  and  100  km  altitude  above  station  (Campbell  and  Matsushita,  1967). 


Table  2 


Average  Amplitude  Ratio  (Quebec/Eights)  for  Dark  Ionosphere 
at  Eights  and  Daylight  at  Quebec  (Campbell,  1967). 


Frequency  c/s 

Kp  index 

0.5 

1.0 

0  to  1+ 

.70 

•  49 

2-  to  2+ 

•‘39 

.45 

3-  to  3+ 

•  53 

.19 
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Figure  3*  Sonagraph  spectral  display  of  conjugate  behavior  of  gi  signals  at  Great 
Whale  River  and  Byrd  on  August  3,  1963  (Campbell,  1965 ). 


Figure  U.  Telluric  current  charts  for  April  7,  1962,  at  Cold  Say,  Alaska,  and  C  maru, 
New  Zealand  (Mather  and  Weacott,  1962} . 


PERIOD  (sec) 


Figure  9-  Example  of  locua  of  rotating  vectors  for  the  morning  pulsations  in  the 
horltontal  plane  (top)  and  in  the  magnetic  plane  (bottom)  observed  at 
Reykjavik  (Re)  and  Syowa  (Sy)  at  0655  to  0705  UT  on  Setpember  18,  1959 
(Hagata  et  al. ,  1963). 


Figure  10.  Slailtaneoua  £c  5  magnetic  field  pulsation,  at  College  and  Macquarie 
Island  on  19  March  1958.  Concurrent,  ionospheric  absorption  measured 
with  the  College  Rlameter  is  shown  on  a  similar  time  scalfe  at  the  top 
(Sato,  1965). 
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Figure  II.  Polarization  of  40  sec  pulsations  received  at  Cold  Bay  and  Oaaaru  for 
a  24  hr  period  on  11-12  April  1967.  Radial  arrows  drawn  from  the 
reference  circle  denote  the  direction  of  rotation  at  various  times 
throughout  the  day.  Long  irrovo  directed  outward  denote  clockwise 
rotation;  directed  inward,  counterclockwise.  Short  bi-directional 
arrows  mean  both  sense  or  linear  polarization  were  present  with  no 
preferred  direction  (Mather  et  al.,  1964). 
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Figure  12.  A  pc  2-3  event  recorded  simultaneously  along  the  Pacific  Ocean  region  on 
12  February  1962  (Campbell,  1963). 
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Figure  13.  Monthly  cross -correlation  diagrams  for  5-mln  scalings  of  pc  1  events  at 
three  sets  of  conjugate  stations  for  lead  and  lag  period  to  two  hours 
(Campbell  and  Stiltner,  1965). 
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Figure  15.  Average  polarization  sense  in  horizontal  plane  for  ££.  1  events  measured  simultaneously  for  the 
periods  indicated  at  Sogra,  USSR,  and  Kerguelen  Islands  on  U  March  1965.  Percentage.-  indicate 
the  proportional  mixture  of  +  and  -  rotation  sense  (Oendrin  et  al. ,  1966). 


Figure  16.  Correlation  function  for  envelopes  of  pc  1  signals  received 
at  Borok  and  Kerguelen  on  February  17 ,  19&  (Borsukov  and 
Ponsot,  196^)  ■ 
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Pi  Events 


J.A.  Jacobs 
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Vancouver  8,  B.C.,  Canada 


1.  Introduction 

The  amount  of  literature  on  irregular  micropulsations 
(Pi)  is  relatively  sparse,  particularly  when  compared  to  that 
on  the  more  regular,  continuous  pulsations,  and  very  little 
work  has  been  done  on  the  characteristics  of  Pi  activity  at 
conjugate  points.  Originally  such  irregular  micropulsations, 
which  tend  to  appear  during  the  initial  stages  of  relatively 
large  scale  geomagnetic  disturbances  known  as  bays,  were 
designated  Pt.  At  the  13th  General  Assembly  of  the  IUGG, 
Jacobs  et  alia  (1964)  proposed  the  present  classification  of 
impulsive  micropulsation  activity  known  as  Pi;  this  was  sub¬ 
divided  into  the  categories  Pi  1  (periods  1-40  sec.)  and  Pi  2 
(periods  40-150  sec.).  Pi  1  are  usually  seen  in  the  initial 
stages  of  a  Pi  2  and  appear  as  riders  on  the  longer  period 
pulsations . 
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2.  PI  2 

2.1  The  frequency  spectrum  of  PI  2's 

As  In  the  case  of  auroral  and  geomagnetic  bay  acti¬ 
vity,  Pi  2  micropulsations  have  a  maximum  occurrence  frequency 
around  local  midnight.  The  relationship  between  the  occur¬ 
rence  of  Pi  2  and  planetary  magnetic  activity  is,  however, 
not  clear.  Yanagihara  (1959)  has  presented  evidence  that, 
while  auroral  and  geomagnetic  bay  activity  increase  towards 
sunspot  maximum.  Pi  2  activity  in  low  latitude  regions 
decreases.  On  the  other  hand,  he  found  that  the  probability 
of  occurrence  of  a  Pi  2  increases  with  increasing  Kp  index 
(Yanagihara,  I960).  At  the  same  time  Pi  1  activity  increases 
towards  sunspot  maximum,  suggesting  that  either  Pi  1  and  Pi  2 
have  different  source  mechanisms  or  that  they  have  different 
modes  of  propagation  over  the  earth's  surface.  Rostoker 
(1967a)  has  shown  that  Pi  2  may  be  considered  to  consist  of 
two  separate  contributions:  the  initial  kick  which  is  a  broad 
band  effect  and  the  actual  damped  oscillation.  Sometimes 
the  initial  kick  may  occur  but  insufficient  energy  is  imparted 
to  the  mechanism  of  generation  to  cause  an  oscillation  to  take 
place.  Such  kicks  appear  often  prior  to  an  auroral  breakup, 
and  Rostoker  has  referred  to  these  occurrences  as  "one  pulse 
effects" . 
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Spectral  analysis  of  Pi  2  trains  bear  out  this 
picture.  Rostoker  (1967b)  found  that,  exclusive  of  Pi  1 
activity,  a  Pi  2  may  contain  anywhere  from  one  to  several 
frequency  components  (see  Pigs.  1  and  2).  None  of  the 
frequency  components  in  the  spectra  he  obtained  were  present 
prior  to  the  initiation  of  the  Pi  2  event.  In  his  analysis 
of  18  Pi  2  events  periods  as  long  as  286  seconds  were  found. 
Thus  it  appears  as  though  the  range  of  Pi  2  should  perhaps  be 
extended  to  300  sec. 


2.2  The  relationship  between  Pi  2  activity  and  the  Kp  index 

Rostoker  also  found  that  the  number  of  frequency 
components  in  a  Pi  2  event  increased  approximately  linearly 
with  Increasing  Kp  index  (see  Pig.  3)»  reminiscent  of  the 
relationship  obtained  by  Yanagihara  (I960)  of  the  probability 
of  occurrence  of  a  Pi  2  as  a  function  of  the  Kp  index.  The 
frequencies  measured  in  each  event  were  usually  not  harmonics 
of  a  fundamental  frequency  and  thus  each  frequency  may  be 
thought  of  as  having  a  separate  generation  from  one  another. 

Troitskaya  (1964)  has  given  evidence  that  the  period 
of  a  Pi  2  decreases  with  increasing  Kp  index;  however,  this 
conclusion  is  a  result  of  the  assignment  of  a  single  frequency 
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to  a  Pi  2  event  by  visual  analysis  of  chart  recordings. 

With  the  aid  of  power  spectra  of  18  Pi  2  events,  Rostoker 
(1967b)  has  shown  that  the  minimum  period  in  a  Pi  2  event 
decreases  with  increasing  Kp  while  the  maximum  period  increases 
with  increasing  Kp  up  to  Kp  2o,  after  which  it  decreases 
slightly  with  further  increase  in  Kp.  Visual  analysis  of  the 
18  events  reveals  a  decrease  in  dominant  apparent  period  with 
Increasing  Kp,  in  accordance  with  the  results  of  Troitskaya. 
Rostoker' s  results  also  indicate  that  the  period  of  a  given 
Pi  2  event  does  not  change  with  latitude  for  a  single  event, 
although  Eleman  (1966)  has  shown  that  the  average  period  of 
many  Pi  2's  at  a  given  station  is  longer  the  closer  the 
observing  station  is  to  the  auroral  zone. 


2 . 3  Polarization  characteristics  of  Pi  2's 

It  is  generally  agreed  that  Pi  2's  are  polarized  in 
the  horizontal  plane  in  mid  and  low  latitude  regions,  although 
there  is  usually  a  relatively  strong  vertical  component  in  the 
vicinity  of  the  auroral  zone  of  an  event.  The  results 
obtained  at  low  latitude  stations  by  Kato  et  alia  (1956)  indi¬ 
cate  that  Pi  2's  are  polarized  in  a  clockwise  sense  before 
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local  midnight.  This  is  in  disagreement  with  the  results  of 
Christoffel  and  Linford  (1966)  who  found  that  Pi  2*s  were 
polarized  predominantly,  although  not  always,  in  a  clockwise 
sense  in  midlatitude  regions  in  the  southern  hemisphere. 

Rostoker  (1967a)  showed  that  the  polarization  of 
Pi  2's  in  midlatitude  regions  in  the  northern  hemisphere  is 
predominantly  in  the  anticlockwise  sense,  although  cases  of 
clockwise  rotation  are  observed  on  occasion.  He  found  that 
the  sense  of  polarization  does  not  depend  on  the  position  of 
the  observing  station  relative  to  local  midnight,  but  that 
its  geographical  location  '-'■'ative  to  the  source  of  the  Pi  2 
directly  influences  the  polarization  characteristics. 

The  implication  of  the  results  of  Rostoker  and 
those  of  Christoffel  and  Linford  is  that  at  least  the  last  part 
of  the  propagation  of  Pi  2’s  may  be  dominated  by  the  Alfven 
mode  of  a  hydromagnetic  wave  propagating  through  the  iono¬ 
sphere  in  a  direction  parallel  to  the  ambient  field  lines. 
However,  the  fact  that  some  clockwise  pulsations  are  observed 
in  the  northern  hemisphere  and  some  anticlockwise  pulsations 
in  the  southern  hemisphere  suggests  that  the  mode  of  propaga¬ 
tion  of  a  Pi  2  is  not  a  pure  Alfven  wave. 

The  observed  polarization  characteristics  indicate 
that  Pi  2  micropulsations  cannot  be  described  by  a  simple 


equivalent  current  system  such  as  that  used  to  describe  a  bay. 
Such  a  current  system  would  provide  an  "in  phase"  relationship 
between  the  H  and  D  components  which  is  not  observed  to  be 
the  case.  Nor  can  a  rotating  line  current  such  as  that  pro¬ 
posed  by  Wilson  (1966)  for  Pc  4,  Pc  5,  and  Pg-type  micropul¬ 
sations  fully  describe  the  source  mechanism  for  Pi  2’s,  since 
for  anticlockwise  rotation  of  the  current  system  all  pulsa¬ 
tions  in  the  northern  hemisphere  should  be  polarized  in  an 
anticlockwise  sense,  which  is  not  the  case.  The  appearance 
on  occasion  of  a  clockwise  rotation  at  a  station  while  the 
sense  of  polarization  at  other  stations  for  the  same  event 
is  anticlockwise,  together  with  the  fact  that  the  dominant 
sense  of  rotation  is  anticlockwise  in  the  northern  hemisphere, 
and  clockwise  in  the  southern  hemisphere,  suggests  that  the 
disturbance  is  subject  to  propagation  effects  and  hence  in 
midlatitude  regions  Pi  2’s  may  be  considered  to  propagate  as 
a  wave  rather  than  being  the  effect  of  an  ionospheric  current 
system. 


Although  Pi  2  micropulsations  as  a  whole  seem  to 
be  unrelated  to  the  bay  type  equivalent  current  system,  the 
initial  kick  on  the  magnetogram  appears  to  be  strongly 
correlated  with  the  polarity  of  the  accompanying  bay.  In 
general  the  direction  of  the  initial  kick  of  either  component 
on  the  micropulsation  record  has  the  same  sign  as  the  polarity 


111-2-6 


of  the  corresponding  component  of  the  bay;  this  effect  is 
strongest  in  the  D  component,  although  it  is  also  prominent 
in  the  H  component. 

It  is  suggested  that  the  initial  kick  of  a  Pi  2 
micropulsation  may  be  represented  by  an  equivalent  overhead 
current  system  in  the  manner  shown  by  Jacobs  and  Sinno  (I960). 
The  current  is  set  up  in  the  Hall  region  of  the  ionosphere  as 
a  consequence  of  the  fact  that  the  ambient  ionospheric  parti¬ 
cles  cannot  respond  sufficiently  rapidly  to  the  first  hydro- 
magnetic  impulse  generated  in  the  breakup. 


2 .  '  r’-,-:jible  source  mechanisms  for  Pi  2  activity 

The  initial  hydromagnetic  impulse  arriving  in  the 
auroral  zone  from  the  tail  region  of  the  magnetosphere  is 
broad  band,  containing  periods  ranging  from  those  of  a  Pi  to 
a  bay.  Near  the  source  all  the  frequency  components  arrive 
at  approximately  the  same  time,  and  thus  the  leading  edge  of 
the  bay  near  the  source  is  very  steep.  However  the  disper¬ 
sive  effect  of  the  ionosphere  on  the  propagation  to  lower 
latitudes  causes  the  higher  frequencies  to  arrive  at  a  field 
point  earlier  than  the  lower  frequencies;  this  effect  results 
in  a  spreading  on  the  time  axis  of  the  Impulse,  and  the  lead- 
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ing  edge  of  the  bay  becomes  less  steep  the  further  it  is 
observed  from  the  source.  In  this  respect  the  bay  and 
initial  pulse  of  the  Pi  2  may  be  regarded  as  part  of  the 
response  of  the  ionosphere  to  a  hydromagnetic  impulse  gener¬ 
ated  in  the  magnetosphere.  Ward  (1963)  has  endeavored  to 
find  the  transfer  function  of  the  magnetosphere,  that  is, 
how  the  incoming  solar  plasma  generates  hydromagnetic  dis¬ 
turbances  which,  when  modified  by  the  ionosphere,  are  received 
on  the  surface  of  the  earth. 

Rostoker  (1967a)  suggests  that  Pi  2's  originate  as 
eigen  oscillations  on  auroral  field  lines,  although  the 

initial  movement  is  probably  due  to  a  current  system  set  up 
in  the  ionosphere  as  an  impulse  response  tc  the  Influx  of 
energy  from  the  tail  of  the  magnetosphere.  That  part  of  the 
Pi  2  generated  by  the  eigen  oscillations  of  the  field  lines 
is  thought  to  propagate  through  the  E  region  of  the  ionosphere 
and  the  observed  polarization  at  a  field  point  is  determined 
by  the  screening  effect  of  the  ionosphere  over  that  point. 

Herron  (1966)  observed  sizeable  phase  shifts  (of 
the  order  of  seconds)  in  the  appearance  of  micropulsations  at 
stations  separated  by  hundreds  of  kms  along  an  east-west 
line.  From  Herron's  results  the  relationship  between  the 
projected  velocity  of  propagation  of  the  disturbance  and  its 
period  appears  to  be  hyperbolic  in  nature.  Rostoker  (1965) 


showed  that  PI  2  mleropulsations  can  he  Described  by  a  weak 
couf  1 1 np;  of  fluid  and  electromagnetic  disturbances  In  which 
the  electromagnetic  contribution  is  dominant.  Any  variation 
in  the  values  of  the  parameters  which  may  occur  in  the  course 
of  propagation  of  the  waves  over  the  earth's  surface  may 
cause  one  or  other  of  the  coupled  mooes  to  become  dominant. 
Thus  alterations  in  ionospheric  conductivity  caused  by  the 
precipitation  of  auroral  particles  will  cause  the  waves  to 
acquire  irregular  phase  velocities  during  the  course  of  their 
propagation  from  the  source  of  the  pulsations  to  the  field 
points.  In  particular,  in  the  case  of  pulsating  auroras 
different  phase  shifts  between  two  stations  would  be  expected 
for  different  sets  of  peaks  in  a  chain  of  pulsations. 

Rostoker  observed  this  to  be  the  case  for  the  horizontal 
components  of  Pi  2's  observed  at  Vancouver,  B.C.,  and  Ralston, 
Alberta. 
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3.  PI  activity  at  Great  Whale  River  and  Byrd  Station 


When  Pi  activity  can  clearly  be  recognized  on  the 
records  as  by  the  abrupt  commencement  of  polar  substorms  they 
occur  within  a  minute  or  two  at  stations  in  and  between  the 
two  auroral  zones  over  large  areas  of  the  earth.  The  changes 
when  sufficiently  abrupt  have  been  called  psc's  (£olar  sudden 
commencement  or  change)  by  analogy  with  the  much  larger 
commencements  (ssc)  of  a  magnetic  storm.  It  is  this  system 
of  psc  events,  the  following  impulsive  Pi  activity  and  the 
preceeding  quiet  interval  that  governs  the  general  pattern  of 
nighttime  activity  at  the  pair  of  conjugate  auroral  zone 
stations.  Great  Whale  River,  Quebec,  and  Byrd  Station, 
Antarctica.  There  is  a  very  persistent  statistical  occurrence 
of  psc’s  and  psi’s  between  0300  and  0*100  UT  at  these  two 
stations.  A  quiet  interval  which  often  preceeds  the  psc  makes 
possible  greater  accuracy  in  the  timing  of  the  event.  Nega¬ 
tive  psc's  without  exception  move  initially  in  the  direction 
corresponding  to  a  decrease  in  H,  in  contrast  to  the  relatively 
few  (in  1963)  daytime  positive  psc’s  which  all  moved  in  the 
direction  of  increasing  H. 

The  fact  that  high  frequency  components  seem  to 
appear  on  the  micropulsation  records  only  when  the  psc  events 
are  sufficiently  abrupt  is  not  new  and  it  has  long  been  known 
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that  the  sudden  commencements  of  magnetic  storms  are  assoc¬ 
iated  with  frequencies  up  to  at  least  3  cps.  It  is  known 
also  that  abrupt  psc's  in  the  auroral  zones  can  be  recognized 
in  the  midlatitudes  of  Canada  almost  simultaneously.  This 
behaviour  suggests  that  there  might  be  a  layer  in  the  iono¬ 
sphere  that  can  be  shock  excited  at  its  natural  frequency  of 
about  1  cps.  It  is  then  natural  to  consider  the  sudden 
event  as  being  caused  by  the  shock  of  a  sudden  change  in 
velocity  of  the  solar  wind,  the  tffect  of  which  is  propagated 
towards  the  earth  by  hydromagnet ic  waves.  The  earlier,  less 
abrupt  event  might  possibly  be  due  to  a  change  in  direction 
of  the  solar  wind  which  probably  occurred  less  suddenly. 

A  comparison  of  the  times  of  individual  negative 
psc's  at  Byrd  and  Great  Whale  has  been  made  by  Jacobs  and 
Wright  (1965).  Examination  of  a  series  of  1^0  days  fron 
September  28.  1963,  to  February  21,  196*1,  yielded  72  cases 
judged  suitable  for  comparison  of  times.  Measurement  was 
made  to  the  nearest  half  minute  from  paper  records  moving  at 
a  speed  of  6  Inches  per  hour.  This  analysis  showed  17  cases 
of  zero  difference,  29  cases  of  positive  difference  and  26 
cases  of  negative  difference  (Byrd  minus  Great  Whale).  Ex¬ 
cluding  the  17  cases  of  zero  difference  the  average  differ¬ 
ences  were  +1.1  minute  and  -1.1  minute.  The  differences  do 
not  conform  to  a  smooth  error  curve  and  personal  errors  in 
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judging  the  times  are  not  wholly  responsible;  instead  we 
believe  that  the  stations  are  not  truly  conjugate  at  the  times 
of  all  negative  psc's. 

Green  (1967)  has  carried  out  a  detailed  study  of 
Pi  1  events  recorded  during  the  IQSY  at  Great  Whale  River, 

Byrd  Station,  College,  Alaska  and  Kiruna,  Sweden.  In  parti¬ 
cular  he  investigated  the  nature  of  the  conjugacy  and  the 
degree  of  similarity  between  Great  Whale  River  and  Byrd 
Station  in  an  effort  to  better  understand  electron  precipi¬ 
tation  from  the  magnetospheric  boundary.  His  study  included 
comparisons  of  dynamic  micropulsation  spectra,  the  relation¬ 
ship  with  concurrent  cosmic  noise  absorption  and  possible 
conductivity  enhancement  effects. 

At  auroral  zone  stations  the  occurrence  of  irregular 
and  burst-like  micropulsations  is  frequently  observed. 

These  are  very  different  from  the  intrinsic  regular  and  con¬ 
tinuous  micropulsations.  Noise  bursts  include  a  wide  range 
of  frequencies  and  their  wave  forms  are  therefore  irregular. 
They  occur  suddenly  and  die  out  after  a  short  time.  They 
are  simultaneously  accompanied  by  the  sudden  enhancement  of 
DS  fields  such  as  negative  bays.  The  main  part  of  a  noise 
burst  decreases  abruptly  away  from  the  central  area  of  activity. 
Thus  frequently  no  evidence  of  a  noise  burst  is  observed  at 
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Churchill  although  the  Great  Whale  records  clearly  show  a 
noise  burst  and  vice  versa.  (The  distance  between  Churchill 
and  Great  Whale  is  about  1100  km.).  At  Victoria,  Vancouver 
and  Ottawa  corresponding  noise  bursts  are  not  found  at  all 
on  records  of  the  same  sensitivity.  On  the  other  hand  Pi  2 
micropulsations  are  sometimes  found  at  the  middle  latitude 
stations  when  noise  bursts  occur  at  Great  Whale  or  Churchill 
(Yanagihara,  1963).  The  corresponding  Pi  2  part  of  the 

t 

noise  bursts  in  the  central  area  of  activity  in  the  auroral 
zone  is  not  clearly  seen  because  it  is  swamped  by  the  higher 
activity.  Generally  speaking  there  is  a  conjugate  relation¬ 
ship  between  Byrd  and  Great  Whale  or  Churchill  for  noise 
bursts.  Sometimes  a  burst  at  Byrd  corresponds  to  one  at 
Great  Whale  and  sometimes  to  one  at  Churchill.  The  conjugacy 
is  thus  rather  loose.  The  loose  conjugacy  can  be  explained 
by  movement  of  the  active  centre  provided  that  the  dimensions 
of  the  active  area  of  the  noise  burst  are  of  the  order  of 


several  hundred  km. 


.  Other  forms  of  PI  activity 

Tepley  (1966)  has  found  that  other  types  of  signals 
may  be  observed,  particularly  during  stormy  periods,  which 
are  less  regular  in  appearance  than  the  continuous  emissions 
but  are  more  regular  than  noise  bursts  (Pi  1).  Signals 
referred  to  as  SIP  and  IPDP  by  Troitskaya  may  possibly  be 
placed  in  this  transition  category.  The  signal  is  character¬ 
ized  by  a  rising  frequency  fine  structure  similar  to  that  of 
hm  emissions  but  the  structural  elements  do  not  repeat 
periodically.  Tepley  refers  to  these  signals  as  storm  time 
emissions  since  all  events  he  has  observed  so  far  have  occurred 
during  magnetically  disturbed  periods,  often  in  conjunction 
with  magnetic  bays.  Less  formally  he  has  called  such  signals 
gurglers  since  when  monitored  aurally  on  time  compressed 
magnetic  tape  with  a  speed  up  factor  of  1000  to  2000,  the 
signals  are  characterized  by  a  sound  similar  to  bubbles  being 
blown  under  water.  Occasionally  rising  frequency  structural 
elements  are  superimposed  on  Pi  emissions  and  a  composite 
signal  may  be  placed  in  the  transition  category,  gurgler-Pi  1. 
Transition  signals  are  characterized  aurally  by  bubbly  sounds 
similar  to  those  made  by  gurglers.  The  occurrence  times  of 
gurglers  and  gurgler-Pi  1  transition  signals  at  Palo  Alto  show 
a  strong  maximum  around  2000  local  time  (0400  UT) .  Helms 
and  Turtle  (1964)  found  an  almost  identical  occurrence  curve 
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for  impulsive  high  frequency  micropulsations  recorded  at 
Byrd  Station,  Antarctica,  and  reported  a  strong  intensifi¬ 
cation  of  visual  auroral  activity  with  each  micropulsation 
burst. 


IPDP  events  frequently  have  structural  elements  in 
the  nolseand  have  been  detected  simultaneously  at  conjugate 
stations  by  Gendrin  et  alia  (1966).  On  the  assumption  that 
these  structural  elements  have  been  generated  at  the  equator 
in  the  magnetosphere  by  an  impulsive  source,  Knaflich  and 
Kenney  (1967)  have  interpreted  the  later  arrival  of  the  higher 
frequencies  as  the  results  of  propagation  in  a  dispersive 
medium.  A  dispersion  analysis  of  structured  IPDP  elements 
observed  during  1965  and  1966  showed  that  they  were  generated 
in  a  narrow  range  of  solar  magnetospheric  longitude  centred 
around  2000  LMT  at  L  values  6 <**13.  The  broad  band  unstruc¬ 
tured  noise  seems  to  originate  in  a  zone  of  magnetic  turbul¬ 
ence  just  outside  the  normal  trapping  region  so  that  the 
structured  elements  could  be  the  dispersed  signal  resulting 
from  a  broad -band  impulse.  That  part  of  the  magnetosphere 
where  turbulence  exists  is  also  the  same  region  where  the 
field  lines  cease  to  be  closed  and  become  transported  out  to 
the  tail  region.  Hence  it  is  possible  that  the  broad-band 
Impulses  that  cause  the  amplitude  enhancements  are  merely 
manifestations  of  the  instabilities  associated  with  opening 
and  closing  of  field  lines  in  a  critical  region. 
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Heacock  (1967)  has  distinguished  two  subtypes  of  Pi. 
The  first  subtype  are  midnight  Pi  bursts.  He  found  (1966) 
that  most  weak  impulsive  noise  bursts  occurring  on  summer 
nights  showed  an  enhancement  in  the  2  to  6  sec.  period  range. 
There  seems  to  be  an  annual  variation  in  the  mid-frequency 
of  the  enhancement  with  the  lowest  value  occurring  in  summer  - 
the  mid-frequency  varying  from  approximately  0.2  cps  in  day- 
t  ime  in  summer  to  0. 3-1.0  cps  during  dark  nights  in  fall  and 
winter.  Bursts  tend  to  be  seen  clearly  on  the  records  only 
when  Kp  is  less  than  or  equal  to  3.  At  times  of  high  Kp, 
bursts  seem  to  be  superimposed  on  other  Pi  activity  and  cannot 
be  clearly  distinguished.  Sudden  commencements  and  sudden 
impulses  that  occur  during  the  nighttime  at  College  rarely 
coincide  with  the  onset  of  micropulsation  bursts  (although 
ssc’s  that  occur  at  College  during  the  daytime  usually  coin¬ 
cide  with  the  onset  of  micropulsation  events). 

In  addition  to  impulsive  bursts,  the  records  of  Pi 
activity  at  College  contain  events  which  have  a  continuous 
non-impulsive  character,  designated  by  Heacock  as  Pi  (c). 

A  subclass  (elementary  Pi  (c))  are  simple,  isolated  events 
characterized  by  the  auroral  electrojet  being  nearly  station¬ 
ary  just  north  of  College.  These  events  are  characterized 
by  simultaneous  negative  bays  in  the  magnetic  H  and  Z  compon¬ 
ents.  Heacock  found  that  every  elementary  Pi  (c)  event  which 
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began  in  the  interval  2300  -  0200  LMT  was  inaugurated  with  at 
least  one  Pi  burst.  Pi  (c)  can  occur  any  time  of  the  day, 
but  are  most  intense  in  early  morning  near  0400  LMT.  It  was 
found  that  at  high  Kp  values,  Pi  occurrence  is  increasingly 
dominated  by  Pi  (c). 

5.  Theories  of  Pi's 

Pc  theories  are  quite  different  from  Pi  theories. 

Pc  theories  generally  depend  on  some  kind  of  resonance  to 
achieve  the  narrow  band  effect.  Although  a  number  of  sources 
for  broad  band  geomagnetic  noise  have  been  suggested,  in¬ 
cluding  auroral  electrojets  and  the  action  of  the  solar  wind 
at  the  surface  of  the  magnetosphere,  the  only  detailed  Pi 
theory  is  that  of  Nishida  (1964).  Nishida  showed  that  the 
presence  of  a  passing  electron  beam  through  the  magnetospheric 
plasma  modifies  the  character  of  hydromagnetic  waves  in  such 
a  way  that  an  instability  takes  place.  The  situation  is 
reminiscent  of  the  hydromagnetic  instability  due  to  run-away 
electrons  encountered  in  Project  Matterhorn.  With  adequately 
chosen  parameters  for  the  beam  and  the  ambient  plasma,  Nishida 
estimated  the  range  of  the  period  of  growing  waves  generated 
in  the  magnetosphere  and  found  close  agreement  with  that  of  a 
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pulsation  burst.  As  the  bombardment  of  the  Ionosphere  by 
the  electron  beam  continues  the  electrics!  conductivity  of  the 
Ionosphere  would  be  enhanced  and  hydromagnetlc  waves  Incident 
from  above  would  be  shielded  to  an  Increasing  degree.  This 
might  be  the  explanation  of  why  pulsation  bursts  are  observed 
only  In  the  early  stages  of  the  development  of  a  magnetic  bay. 

It  might  also  explain  why  the  longer  period  part  of  the  burst 
lasts  for  a  longer  time.  Nishida  suggests  two  mechanisms  as 
possible  causes  of  PI  1  and  PI  2  in  middle  and  low  latitudes. 

The  first  Is  the  oblique  propagation  of  unstable  hydromagnetlc 
waves  to  lower  latitudes.  The  second  Is  the  leakage  of  the 
screening  current  which  flows  In  the  auroral  zone  when  the 
pulsation  burst  passes  through  the  ionosphere  to  lower  latitudes. 
The  occurrence  of  PI  2  Is  less  frequent  In  solar  active  years 
than  in  quiet  years  although  the  activity  of  PI  1  is  higher 
In  solar  active  years.  Because  the  Ionospheric  conductivity 
would  be  higher  and  the  degree  of  ionospheric  screening  and 
the  Intensity  of  the  screening  current  would  also  be  higher 
In  solar  active  years  It  might  be  possible  to  account  for 
PI  2  by  oblique  propagation  and  PI  1  by  the  leakage  of  the 
screening  current. 

Although  there  are  many  examples  of  good  correlation 
between  micropulsations  at  conjugate  points,  it  Is  not  obvious 
that  there  should  be  perfect  correlation  since  conjugate 
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points  are  defined  by  the  geometry  of  the  static  field. 

There  are  many  Important  questions  requiring  answers,  some  of 
which  are:  How  much  of  the  lack  of  conjugacy  Is  due  to  mag¬ 
netic  activity?  How  much  to  the  magnitude  of  the  disturbance, 
which  we  might  assume  determines  the  area  of  conjugate 
relationship  and  how  much  to  the  positions  of  the  stations 
vls-a-vls  the  angle  of  the  solar  wind?  To  settle  these 
questions  more  statistical  data  are  required,  particularly 
at  different  seasons.  Limited  data  may  give  quite  mislead¬ 
ing  results. 
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Figure  3.  Number  of  frequency  components 
in  Pi  2  events  as  a  function  of 
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Conjugate  Observations  of  the  pi  1  Micropulsation  Event 
of  May  25,  1964 

by 

R.  G.  Green 

Institute  for  Telecommunication  Sciences  and  Aeronomy 
Environmental  Science  Services  Administration 
Boulder,  Colorado 

1.  Introduction 

This  is  a  brief  account  of  an  IQSY  event  that  probably  originated 
with  the  108  MHz  burst  of  May  21^11^46m  U.T.  Conjugate  ULF  pulsations 
of  irregular  amplitude  and  period  (l-40  sec)  were  recorded  as  a  micro¬ 
pulsation  storm  at  Baie  St.  Paul,  Quebec,  and  Eights,  Antarctica  for 
several  hours  on  May  25.  These  stations  comprise  a  pair  very  near  to 
the  exact  computed  conjugacy  (L  =  4.0)  as  well  as  being  very  close  to 
each  other  in  local  time  (U.T.  -  5h) •  It  is  the  purpose  of  this  paper 
to  show  the  nature  of  conjugacy  between  this  pair  of  stations.  Comparisons 
of  dynamic  micropulsation  spectra,  simultaneous  amplitude  spectra,  and 
concurrent  cosmic  noise  absorption  are  given. 

2.  Method  of  Analysis 

For  the  measurements  the  magnetic  north-south  component  of  the  magnetic 
flux  density  vector  was  sampled  using  a  two  meter  diameter  antenna  of 
16,000  turns.  The  aniplified  signal  from  the  antenna  was  recorded  with  a 
frequency  modulated ,  magnetic  tape  recording  system.  Sonagrams,  or  fre¬ 
quency  (and  intensity) -time  display,  are  then  produced  for  study. 


3.  Results 


Figure  1  is  such  a  sonagram  of  the  entire  event.  There  is  the  usual 
gt>od  delineation  of  gross  features.  However  it  should  be  noted  that  the 
BSP  record  is  higher  in  sensitivity  than  the  El.  Note  the. general  impulsive 
behavior  and  broad  spectrum  of  the  frequency  components.  Also  shown  are 
the  concurrent  changes  in  the  30  MHz  riometer  ionospheric  absorption  in 
decibels  above  the  quiet  day  values.  Note  the  degree  of  dependence  between 
frequency  span  and  absorption.  Figure  2  shows  the  event  on  the  same  time 
scale  but  with  an  expanded  sonagraph  frequency  display  of  the  portion  below 

0.12  cps.  In  quite  another  way,  through  ionospheric  conductivity  enhance- 

* 

ment  effects,  certain  features  may  not  be  at  all  comparable.  For  practically 
the  entire  event  BSP  was  in  sunlight,  and  El  was  in  darkness.  There  being 
no  transitional  sunlight  period  for  e: 4  uer  station  during  the  storm  it  is 
not  possible  at  this  time  to  point  out  any  obvious  conductivity  enhancement 
effects. 

< 

One  may,  upon  examination  of  figures  1  and  2,  note  many  groups  of 
spectral  patches  concurrent  at  the  two  sites.  Figure  3  is  provided  to  show 
one  such  group. 

A  display  of  simultaneous  amplitude  spectra  is  useful  to  determine 
relative  frequency  enhancement.  Such  a  display  below  0.5  cps  at  30® 
intervals  during  the  event  is  shown  in  figure  k. 
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Figure  1.  Scuagraph  apectral  dlaplagr  of  conjugate  behavior  of  signals  at  Bat*  t.  Pul  nd  Eight*  ststlon* 
Tram  06l6  to  llM  U.I.  on  Nay  25,  1964.  Al*o  aboun  (above  and  brio*)  ar*  tba  ooocmrrant  change 
la  30  Mb  rlcnatsr  bUnflicc  abort  th*  quiet  day  vain*.  On  the  BSP  aonagrau  note  lha  baavy 
sinusoidal  nolaa  band*  occurring  at  2CR  Intervals  abova  0.5  cpa  at  vail  a*-  tba  artlflcally 
produced  vertical  atriationa  elaaabar*.  Slntlarly  on  the  SI  portion  nolaa  bandi  ar*  aaan  ait 
0.Q3,  0.1,  0,3,  and  0.6  cpa. 
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Hfbm  8.  A  portion  of  tba  pravloua  figure  expanded  to  ahov  detail  lalow  0.12  cpa.  On  tba  BSP  mum 

*"  artlflcally  produced,  (hi  the  El  aonagran  nolaa  band*  ar*  rrldant 
at  0.03  and  0.1  cpa.  , 
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Abstract 

Magnetic  micropulsation  and  energetic  electron  precipitation  data 

^  ’ 

recorded  simultaneously  in  the  auroral  zone  have  been  used  to  charac¬ 
terize  distinct  types  of  pulsation  and  precipitation  phenomenal  'It  is 
found  that  certain  types  of  pulsation  and  precipitation  are  invariably 
associated  with  each  other.  Further,  these  phenomena  occur  only  in  cer¬ 
tain  limited  regions  of  local  time  and  are  always  associated  with  an 
auroral  substorm.  Correlated  phenomena  include: 

Local  Time  2200  -  0200,  Micropulsation  noise  bursts  and  impulsive 
electrons,  associated  with  auroral  breakup  and  sudden  changes 
in  auroral  electro,}  et. 

Local  Time  0200  -  1000,  Band  limited  irregular  pulsations  and 
modulated  electrons,  associated  with  diffuse,  patchy  aurora 
and  negative  bays. 

Local  Time  1000  -  1500,  Quasi-sinusoidal  pulsations  and  modulated 
electrons,  associated  with  large  amplitude  dayside  magnetic 
disturbances  occurring  during  auroral  cubstorms. 
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Local  lime  lUOO  -  11)00,  Magnetic  impulses  and  microburst  elec¬ 
trons,  associated  with  dayside  substorm  disturbances. 

Local  Time  1500  -  2200,  Micropulsation  sweepers  and  pre-midnight 
pulsations,  associated  with  no  type  of  structured  electron 
precipitation,  but  occurring  only  during  auroral  substorms. 

Evidence  is  given  substantiating  the  cc /'relation  between  these  micro¬ 
pulsation  and  energetic  electron  precipitation  phenomena,  and  linking  them 
to  the  auroral  substorm.  This  evidence  suggests  a  causal  relation  between 
modulated  electron  precipitation  are  associated  micropulsation  types  with 
micropulsations  playing  the  more  fundamental  role. 


Introduction 


Several  of  us  from  the  University  of  California  are  studying  the 
characteristics  of  auroral  zone  micropul  cations  and  energetic  electron 
precipitation.  We  have  been  motivated  by  the  desire  to  establish  two 
concepts.  First,  we  believe  a  causal  relationship  exists  between  certain 
types  of  micropulsations  and  certain  forms  of  electron  precipitation. 
Second,  we  feel  that  both  micropulsations  and  electron  precipitation  can 
be  included  in  a  comprehensive  model  of  auroral  zone  disturbances.  To 
demonstrate  the  truth  of  these  concepts  we  have  performed  simultaneous 
measurements  of  micropulsations  and  electron  precipitation  at  the  auroral 
zone  station,  Flin  Flon,  Manitoba.  Data  were  acquired  in  August  and 
September  of  1965  and  1966.  Our  analysis  of  these  data  supports  the 
following  conclusions: 

1.  Auroral  zone  micropulsations  and  energetic  electron  precipitation 
can  be  divided  into  distinct  types  which  usually  occur  only  in 
limited  local  time  intervals. 

2.  Certain  types  of  micropulsations  and  electron  precipitation 
always  occur  together.  Some  of  these  appear  to  be  causally 
related. 

3.  Most  types  of  micropulsations  and  electron  precipitation  appear 
at  their  respective  local  times  only  when  an  auroral  substorm 
occurs  in  the  midnight  sector. 

h.  Dynamical  processes  associated  with  auroral  substorms  are  not 
limited  to  the  auroral  ionosphere  in  the  midnight  sector.  They 
appear  completely  around  the  euroral  zone,  and  throughout  the 
magnetosphere . 
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The  date  and  arguments  which  justify  these  conclusions  in  detail 
will  nppear  in  a  sequence  of  three  papers:  Parks  et  al,  ( 19&7 ) »  McPherron 
et  al.  (1967)  and  Coroniti  et  al.  (1967).  In  this  paper  we  wish  to  illus¬ 
trate  the  type  of  analysis  we  have  carried  out  and  to  summarize  our  con¬ 
clusions. 

Instrumentation  and  Analysis 

X-ray  detectors  carried  by  balloons  to  very  high  altitudes  have  pro¬ 
vided  us  with  precise  data  on  the  temporal  characteristics  of  energetic 
electron  precipitation.  Various  geometrical  arrangements  of  detectors 
and  pulse  height  analyzers  have  given  information  about  the  spatial  di¬ 
mensions  and  energy  spectra  of  tv  precipitating  electrons.  A  pair  of 
induction  coils  and  associated  electronics  at  the  balloon  launch  site 
detected  magnetic  variations.  Both  micropulsation  and  precipitation  data 
were  recorded  continuously  on  an  eight  channel  recorder  and  also  on  FM  tape. 

3y  visual  examination  of  these  records  we  identified  a  number  of  dis¬ 
tinct  types  of  micropulsaticns  and  precipitation.  Interesting  events  were 
replayed  from  the  FM  tapes  ana  digitized.  Fourier  and  dynamic  spectral 
analysis  were  carried  out  on  the  computer  to  study  the  frequency  charac¬ 
teristics  of  pulsating  signals.  Numerical  calculations  were  performed 
tc  determine  source  size  and  e-folding  energy  of  the  incident  electrons. 

Using  riometer  and  earth  current  data  from  permanent  auroral  zone 
stations,  we  extended  our  electron  precipitation  and  micropulsation  data 
to  cover,  simultaneously,  half  of  the  auroral  zone.  Standard  magneto- 
~rcms  from  a  number  of  auroral  zone  observatories  demonstrated  the  occur¬ 
rence  of  auroral  substorms. 
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Experimental  Results 


(a)  Typos  of  Micropulsations  and  Electron  Precipitation 

Visual  examination  of  our  records  shows  both  micropulsations 
and  electron  precipitation  have  a  well  defined  diurnal  occurrence  pattern. 
We  find  also  certain  forms  of  electron  precipitation  always  accompany 
definite  micropulsation  types.  A  sumnary  of  our  results  is  presented  in 
figure  1.  In  a  paper  as  short  as  this  one  we  shall  not  be  able  to  prove 
in  detail  the  truth  of  this  picture.  Detailed  presentations  were  given 
at  the  April  1967  meeting  of  the  American  Geophysical  Union  in  Washington, 
D.  C.  Our  results  are  consistent  with  previous  work  which  has  appeared^, 
in  the  literature. 

The  region  around  local  midnight  is  dominated  by  the  most  violent 
effects  of  the  auroral  substorms.  The  micropulsations  are  predominantly 
in  the  form  of  broad  band  noise  bursts  which  are  accompanied  by  impulsive 
electron  precipitation.  Fifty  millisecond  bursts  of  electrons  are  also 
observed  at  this  time.  Progressing  in  local  time  from  0200  to  1000, 
micropulsaticn  activity  is  in  the  form  of  band  limited  irregular  pulsations 
and  is  associated  with  modulated  electron  precipitation  of  5  to  10  second 
period.  In  the  region  of  local  time  from  1000  to  1500,  the  characteristic 
type  of  micropulsation  is  the  quasi-sinusoidal  oscillation  from  15  to  40 
second  period.  The  electron  precipitation  is  modulated  at  about  the  same 
period.  Also,  around  local  noon  magnetic  impulses  are  associated  with 
microbursts.  From  1500  to  2200  local  time  three  types  of  raicropulsations 
are  seen.  Auroral  zone  pearl  pulsations  peak  around  1500.  Around  the 
dusk  meridian,  the  so-called  sweeper  of  I.P.D.F.  events  are  seen.  Some¬ 
what  later  in  local  time  pre-midnight  pulsations  occur.  None  of  these 
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throe  micropulnation  types  are  associated  with  electron  precipitation. 

(b)  Causal  Relation  between  Micropulsations  and  Electron  Precipitation 
Wc  have  indicated  in  figure  1  that  micropuluations  and  electrons 
are  correlated.  Our  analysis  suggests  that  in  some  cases  a  causal  rela¬ 
tionship  exists  as  well.  To  illustrate  this  point  we  would  like  to 
present  the  most  convincing  example  we  have  studied  so  far.  Figure- 2 
shows  micropulsations  of  20-U0  second  period  beginning  at  1216  (90°  WMT), 
September  8,  1966.  Modulated  electron  precipitation  of  similar  period 
accompanies  the  micropulsations.  This  particular  event  was  initiated 
by  two  5 -minute  sequences  of  nearly  pure  ^0  second  pulsations  at  1216 
and  1220.  The  main  pulsation  ever+  oegan  at  1232,  increasing  in  20 
minutes  to  a  maximum  amplitude  of  about  3  gammas.  Thereafter,  the  signal 
amplitude  decreased  slowly  until  a  new  micropulsation  event  began  at  1630. 
Obvious  modulated  electron  precipitation  of  20-40  second  period  began  at 
the  same  time  as  the  main  micropulsation  event  and  persisted  throughout. 

Figure  3a  and  b  shows  the  results  of  Fourier  analysis  of  both  the 
X-ray  and  micropulsation  signals.  Part  a  contains  the  entire  frequency 
spectra  from  DC  to  10  second  period  plotted  on  a  log  scale.  The  micro- 

0 

> 

pulsation  spectrum  is  made  up  of  a  number  of  well  defined  spectral  peaks 
superimposed  on  a  falling  spectrum.  The  X-ray  spectrum,  however,  falls 
very  rapidly  from  DC  to  about  100  seconds  and  then  appears  to  be  made  up 
of  less  well  defined  peaks  on  a  more  slowly  decreasing  background. 

Spectral  peaks  which  appear  to  correlate  have  been  connected  by  vertical 
dashed  lines.  The  period  of  significant  peaks  in  the  micropulsation 
spectrum  is  indicated  above  each  peak.  Portions  of  the  spectra  including 
the  40  second  peaks  have  been  plotted  on  a  linear  scales  shown  in  part  b 
of  figure  3. 
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Visual  examination  of  the  chart  records  suggests  that  the  micro- 
pulsations  X-rays  are  correlated  at  a  period  of  25  seconds.  The 
spectra  confix*  this  observation  but  also  indicate  correlation  at  20, 

37  and  170  seconds.  It  should  be  noted  that  the  nicropulsation  peak 
near  40  seconds  has  fine  structure.  Furthermore,  the  uain  peak  at  4l 
seconds  does  not  appear  to  be  correlated  with  X-rays,  while  the  subsi¬ 
diary  peak  at  37  seconds  does. 

Figure  4  illustrates  the  results  of  dynamic  spectral  analysis  of 
both  signals.  All  the  features  described  above  are  clearly  evident  in 
the  sonagram. 

(c)  Relation  of  Micropulsations  and  Electron  Precipitation  to  the 
Auroral  Sub  storm 

Data  obtained  at  one  station  are  sufficient  to  determine  the 
diurnal  occurrence  pattern  of  the  auroral  zone  ai  cropulsat ions  and  elec¬ 
tron  precipitation.  However,  these  data  cannot  give  any  indication  of 
their  spatial  extent.  Neither  can  data  at  a  single  station  always  show 
the  relation  of  these  phenomena  to  auroral  substorms.  The  auroral 
breakup  and  magnetic  bay  disturbances  which  characterize  an  auroral 
substorm  do  not  generally  occur  in  the  day  hemisphere.  To  extend  our 
data  coverage  we  have  obtained  riometer,  magnetometer  and  earth  current 
records  from  two  stations  east  and  west  of  Flin  Flon.  Kiruna,  Sweden, 
eight  hours  ahead,  and  College,  Alaska,  three  hours  behind  Flin  Flon, 
include  eleven  hours  in  local  time.  The  relative  positions  of  these 
stations  at  1500  UT  are  shown  in  figure  5.  College  is  at  0500,  Flin 
Flon  at  0800  and  Kiruna  at  1600  local  time. 
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Observations  maiU>  at  these  three  stations  on  the  19th  of  August,  19&5, 
are  shown  in  figure  (i.  The  top  three  curves  indicate  the  occurrence  of 
energetic  electron  precipitation  at  these  stations  as  detected  by  rio- 
meters  or  balloons.  The  fourth  block  contains  a  plot  of  five  minute 
scalings  of  micropulsation  amplitude  and  frequency  at  Flin  Flon.  The 
bottom  set  of  curves  are  the  horizontal  magnetic  field  variations  at  the 
three  stations.  The  occurrence  of  at  least  six  substorms  in  this  univer¬ 
sal  day  is  obvious  from  the  characteristic  bay  disturbances.  Comparison 
of  the  Flin  Flon  micropulsation  and  precipitation  data  with  the  magneto¬ 
meter  traces  leads  to  the  obvious  conclusion:  every  major  micropulsation 
or  precipitation  event  occurs  during  an  auroral  substorm  in  the  midnight 
sector.  Although  it  is  not  evident  from  this  figure,  ve  find  that  there 
is  one  exception  to  this  rule.  Pearl  pulsations  do  not  appear  to  be  re¬ 
lated  to  substorms.  All  other  micropulsation  precipitation  types  mentioned 
above  are  related. 

These  data  support  a  second  conclusion  as  well.  Consider  the  substorm 
which  occurs  over  central  Russia  at  1?00  UT  as  indicated  by  a  negative 
bey  at  College.  Electron  precipitation  begins  first  at  College,  next  at 
Flin  Flon,  and  later  at  Kiruna.  By  1"00  UT,  precipitation  is  occurring 
simultaneously  over  at  least  half  the  aurcral  zone.  By  inference  from 
our  midnight  data  we  suggest  that  it  is  also  occurring  over  Russia. 

Thus,  we  conclude  that  electron  precipitation  appears  nearly  simultan¬ 
eously  throughout  the  entire  auroral  zone  when  an  auroral  substorm  occurs 
in  the  midnight  sector. 

Micrcpulsations  of  the  band  limited  irregular  type  occur  at  Flin 
Flcr.  during  this  subetorm.  College  earth  currents  also  show  the  presence 
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of  these  pul  nations.  Kiruna  is  beyond  the  appropriate  local  time  in¬ 
terval  and  no  significant  activity  is  observed.  Thus,  aicropulsations 
also  occur  simultaneously  over  a  more  limited  region  during  the  course 
of  an  auroral  substorm. 

Conclusions 

Before  concluding,  we  summarize  again  the  results  of  our  work: 

1.  Auroral  zone  micropulsations  and  energetic  electric  precipi¬ 
tation  can  be  divided  into  distinct  types  which  usually  occur 
only  in  limited  local  time  intervals. 

2.  Certain  types  of  micropulsations  and  electron  precipitation 
always  occur  together.  Some  of  these  appear  to  be  causally 
related. 

3.  Most  types  of  micropulsations  and  electron  precipitation 
appear  at  their  respective  local  times  only  when  an  auroral 
substorm  occurs  in  the  midnight  sector. 

U.  Dynamical  processes  associated  with  auroral  substorms  are 

not  limited  to  the  auroral  ionosphere  in  the  midnight  sector. 

j 

They  appear  completely  around  the  auroral  zone  and  through¬ 
out  the  magnetosphere. 

The  fact  that  the  occurrence  of  an  auroral  substorm  is  present  not 
only  in  auroral  disturbances  but  also  in  electron  precipitation  and 
micropulsation  activity  throughout  the  auroral  zone,  indicates  that  it 
is  the  dynamical  processes  occurring  in  the  magnetosphere  which  deter- 
.  mine  the  local  time  characteristics  of  the  substorm.  Recent  satellite 

observations  have  shown  that  many  magnetospheric  phenomena  are  correlated 
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with  substorms.  In  order  to  generalise  the  concept  of  the  auroral  sub- 
(jorm  to  include  the  world-wide  disturbance  characteristics  and  to  em-. 
phasize  the  importance  of  the  magnetosphere  in  auroral  zone  observations, 
we  suggest  a  new  terminology:  magnetospheric  substorm. 
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Figure  Captions 


1.  Diurnal  occurrence  pattern  for  various  types  of  micro¬ 
pulsations  and  electron  precipitations  showing  correlated 
phenomena. 

2.  Quasi-sinusoidal  micropulsations  and  modulated  electron 
precipitation  of  25-40  second  period,  1216  (90°  WMT), 
September  8,  1966. 

3.  Fourier  spectra  for  mi rrv -pulsations  and  modulated  electron 
precipitation,  1240-1440  (90°  WMC),  September  8,  1966. 

4.  Digital  dynamic  spectra,  power  spectra  and  computer  plot 

of  digitized  data  for  micropulsations  and  modulated  electron 
precipitation,  1240-144-0  (90°  WMT),  September  8,  1966. 

5.  Locations  of  various  auroral  ocr.c  stations  in  '’oral  time 
for  substorm  of  1500  UT,  August  19,  1965. 

c.  Plot  showing  world-wide  extent  and  relation  of  electron  pre¬ 
cipitation  and  micropulsations  to  auroral  substorms,  August 

19,  1965. 
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POLARIZATION  ANALYSIS  OF  NATURAL  AND 
ARTIFICIALLY  INDUCED  GEOMAGNETIC  MICROPULSATIONS 

by 

R.  A.  Fowler 
B.  J.  Kotlck 

The  polarization  characteristics  of  geomagnetic  micropulsations  are 
investigated  using  quasi-monochromatic  wave  train  theory  of  physical 
optics.  The  polarization  of  a  wave  train  is  analyzed  by  considering 
the  coherency  matrix  of  the  wave  field.  The  elements  of  the  matric 
are  calculated  over  the  desired  frequency  band  and  the  polarization 
parameters  evaluated.  Since  the  determinant  of  the  coherency  matrix 
is  zero  for  a  totally  polarized  wave  field,  this  property  can  be  used 
to  separate  a  partially  polarized  wave  field  into  its  polarized  and 
unpolarized  parts.  The  ellipticity,  orientation  and  sense  of  polar¬ 
ization  for  the  polarized  portion  are  then  calculated  in  terms  of  the 
martix  elements  using  power  and  cross-power  spectral  analysis  techniques. 
Based  on  this  approach,  the  polarization  of  a  number  of  different  signals 
has  been  digitally  determined.  Test  cases  include  a  random  number  time 
series,  a  constructed  signal  with  a  predetermined  polarization,  and  an 
example  typical  of  geomagnetic  micropulsations  resulting  from  high 
altitude  nuclear  weapons  tests  superimposed  on  high  natural  backgrounds 
The  results  of  these  analyses  are  then  compared  with  those  from  the 
standard  hodogram  technique  and  the  effectiveness  of  both  methods 
discussed  in  terms  of  the  signal  characteristics. 
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Frequency  Analysis  of  the  (.Vomagnetlc  Micropulsations 
Observed  at  the  Magnetic  Equator 

Yoshlo  Kato 

Geophysical  Institute  of  Tohoku  University 
Sendai,  Japan 

The  conjugate  phenomena  concerning  geomagnetic  micropulsations  have 
been  investigated  by  many  workers  and  have  been  reported  by  Campbell  and 
others  in  this  symposium.  However,  at  the  present  time  characteristics 
of  the  geomagnetic  micropulsations  at  the  magnetic  equator  have  not  been 
so  thoroughly  investigated. 

Prince  and  Bostick  (1964)  calculated  theoretically  the  power  spectrum 
of  micropulsations  of  the  extraordinary  mode,  propagating  at  90°  with  respect 
to  the  earth's  magnetic  field  from  the  outer  boundary  of  the  magnetosphere. 
Field  and  Greiflnger  (1965)  calculated  the  transmission  of  micropulsations 
of  ordinary  mode  which  propagate  parallel  to  the  geomagnetic  field. 

Last  year  we  observed  the  magnetic  micropulsations  at  Huancayo,  Peru 
whose  location  is  near  the  magnetic  equator.  This  enabled  us  to  compare 
the  observed  power  spectrum  of  the  micropulsations  at  the  magnetic  equator 
with  that  of  Maui  Island,  Hawaii,  and  Onagawa,  Japan,  in  the  middle  latitude. 
We  also  compared  the  same  magnetic  micropulsations  with  the  theoretical 
power  spectrum  of  the  results  by  the  above-mentioned  workers. 

Before  talking  about  these  results  the  occurrence  mechanism  of  pi  2 
using  the  data  of  its  power  spectrum  is  proposed.  Figs.  1  and  2  show 
examples  of  records  of  the  frequency  analysis  of  pi  2.  These  are  reproduced 
through  several  band  pass  filters. 
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It  is  Interesting  to  note  that  pi  2  has  a  broad  band  of  periods  from 
1  sec  to  210  sec,  and  it  seems  as  if  the  maximum  amplitude  of  the  individual 
bands  nearly  form  a  series  of  harmonic  periods  as  shown  in  the  power  spectrum 
of  pi  2  in  Figs.  3  and  4. 

Concerning  the  mechanisms  of  the  excitation  of  pi  2,  the  following  were 
proposed  by  the  author.  We  can  expect  a  kind  of  instability  (a  neutral  point 
instability)  in  the  transition  region  of  the  neutral  sheet  of  the  magnetic 
tail  (Fig.  5).  The  energy  due  to  the  deformation  of  a  dipole  magnetic  field 
is  stored  within  the  magnetosphere,  and  the  growth  of  the  instability  in  the 
transition  region  is  necessarily  accompanied  by  the  release  of  stored  energy. 
After  releasing  the  energy  for  thei  deformation,  the  magnetic  lines  of 
force  around  the  neutral  point  would  then  be  compressed  toward  the  earth. 

The  energy  released  is  then  used  to  accelerate  plasma  particles  trapped 
there  and  to  generate  the  hydromagnetlc  shock  wave. 

This  hydromagnetlc  shock  wave  will  reach  the  outermost  side  of  the 
magnetic  cavity  where  the  hydromagnetlc  waves  will  be  excited  by  this  stress, 
with  not  only  its  fundamental  mode  but  also  together  with  its  higher  harmonic 
periods. 

This  ordinary  wave  will  propagate  parallel  to  the  magnetic  field  line 
and  we  can  observe  the  harmonic  power  spectrum  on  the  earth's  surface. 

Fig.  6  shows  the  record  of  pc  5,  observed  at  Maui  Island  on  the  dayside 
when  typical  pi  2  was  observed  at  Onagawa  on  the  tailside. 

On  the  other  hand,  the  extraordinary  wave  will  propagate  to  the  equator 
and  the  power  spectrum  of  pi  2  is  quite  different  at  the  magnetic  equator  as 
is  shown  in  Fig.  7  or  8  and  Fig.  9  or  10,  showing  that  the  powers  of  short 
period  rapidly  decay  (Fig.  11). 
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Next  the  results  of  frequency  analysis  of  pc  pulsation  at  Huancayo 
are  compared  with  those  at  Maul.  One  of  the  very  clear  characteristics  Is 
that  the  polarization  of  micropulsations  at  the  magnetic  equator  is  linear, 
while  that  at  middle  latitudes  is  circular. 

The  frequency  analysis  and  the  power  spectrum  of  pc  at  Huancayo  are 
shown  in  Figs.  12  and  13.  As  the  figures  show,  the  power  spectrum  of  pc 
at  Huancayo  is  different  from  that  at  Maui.  The  slope  of  power  spectrum  in 
daytime  at  Huancayo  decays  rapidly  and  the  frequency  of  the  most  resonant 
peak  occurs  in  the  short  period  at  night  but  in  the  long  period  during  the 
day. 

This  characteristic  is  consistent  with  the  theoretical  result  presented 
by  Prince  and  Bostick. 

The  frequency  analysis  and  power  spectrum  at  Maui  is  shown  in  Figs.  14 
and  15.  The  slope  of  the  power  spectrum  is  not  as  rapid  as  that  at  Huancayo. 
Also  the  mode  of  the  power  spectrum  during  the  day  does  not  decay  as  rapidly 
as  that  at  night  (Fig.  16).  This  characteristic  seems  also  consistent  with 
Field  and  Greifinger's  result,  but  the  resonant  frequency  peak  is  not  always 
consistent  with  their  theory. 

More  detailed  investigation  not  only  of  theory  but  also  of  simultaneous 
observations  for  the  northern  and  southern  hemispheres  and  the  equatorial 
region  at  high  and  middle  latitudes  to  include  the  day-  and  nightslde  of  the 
earth  are  highly  desirable. 
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SIMULTANEITY  OF  Pc  3  MICROPULSATIONS  AT  CONJUGATE  POINTS 


A.  C.  Fraser-Smith 

Lockheed  Palo  Alto  Research  Laboratory 
Palo  Alto,  California  94304 

In* roduction  Pc  3  micropulsations  are  small  amplitude  (~  gamma) 
oscillatory  geomagnetic  variations  with  periods  in  the  range  of  10-45  seconds. 
In  auroral  and  middle  latitudes  they  occur  predominantly  during  the  day 
(Campbell,  1963),  and  they  are,  therefore,  sometimes  called  "daytime  activity". 
Although  they  have  been  recognized  as  a  distant  class  of  micropulsations  for 
a  number  of  years,  their  origin  is  still  uncertain  and  there  are  few  quan¬ 
titative  results  concerning  their  properties.  Unlike  Pc  1  micropulsations 
they  have  little  apparent  frequency  structure. 

One  property  of  the  Pc  3  micropulsations  which  is  of  importance  in 
distinguishing  between  various  theories  of  their  origin,  and  in  particular 
between  those  involving  poloidal  and  toroidal  oscillations  (Jacobs  &  Westphal, 
1964;  Carovillano  &  Radoski,  1967)  and  other  theories,  is  the  closeness  of 
their  simultaneity  at  conjugate  points.  This  simultaneity  is  difficult  to 
measure  accurately  from  chart  records  because  of  the  nature  of  the  waveforms 
and  the  common  presence  of  longer  period  waves  (Pc  4,  for  example).  The 
overlay  technique  to  be  described  in  this  paper  is  tedious  but  can  in  principle 
measure  Pc  3  simultaneity  between  two  stations  to  an  accuracy  of  a  few  seconds. 

The  micropulsation  data  analyzed  for  this  study  were  recorded  during 
1963  on  magnetic  tape  at  the  four  Lockheed  Pacific  geomagnetic  observatories 
at  Palo  Alto,  Kauai,  Canton  Island  and  Tongatapu  (figure  l).  The  three 
island  stations  are  located  at  approximately  the  same  geomagnetic  longitude, 
and  Kauai  and  Tongatapu  are  also  located  at  approximately  conjugate  geomagnetic 
latitudes.  All  observatories  had  identical  recording  equipment.  Pc  3  signals, 


which  were  recognizably  the  cane,  were  found  to  occur  at  all  four  stations. 
However,  those  at  Canton  Island,  which  lies  close  to  the  geomagnetic  equator, 
were  comparatively  very  weak  and  difficult  to  compare  in  detail  with  the 
corresponding  signals  at  the  other  three  stations.  The  Kauai  and  Tonga 
signals  were  always  closely  similar  when  compared  with  those  at  Palo  Alto, 
an  indication  of  the  significance  of  conjugacy  in  P>;  3  propagation  (figure 
2).  Two  periods  of  Pc  3  activity  were  analyzed  in  detail.  These  occurred 
during  August  1-2,  1963,  and  August  27-23,  19-0,  when  there  were  approximately 
seven  and  nine  hours,  respectively,  of  comparable  activity.  Altogether,  a 
total  of  nearly  100  hours  of  recorded  acta  from  the  four  stations  were 
analysed.  The  results  to  be  presen  *  here  mostly  concern  the  two  conjugate 
stations .  Some  comparison  will  lq  made,  however,  with  the  non-conjugate 
i.-'a.o  Alto  station. 

Method  of  Analysis  Following  a  speed-up  of  about  2000,  the  micropulsation 
data  were  frequency  analyzed  using  a  Kay  Missilyzer  and  Scale  Expander.  This 
enabled  the  Pc  3  signals  to  be  separated  from  other  micropulsation  activity. 
Time  marks  were  recorded  on  one  channel  of  the  Missilyzer  and  the  data  on 
the  other,  a  procedure  enabling  a  close  connection  between  the  Pc  3  signals 
and  their  times  of  occurrence  to  be  maintained.  Only  North-South  component 
data  were  used,  and  all  these  data  were  analyzed  in  exactly  the  same  way. 

Frequency-time  plots  were  made  for  each  hour  of  the  periods  of  Pc  3 
activity  studied,  and  the  Missilyzer  speed  was  adjusted  using  a  frequency 
control  system  to  give  the  same  distance  between  hour  marks  for  each  station 
(figure  3).  The  hour  was  then  divided  into  five  sections  of  approximately 
twelve  minutes  each  which  were  then  photographed  with  polaroid  transparency 
film.  Time  differences  between  corresponding  signals  at  the  different 
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stations  could  then  be  measured  by  overlapping  and  adjusting  the  transparencies 
for  best  fit  (figure  4).  In  general,  the  two  matching  Tongatapu  and  Kauai 
transparencies  were  compared  first,  after  which  the  Palo  Alto  transparency 
(which  was  never  as  good  a  match)  was  compared  separately  with  both  Tongatapu 
and  Kauai.  In  this  way,  a  double  check  of  the  time  differences  was  possible. 
Once  the  transparencies  were  positioned  the  time  differences  could  be 
measured  to  one  second  accuracy.  However,  since  positioning  of  the  trans¬ 
parencies  was  uncertain  within  small  limits,  measurement  accuracy  was 
effectively  2-3  seconds. 

Results  The  results  of  the  measurements  for  the  two  periods  of  Pc  3 
activity  studied  are  shown  in  figures  5  and  6.  All  time  differences  are 
measured  with  respect  to  Kauai.  As  was  to  be  expected,  no  time  differences 
greater  than  a  minute  were  found,  and  on  this  scale  of  time  measurement  the 
Pc  3  micropulsations  are  simultaneous  at  Kauai,  Tongatapu  and  Palo  Alto. 

For  time  measurements  on  the  order  of  seconds,  however,  the  micropulsations 
may  occur  at  different  times  at  the  three  stations.  On  comparing  the 
Palo  Alto-Kauai  (non-conjugate)  results  with  those  for  Tongatapu-Kauai 
(conjugate),  it  appears  that  conjugacy  reduces  but  does  not  eliminate  the 
time  differences  between  stations.  The  following  three  points  are  also 
noted: 

1)  There  is  little  if  any  time  difference  between  Kauai  and  Tongatapu  at 
either  the  beginning  or  end  of  a  period  of  activity. 

2)  Kauai  predominantly  leads  Tonga  in  both  figures  5  and  6.  This  is 
considered  to  be  a  seasonal  effect  and  can  be  expected  to  change  as  the 


year  progresses. 


3)  The  approximate  times  of  local  noon  at  Palo  Alto,  Kauai  and  Tongatapu 
are  as  follows: 

Palo  Alto  noon  ~  2015  UT 

Kauai  "  ~  2245  UT 

Tonga  "  ~  2345  UT 

While  there  are  insufficient  data  to  draw  any  general  conclusions,  there 
appears  to  be  a  connection  between  these  times  and  the  general  trend  of 
time  differences.  For  example,  in  figure  6.  '  e  Palo  Alto  and  Tongatapu 
local  noontimes  just  precede  the  times  of  maximum  time  difference.  Further¬ 
more,  in  both  figures,  there  is  a  small  transition  period  in  the  Tongatapu- 
Xauai  results  preceding  the  Kauai  r.  ..  Some  local  time  control  is 
consistent  with  earlier  observ.v-j.ons  (Troitskaya,  1964). 
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(luring  August  27-23,  1963. 


Figure  3. 


Figure  It. 


Detail  of  coe  Hour  (2200-2300  OT,  August  27(  1963)  of  Ft  3 
sdcrepulsatlcu  activity  at  thr  two  conjugate  etatlcna  Kauai  and 
Toogstapu.  Note  the  close  slnilaiity  between  the  Individual 
signal  bursts. 


Two  over  lapped  polarold  transpareneles  of  the  Pc  3  alero- 
pulsatlans  seen  in  figure  3  araal  2230  UT.  It  is  difficult  to 
reproduce  the  transparency  patterns  and  a  negative  is  shown  for 
greater  clarity. 
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TIME  DIFFERENCE  (SECONDS)  TIME  DIFFERENCE  (SECONDS) 


August  1  - 2, 1963.  Comparison  of  the  times 
of  occurrence  of  Pc  3  activity  at  Tonga 
and  Kauai.  The  time  differences  shown 
are  relative  to  Kauai. 


THE  BEHAVIOR  OF  MAGNETOSPHERIC  PROTON  DENSITIES  DERIVED  FROM 


PEARL  DISPERSION  MEASUREMENTS* 


H.  B.  Liemohn,  J.  F.  Kenney  and  H.  B.  Knaflich 
Boeing  Scientific  Research  Laboratories 
Seattle,  Washington  98124 


The  dispersion  characteristics  of  116  pc  1  micropulsation  pearls 
are  described  and  analyzed  using  various  models  of  the  magnetosphere. 
Theoretical  analyses  have  shown  that  the  dispersion  of  the  signals  de¬ 
pends  primarily  on  the  plasma  density  distribution  along  the  field 
line,  and  the  geocentric  distance  measured  in  earth  radii  to  the  place 
where  the  field  line  crosses  the  equatorial  plane.  Thus  measurements 
of  experimental  dispersion  can  be  used  to  estimate  the  plasma  density 
at  various  distances  out  in  the  magnetosphere.  The  geomagnetic  field 
B  is  assumed  to  be  a  dipole  so  that  the  field-line  paths  are  given  by 
R  *  L  cos2X  where  X  is  the  geomagnetic  latitude  and  L  is  the  geocentric 
distance  to  the  equatorial  intercept  of  the  field.  The  plasma  distri¬ 
bution  is  assumed  to  have  the  general  form  N  *  NQ(L)P(L,X)  where  P(L,0) 
is  normalized  to  unity  so  that  NQ  can  be  interpreted  as  the  equatorial 
density  distribution.  The  different  assumptions  concerning  the  distri¬ 
bution  of  plasma  along  the  field  line  will  change  the  determination  of 
L  only  a  small  amount  while  making  a  larger  change  in  N-, 

* 

To  be  published  in  Earth  and  Planetary  Science  Letters,  July  1967. 
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Figure  1  shows  that  pearl  emissions  are  generated  on  field  lines 
with  equatorial  intercepts  between  four  and  nine  earth  radii.  The 
upper-frequency  limit  rarely  exceeds  0.6  of  the  minimum  ion  cyclotron 
frequency,  f ^  ,  and  the  detectable  bandwidth  is  about  0.1  . 

Figures  2  and  3  show  that  there  is  a  diurnal  variation  in  the  lo¬ 
cation  of  the  propagation  paths  which  constitute  a  pearl  oval  with  high 
latitude  or  large  L-values  on  the  dayside  of  th  ■  earth  and  small  L-val- 
ues  or  low  latitude  on  the  nightside.  Figure  2  shows  that  in  general 
data  taken  during  magnetically  quiet  periods  occur  farther  out  in  the 
magnetosphere  than  those  taken  at  the  ame  local  time  during  the  most 
disturbed  periods.  Figure  3  inch.  „<_es  no  obvious  seasonal  shift  in  the 
location  of  the  pearl  zone. 

Figure  4  shows  the  equatorial  plasma  density,  NQ(L),  for  the 
field-line  distribution  P~R”5.  Determinations  of  N0  from  pearls  which 
are  detected  either  simultaneously  or  within  a  short  time  interval  are 
connected  with  straight-line  segments.  These  straight  line  segments 
indicate  that  pearls  which  are  detected  essentially  at  the  same  time 
may  have  a  spread  of  L-value  in  excess  of  1  earth  radius,  but  usually 
occur  on  a  very  steep  density  gradient  that  may  coincide  with  the  plas- 
mapause  already  detected  by  VLF  techniques. 
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Pearl  frequency  bands .  Fig.  2  —  Diurnal  variation  of  propa1 

gation  paths  Ky  is  local 
magnetic  index. 


Diurnal  variation  of 
propagation  paths. 


Fig.  4  -  Typical  plasma  density 
distribution. 
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Recent  Investigation  on  liydromagnetic  Emissions 
by  Tomiya  Watanabe 


Institute  of  Earth  and  Planetary  Sciences 
University  of  British  Columbia 
Vancouver,  Canada 


1.  INTRODUCTION 

The  words  "hm  emissions"  in  the  title  are  used  in 
a  very  broad  sense.  It  is  implied  that  hm  emissions  are 
synonymous  with  geomagnetic  micropulsations  with  short 
periods,  say,  from  about  0.2  sec  to  5  sec,  that  is,  the 
range  of  the  shortest  periods  ever  observed  on  the  earth's 
surface  for  geomagnetic  micropulsations.  If  hm  emissions 
are  defined  like  this,  all  the  pc  l1  type  of  micropulsations 
are  contained  in  this  category.  Besides,  some  pi  l1  type 
of  micropulsations  are  also  included. 

There  are  various  types  of  signals  in  hm  emissions. 
Signals  of  the  simplest  type  are  hm  whistlers.  Figure  2  of 
Ref.  2  shows  a  typical  example  of  hm  whistlers  displayed  on 
a  spectrogram.  They  appear  as  a  series  of  rising  tones  which 
form  a  fan-shaped  structure.  An  hm  whistler  is  due  to  dis¬ 
persion  of  an  hn  wave  packet  which  propagates  in  the  magneto¬ 
sphere  and  bounces  between  the  twc  hemispheres  along  a  magnetic 
line  of  force  of  the  earth3  4  .  In  general,  there  are  two 
modes  of  waves  in  a  magnetoactive  plasma.  The  mode  relevant 
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to  hm  whistlers  seems  to  be  the  L  mode  (or  the  slow  mode, 
more  exactly).  In  the  hm  frequency  range,  only  the  L  mode 
of  waves  is  likely  to  be  guided  by  the  ambient  magnetic 
field4.  The  R  mode  (the  fast  mode)  of  hm  waves  is  not 
field-guided  in  the  hm  regime.  The  group  delay  of  the  L 
node  of  an  hm  wave  is  larger  for  a  higher  frequency.  This 
explains  the  fan-shaped  structure  of  hm  whistlers. 

Tepley  ani  Wentworth5  mentioned  that  an  hm  emission 
often  starts  as  a  weak  signal  and  as  time  passes  the  signal 
intensity  becomes  continually  stronger  and  in  the  final  stage 
it  dies  out  gradually.  McNicol  and  Johnson6  obtained  an 
example  v.’hich  clearly  shows  such  a  trend  (their  spectrogram 
is  cited  in  Ref.  (15)  ).  Tepley  and  Wentworth5  also  mentioned 
that  in  such  an  hm  emission  event  the  hm  wave  does  not  follow 
the  dispersion  process  simply  out-  that  it  should  be  amplified 
while  propagating  in  the  magnetosphere.  Furthermore,  it 
should  be  noticed  that  there  is  no  fan-shaped  structure  in 
the  event  obtained  by  Johnson  and  McNicol,  unlike  hm  whistlers. 
The  emission  signals  are  more  or  less  parallel  to  each  other. 
This  also  indicates  that  in  this  event  there  is  something  else 
to  be  considered  besides  the  dispersion  process.  Emission 
signals  such  as  this  example  are  quite  common,  and  it  shows 
that  there  is  a  different  group  of  signals  from  hm  whistlers. 
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Those  signals  might  be  called  hm  emissions  in  the  narrow 
sense.  By  adopting  the  nomenclature  this  way,  it  is  implied 
that  hm  emissions  should  be  discriminated  from  hm  whistlers 
as  radio  emissions  in  the  VLF  range  are  differentiated 
from  whistling  atmospherics. 

2 .  Some  Problems  about  the  Cyclotron  Instability  Process 
as  the  Triggering  Mechanism  of  hm  Emissions 

2.1  Model  of  Periodic  Emissions  and  Principal  Results  of 
the  Linear  Theory  of  the  Cyclotron  Instability  Process 

In  1965,  Brice7  mentioned  that  in  hm  emissions  In 

the  narrow  sense  mentioned  above,  there  is  a  type  which  is 

similar  to  periodic  emissions  in  the  VLF  range.  One  might 

km. 

call  it^periodic  emissions.  According  to  Brice8,  the 
mechanism  of  periodic  emissions  is  as  follows.  (See  Fig.  1.) 
Let  us  suppose  that  a  wave  packet  is  propagating  along  a 
magnetic  line  of  force.  In  the  case  of  a  whistler,  this 
wave  packet  is  simply  dispersed  while  propagating  along  the 
magnetic  line  of  force.  What  is  different  in  periodic 
emissions  is  that  somewhere  on  the  propagation  path  there  is 
a  black  box,  a  region  of  triggering  of  hm  emissions.  V.Tien 
the  initial  wave  packet  comes  back  as  an  echo  and  enters  the 
black  box,  it  triggers  a  new  emission.  This  new  emission 
subsequently  triggers  the  next  one  when  it  comes  back  to  the 
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triggering  region.  In  the  cnse  of  VLF  emissions,  the 
triggering  mechanism  is  believed  to  be  the  cyclotron 
instability  process  caused  by  nonthermal  electrons9.  In 
the  case  of  hm  periodic  emissions,  the  triggering  mechanism 
is  also  believed  to  be  the  cyclotron  instability  process. 
However,  nonthermal  protons  rather  than  electrons  are 
believed  to  be  responsible  for  the  instabilities  in  the 
case  of  hm  emissions.  The  same  conclusions  have  been 
obtained  by  some  other  authors3  5  10-1 2  who  Investigated 
the  triggering  mechanism  independently.  According  to  an 
estimate  by  Jacobs  and  the  author  11  12 ,  the  energy  of 
each  nonthermal  proton  is  required  to  be  from  about  a 
few  kev  to  a  few  tens  of  kev,  if  the  triggering  region 
exists  at  a  typical  place  in  the  outer  magnetosphere, 
for  example,  at  six  earth-radii  distance  from  the  centre 
of  the  earth.  On  the  other  hand,  the  required  energy  for 
an  electron  must  be  relativistic  to  cause  an  instability. 

Mathematical  investigation  of  the  cyclotron 
instability  process  is  done  by  looking  for  solutions  of  growing 
waves  in  an  infinite  and  uniform  magneto-active  plasma.  In 
the  case  of  a  plane  wave  where  the  electromagnetic  field 
varies  in  space  and  time  as  ex  jo  i  , 

the  wave  is  said  to  be  a  "growing  wave"  if  the  wave  frequency 
has  a  positive  imaginary  part  while  the  wave  number  % 

9* 

is  taken  real: 


S?n  C  ^  > 
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>  0. 


(l) 


This  type  of  mathematical  analysis  concerns  monochromatic 
waves  only.  In  actuality,  however,  one  has  to  deal  with  a 
wave  packet  that  is  a  superposition  of  many  monochromatic 
waves.  According  to  Sturrock13,  a  growing  wave  packet 
takes  either  of  the  following  two  behaviours.  In  one  case, 
it  is  like  a  travelling  amplifying  wave.  In  the  other  case, 
the  wave  packet  is  like  a  self-exciting  oscillator.  When 
considering  the  model  of  periodic  emissions,  it  is 
important  to  find  out  which  of  the  two  possible  cases  is 
relevant.  If  a  wave  packet  triggered  in  the  black  box  behaves 
as  in  the  first  case,  it  keeps  its  identity  as  a  wave  packet 
while  propagating  through  the  black  box.  However,  if  the 
triggered  wave  packet  behaves  as  in  the  second  case,  it 
fills  the  whole  space  of  the  black  box,  stays  there  and  grows 
as  long  as  nonthermal  particles  supply  their  energy  to  the 
wave  packet .  The  appearance  of  an  hm  emission  signal  on  a 
spectrogram  could  be  quite  different  for  the  two  cases. 

In  the  case  of  VLF  emissions,  it  is  suggested  by  Brice14 
that  the  instability  is  of  the  second  type,  namely,  of  the 
self-exciting  oscillator  type.  For  hm  emissions,  Jacobs 
and  the  author11  12  found  that  the  instability  process  is 
again  of  the  self-exciting  oscillator  type. 
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2.2  Non-linear  Evolution  of  the  Signal  Strength 


This  result  brings  a  difficulty  to  the  proposed 
emission  mechanism,  because  if  this  result  is  correct  the 
appearance  of  an  hm  emission  signal  on  a  spectrogram 
rhould  be  a  continuous  line  as  long  as  the  instability 
process  goes  on.  This  is  contrary  to  what  we  have  observed. 
Jacobs  and  the  writer15  are  of  the  opinion  that  this  difficulty 
can  be  avoided  by  considering  non-linearities  in  the 
cyclotron  instability  process.  Mathematically,  the  non- 
linearities  arise  from  the  Boltzmann  equation  which  is 
needed  to  describe  the  particle  distribution  in  a  plasma 
gas.  According  to  the  linear  theory  of  the  cyclotron 
instability  process,  the  wave  amplitude,  for  instance,  the 
strength  of  the  wave  magnetic  field  b  changes  with  time  as 

7t-rk  (2) 

where  r  is  the  growth  rate  and  becomes  equal  to  the 
Imaginary  part  of  the  wave  frequency  according  to  analysis  in 
the  complex  form: 


r  =  (3> 
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In  the  linear  approximation,  the  growth  rate  is  constant 
with  time,  and  if  it  is  positive,  the  amplitude  b 
grows  indefinitely.  However,  the  growth  rate  r  is  no 
longer  constant  with  time  in  the  non-linear  theory.  Since  the 
wave  is  able  to  grow  only  by  taking  energy  out  of  the 
nonthermal  particles,  th*.'  growth  rate  r  should  decrease 
with  time  as  the  nonthermal  particles  give  up  their  energy. 

As  the  growth  rate  r  tends  to  zero,  the  wave  amplitude 
b  must  tend  to  a  certain  constant  value.  This  is  true, 
however,  only  if  the  temperature  of  the  ambient  plasma 
(excluding  the  nonthermal  particles)  is  zero.  If  the 
temperature  is  not  zero,  the  quasi-linear  theory15  shows 
that  the  above  equation  should  be  rewritten  as  follows. 

^  =  rt>  -  Ab*  ('0 

where  tne  coefficient  A  is  a  positive  constant.  This 
equation  shows  that  the  wave  grows  in  the  initial  stage 
but  as  time  passes  it  decreases  to  zero. 

2.3  Non-linear  Drift  of  the  Emission  Frequency 

Concerning  the  model  of  periodic  emissions,  Brice7 
presented  another  problem  of  why  each  emission  signal  becomes 
a  rising  tone  on  a  spectrogram.  Brice  surmised  that  the 
frequency  rise  might  occur  due,  at  least  partly,  to  the 
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emission  mechanism  itself.  This  conjecture  seems  credible, 
particularly  in  the  case  of  VLP  emissions.  Because  the 
relevant  mode  for  VLF  emissions  is  the  R  mode,  the 
emission  frequency  would  decrease  with  time  if  the  frequency 
change  with  time  were  due  to  dispersion  only.  In  actuality, 
an  individual  signal  in  a  series  of  VLP  periodic  emissions 
is  a  rising  tone  more  often  than  a  falling  tone.  Also  in 
the  case  of  hm  periodic  emissions,  there  is  evidence  to 
support  Brice's  idea.  A  series  of  hm  periodic  emissions 
often  begins  with  a  dim  dot  on  a  spectrogram.  More  often, 
the  starting  point  of  a  series  is  barely  defined.  In  such 
a  case,  there  are  several  or  more  dim  dots  in  the  early 
stage  of  a  periodic  event.  As  t? *  asses,  a  dot  tends  to 
extend  horizontally  mere  and  more,  becoming  more  like  a  line 
segment  rather  than  a  dot.  As  time  goes  on  further,  each 
signal  becomes  a  rising  tone,  and  its  gradient  is  steeper 
as  the  signal  strength  becomes  stronger.  When  a  series 
approaches  the  end,  the  signal  strength  becomes  weaker  and 
weaker  and  subsequent  signals  tend  to  incline  towards  the 
time  axis  more  and  more,  forming  a  fan-shaped  structure. 

An  example  of  an  hm  emission  obtained  by  Gendrin  and  Stefant16 
clearly  shows  such  a  trend. 

Now  the  idea  by  Brice  gives  rise  to  the  question 
as  to  how  the  cyclotron  instability  process  makes,  an  emission 
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frequency  increase  with  time.  Regarding  this  problem, 
Jacobs  and  the  present  author15  are  again  of  the  opinion 
that  non-linearities  in  the  instability  process  might 
give  the  answer.  Suppose  we  are  considering  a  wave 
propagating  along  the  background  magnetic  field.  If  the 
medium  is  a  cold  plasma,  the  wave  is  non-attenuating 
and  circularly  polarized  expressed  as  follows. 


and 


b*  -  b  CM  c  KZ—  OJt  -  6)  f 


b^  —  b  am  txz  —  wt  -  Q ) 


b2  =.  O 


(5) 


The  linear  theory  gives  the  dispersion  equation  by  which 
the  wave  frequency  u  and  the  wave  number  K  are  con¬ 
nected  to  each  other.  For  the  L  (R)  mode,  the  frequency 
u  is  positive  (negative).  The  amplitude  b  as  well  as  the 
phase  angle  e  are  constant  in  time.  If  the  medium  contains 
nontherraal  particles,  however,  the  amplitude  b  is  no 
longer  constant  with  time  but  gradually  changes  with  time 
as  shown  in  Eq.  (2).  The  coefficient  r  is  positive  for 
a  growing  wave.  The  above  equation  can  be  derived  in  the 
linear  approximation.  The  quasi-linear  theory  tells  us 
how  this  equation  must  be  modified,  as  already  mentioned 
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(Eq.  (4)  ).  This  is  not  the  only  result  of  the  quasi- 
linear  theory,  however.  It  also  shows  that  the  phase 
angle  6  is  not  constant  with  time.  This  means  that 
the  effective  frequency  u>'  is  different  from  the  wave 
frequency  w  in  the  linear  approximation,  and  not  only 
that,  the  effective  frequency  might  ’drift'  with  time.  In 
the  case  where  the  ambient  plasma  (excluding  the  non- 
thermal  particles)  is  at  zero  temperature  so  that  the 
non-linear  term  on  the  right  hand  side  of  Eq.  (4)  vanishes, 
according  to  calculation  in  Ref.  (17),  the  rate  of  change 
of  the  effective  frequency  is  given  as  follows 


In  this  expression,  3o  is  the  strength  of  the  ambient 
magnetic  field,  and  C  is  a  constant  which  depends  upon 
the  linear  wave  frequency  u>  and  is  positive  (negative) 
for  the  L  (R)  mode.  This  expression  shows  that  the 
effective  frequency  w'  increases  with  time  as  the  wave 
amplitude  grows  with  time.  The  rate  of  the  frequency 
increase  is  greater  as  the  wave  amplitude  b  and  the  growth 
rate  r  are  greater. 
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3.  Comments  on  an  Unsolved  Problem  in  hm  Emissions 


In  hm  emissions,  there  are  many  other  unsolved 
problems.  The  author  would  like  to  mention  one  of  those 
problems.  (Please  refer  to  Pig.  7  of  Ref.  (2).)  Looking 
at  the  figure,  one  might  wonder  if  this  event  is  a  periodic 
emission,  an  hm  whistler,  or  something  else.  One  way 
to  look  at  this  event  is  to  group  those  dots  lined  up 
horizontally  into  one  series  and  to  regard  it  as  a  series 
of  periodic  emissions.  According  to  this  view,  this 
event  is  an  ensemble  of  several  series  of  periodic  emissions, 
although  these  series  are  not  independent  of  each  other 
but  are  somehow  related.  On  the  other  hand,  one  could 
regard  this  event  as  a  kind  of  whistler,  say,  a  "dotted” 
whistler.  By  connecting  four  or  five  dots,  one  can  trace 
out  a  curve  along  which  the  frequency  increases  with  time. 

A  series  of  curves  traced  out  like  this  look  Just  like 
signals  of  an  hm  whistler. 

If  one  takes  the  second  point  of  view,  an  immediate 
question  is  why  a  whistler  signal  is  not  smooth  but  dotted. 

An  answer  suggested  for  this  question  is  that  the  "dotting" 
occurs  due  to  the  propagation  mechanism  of  hm  waves^. 
Considering  a  field-guided  hm  wave  which  propagates  from  a 
point  deep  in  the  outer  magnetosphere  to  the  earth's  surface. 
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the  transmission  coefficient  of  the  wave  as  a  function 
of  frequency  is  said  to  have  several  maxima  and  minima 
in  a  narrow  frequency  hand.  Whether  this  suggestion  is 
correct  or  not  might  be  checked  by  analysing  the  signals 
of  this  type  with  a  high  resolution  sonograph  or, 
especially,  by  taking  the  signals  on  a  three  dimensional 
spectrogram  where  the  signal  strength  is  displayed  on  the 
third  dimension^  0n  a  three  dimensional  spectrogram  like 
that,  a  dot  must  look  like  a  mountain.  Now,  if  the 
suggestion  given  above  is  correct,  the  edge -line  of  the 
mountain  should  ,run  along  the  whistler  trace.  To  obtain 
a  three  dimensional  spectrogram  of  an  hm  emission  event, 
one  may  employ  a  digital  spectral  analyser  or  a  new 
experimental  technique  which  has  been  developed  by  Kenney19. 
Kitamura  has  been  analysing  a  number  of  hm  emission  events 
using  a  three  dimensional  spectral  analyser.  His  investi¬ 
gation  is  still  under  way.  However,  some  of  the  preliminary 
results  show  that  the  edge-line  of  a  "dot"  mountain  does 
not  necessarily  run  along  the  whistler  trace:  More  often, 
it  runs  oblique  to  the  whistler  trace2?  Prom  this  obser¬ 
vation  Kitamura  and  the  author  agree  that  due  to  this 
a  dot  seems  more  likely  to  be  an  emission  rather  than  a 
segment  of  a  whistler  signal.  However,  why  several  series 
of  periodic  emissions  take  place  simultaneously  is  an 
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unsolved  problem.  It  may  be  that  hm  emissions  of  the 
type  under  consideration  are  something  analogous  to 
multiple  emissions  or  a  chorus  in  the  VLF  range.  Since 
dotted  structures  often  appear  in  hm  emissions,  they 
present  an  important  problem  in  theoretical  investigation 
of  the  emission  mechanism. 
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NOTES 

Concerning  the  problem  on  the  non-linear 
evolution  of  the  signal  strength  (raised  in  the  section  2.2), 
It  is  worthwhile  mentioning  the  recent  result  obtained 
by  I.B.  Bernstein  and  F.  F’  ^elmann  (Phys.  of  Fluids  9, 

937,  1966)  that  in  the  case  of  three  dimensional  electro¬ 
static  plasma  instabilities  in  a  non-magnetic  electron 
plasma  the  fluctuating  electric  fields  decay  effectively 
to  zero  in  a  time  of  the  order  of  the  reciprocal  of  the 
maximum  initial  growth  rate.  It  seems  interesting  to 
extend  the  analysis  by  Bernstein  and  Engelmann  to  the  case 
of  cyclotron  instabilities. 
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Fig.  1  Model  of  HM  Periodic  Emissions 
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ABSTRACT 

Radial  diffusion  coefficients  for  low-energy  protons  in 
the  magnetosphere  are  estimated  for  processes  involving  ?c  1,  2, 
and  3  hydromagnetic  waves •  These  diffusion  coefficients  are  com¬ 
pared  to  rough  estimates  of  diffusion  due  to  magnetic  impulses  und 
hays,  to  arrive  at  an  overall  diffusion  rate*  An  effective  loss 
rate  for  pitch-ungle  diffusion  (caused  by  Pc  1)  i3  developed,  and 
combined  with  radial  diffusion  to  yield  a  nonlinear  differential 
equation,  done  estimates  of  expected  proton  flux  levels  are  aude, 
on  the  basis  of  this  equation. 


IH-10-1 


1.  XSTRGPUCTIOiJ 


Xn  this  paper,  we  discuss  seml-quantitatively  the  ways  In 
which  pitch-angle  scattering  caused  by  hydromagnetic  waves  can  Influence 
the  l.TMr  i  rr.^iul  diffusivi  v  »  j. '-energy  protons.  The  problem  of  radial 
diffusion  is  an  extremely  complicated  one,  with  pitch-angle  diffusion 
taken  into  account,  and  no  definitive  conclusions  can  be  drawn  without 
extensive  computer  studies.  The  results  of  the  present  work  are  neces¬ 
sarily  tentative,  but  we  hope  they  are  sufficiently  reliable  to  use  as 
guidepoits. 

A  number  of  authors  have  discussed  the  gsneration  of  Pc  1 
alcropulsations  (frequencies  0.2  -  5  da)  [Cornwall,  1965*  Obayashi,  196$; 
Brice,  196U}  Jacobs  and  iatanabe ,  196** ].  When  particle  of  the  appro¬ 
priate  energy  have  an  anisotropic  pitch-angle  diutri nation,  an  instability 
results,  with  frequencies  in  the  Pc  1  band  for  Davis-Williamson  [1963, 
1«6U]  protons.  (Electrons  of  £l00  kev  produce  a  similar  instability 
with  whistlers.)  Because  hydromagnet ic  waves  can  bounce  between  conjugate 
points  of  a  magnetic  field  line,  a  single  wave  packet  can  be  aisplified 
many  times.  The  effect  is  reminiscent  of  a  ^-spoiled  laser.  The  micro- 
pulsations  scatter  the  protons  in  pitch  angle,  which  is  best  described 
oe  a  diffusion  process  [Cornwall,  1966;  Kennel  and  Petschek,  1966].  This 
pitch-angle  diffusion  is  an  effective  I033  mechanism  for  Pavi3-Williaoson 
protons . 

We  are  interested  in  the  effect  of  pitch-angle  diffusion  on 
radial  diffusion  [see,  e.g.,  Ilakada  and  'lead,  1(X*5,  which  contains  further 
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references]'  Magnetic  and  electric  disturbances  associated  with 
the  solar  wind  and  other  phenomena  cause  lov-snorgy  protons  (and 
electrons)  to  diffuse  radially  inward  from  the  magnetosphere.  These 
disturbances  act  on  a  time  scale  which  is  long  compared  to  the 
cyclotron  and  bounce  periods  of  the  particles,  and  therefore  violate 
only  the  third  adiabatic  invariant.  Zn  addition,  there  are  processes 
which  act  on  much  shorter  time  scales  to  cause  radial  diffusion:  1)  Pc  1 
micropulsations  can  cause  radial  diffusion  in  an  axially  asymmetric 
magnetosphere,  as  pointed  out  by  Roederer  [1967];  2)  longer-period 
micropulsations  (Pc  2  or  3)  can  resonate  with  proton  bounce  time,  violate 
the  longitudinal  invariant,  and  produce  radial  diffusion  [Cornwall,  i960, 
has  evaluated  the  some  mechanism  for  electrons].  It  is  important  to 
evaluate  the  effect  of  these  last  two  mechanisms  quantitatively,  since 
the  diffusion  coefficients  proposed  by  .lukada  and  Mead,  based  on  long¬ 
time-scale  disturbances,  i3  much  too  small  to  e\plaL  t'.e  observed  ruiU.l 
iitrioution  of  0avi3-VJillianson  protons. 

Probably  the  most  important  effect  of  all  is  the  loss  rate 
induced  by  pitch-angle  diffusion,  wriich  competes  strongly  with  radial 
diffusion.  The  loss  rate  is  nonlinear  vnen  saturation  effects  set  in 
[Cornwall,  1966;  Kennel  and  Petschek,  i960].  When  the  proton  flux  tries 
to  exceeu  a  certain  level,  the  loss  rate  becomes  very  large  and  limits 
the  proton  flux.  W hen  this  loss  rate  is  combined  witn  radial  diffusion, 
a  nonlinear  differential  equation  results. 

In  this  paper,  ve  estimate  the  strength  of  various  radial 
diffusion  processes,  and  combine  this  information  with  nonlinear  pitch-angle 
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diffusion  to  estimate  t.ie  proton  :'l.i>:  nt  various  positions  in  the 
magnetosphere.  It  turns  out  tnat  the  proton  flux  level  at  saturation 
depends  only  logarithmically  on  the  radial  diffusion  constant (  which 
is  fortunate  since  this  diffusion  coefficient  is  very  poorly  known. 
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2.  RADIAL  DIFFUSION  PROCESSES 


The  reader  should  consult  Table  I  for  a  quick  survey  of  the 
results  of  thla  section,  which  is  partly  a  review  of  other  authors'  work* 
Certain  diffusion  processes  are  associated  with  electromagnetic  fields 
which  are  virtually  pure  magnetic ,  and  thus  energy-conserving;  Table  I 
carries  a  column  for  thio  Information. 

A.  Long-tliae-scalt  radial  diffusion  processes  *  These  are 
revieved  by  Dungey  [1965],  and  important  work  has  been  done  by  F&lthammar 
[1965,  1966].  We  define  long  time  as  a  time  long  compared  to  the  bounce 
period  and  cyclotron  period  of  Davis-WiUiamson  protons;  thus  the  primary 
effect  of  these  processes  is  to  violate  the  third  adiabatic  invariant, 
while  the  first  two  are  conserved.  If  these  processes  were  the  only  ones 
acting,  then  the  energy  E  of  an  equatorial  particle  should  be  proportional 
to  the  earth's  field  strength  3~L  ,  by  invariance  of  the  magnetic  moment. 

Many  people  have  discussed  diffusion  by  magnetic  impulses  and 
sudden  commencements,  which  have  a  short  ri3e  time  and  long  decay  time 
[e.g.,  F&lthammar;  tfakada  and  Mead] .  The  diffusion  coefficient  as  given  by 
:.rakada  and  Mead  is  D~3(L/10)^  R^/dey  [Rg  is  the  earth's  radius].  The 
steep  radial  dependence,  as  F&lthammar  has  shown,  comes  from  the  very  long 
decay  time  of  the  magnetic  impulses.  This  has  the  consequence  that,  at 
L  ■  6,  this  type  of  diffusion  is  much  too  small  to  explain  the  Davis- 
Williamson  protons. 

A  potentially  much  more  important  mechanism  is  related  to  mag¬ 
netic  bays,  and  their  associated  electrostatic  potential  fields  [Dungey, 


£> 

F&ltharamar  ];•  Magnetic  oays  la  >t  some  10  to  30  minutes  (  and  are  caused  by 
ionospheric  currents  which  are  driven  by  electrostatic  fields,  Because 
magnetic  field  linos  are  very  nearly  equipotential  lines,  these  ssae 
electrostatic  fields  persist  to  large  radial  distances,  vnere  they  energize 
charged  particles.  The  |ime  scale  (10-30  minutes)  is  short  compared  to 
proton  drift  poriods,  hence  (as  Fdlthsmmar  has  shovn)  the  radial  diffusion 
coefficient  goes  roughly  like  L^,  at  large  (L  t  5)  distances.  For  smdller 
distances,  this  kind  of  diffusion  is  negligible,  because  tne  electrostatic 
potential  vanishes  rapidly  below  subauroral  latitudes. 

F&lthammar's  expression  for  the  diffusion  coefficient  is 


where 


P 

n 


En(L,t)  EQ(L,t+T)  >  cos  nwDT 


(2) 


is  the  power  spectrum  at  the  frequency  nw^T  of  the  azimuthal  electric 
field  E(L,t): 


E(L,t) 


E 


En(L,t)  cos  (nwDt+$n) 


(3) 


In  these  expressions,  B  is  the  earth's  magnetic  field,  is  the  drift 
frequency,  $n  a  phase  angle  presumably  randomly  distributed,  and  the 
quantity  in  angular  brackets  is  the  autocorrelation  function  of  the 
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Fourier  components  of  the  electric  field.  For  fields  generated  b/  mag- 

•  ;V  u 

ne tic  hays,  E  points  west  on  the  night  side,  east  on  the  day  s1£a£  and 
so  the  term  with  n=l  in  equations  (1-3)  is  doubtless  the  most  "•important,  '~,?-y‘p3 
^  i .1  in  any  event  onlv  Keep  the  n— 1  term*  For  the  brackets  in  (2*)  . 

<  >  =  E(L)2  e_2T/T  (*«)  " 


where  T  is  a  characteristic  tine  for  magnetic  bays,  say,  2X10^  sec. 


It  is  then  easy  to  find  D: 


C2L^TE(L)2 


UB*  1  + 
o 


$ 


,  Bq  =  0.31  gauss 


Little  is  known  about  E(L) ,  but  extrapolation  of  ionospheric  potentials 

snows  that  the  associated  potential  over  several  earth  radii  night  be 

-3  -9 

1-10  kv.  Ve  therefore  set  E(L)  =10  -  10  in  cgs  units,  and  get 

6 

1  =  (1  -  103)(^)  H2/day  I 


(This  expression  does  not  include  tne  relative  variation  of  E(L)  with 
and  equation  (o)  should  not  be  used  below  L  =  It  or  5.)  This  kind  of 
diffusion  can  easily  dominate  magnetic  impulse  diffusion.  In  going  from 
(5)  to  (6),  we  nave  set  1  ♦  (u^T /2)2  =  1,  since  T  is  considerably 
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shorter  than  the  drift  period*  In  the  opposite  limit,  Twjj»1,  ve  find 

10  -2 
D  ~  L  ,  since  ~  L  ■ 

B.  Short-time-scale  processes*  We  begin  by  discunsing  very 
briefly  certain  electrostatic  instabilities  which  can  promote  the  fast 
diffusion  of  plasma  across  a  magnetic  field.  (Such  fast  diffusion  is  known 
as  iiohm  diffusion  in  laboratory  plasmas.)  An  example  of  the  type  is  the 
universal  instability  [Kr&ll  and  Rosenbluth,  1963].  In  a  plasma  with  a 
density  gradient,  helical  electrostatic  waves  can  propagate  along  a  mag¬ 
netic  field  line  and  be  amplified  bv  the  instability.  The  growth  rate  of 
the  instability  depends  on  the  ratio  of  plasma  scale  length  to  ion  cyclotron 
radius;  the  smaller  this  ratio,  the  more  severe  the  instability.  For  a 
ratio  of  about  10,  one  expects  growth  times  on  the  order  of  minutes. 
Unfortunately ,  this  ratio  is  liable  to  be  about  10^  in  the  magnetosphere, 
unless  sharp  spikes  of  enhanced  (or  depleted)  field-aligned  plasma  exist. 

But  such  sharp  spikes  would  be  rapidly  dissipated  by  drift  instabilities 
in  the  thermal  plasma.  Our  tentative  conclusion  is  that  this  process  is 
not  liable  to  dominate  radial  diffusion. 

Secondly,  we  consider  processes  that  violate  the  longitudinal 
invariant,  i.e.,  processes  with  cnaracteristic  periods  near  the  proton 
bounce  period  of  about  30  sec  (for  a  75-kev  proton  at  L  *  6).  Pc  2-3 
(5  -  1*5  sec)  fits  ideally.  The  formulas  needed  to  calculate  the  diffusion 
coefficient  have  already  been  giver.  [Cornwall,  1966],  in  an  application 
to  electrons  interacting  with  Pc  1.  The  reader  is  referred  to  this  paper 
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for  details;  the  result  is 


D .  t7) 

where  b  is  the  hydromagnetic  wave  strength,  and  3  is  the  earth's 
field  strength.  For  L  *  4,  we  might  guess  b/D  =  10“2  (b  s  5y  at  L  ■  4), 
and  thus 

n  -  o.k  (fe)'1'2  ih  •  (8) 

Of  course,  if  b/B  changes  with  L  very  much,  formula  (7)  instead  of 
(3)  should  be  used. 

Finally,  Pc  1  hydromagnetic  waves  can  cause  radial  diffusion  in 
an  axially  asymmetric  magnetosphere  [Roedcrer,  1967],  These  waves,  of 
frequency  0.2  -  5  Hz,  are  thought  to  he  generated  by  Davis-Williamson 
protons,  so  in  a  sense  the  protons  are  causing  their  radial  diffusion. 

It  is  difficult  to  estimate  the  diffusion  coefficient  without  extensive 
numerical  analysis,  but  we  can  argue  as  follows:  A  Davis-Williamson  proton 
drifts  around  the  earth  on  a  time  scale  measured  in  hours.  During  this 
tine,  strong  pitch-angle  diffusion  may  change  the  pitcn  angle  by  a  few 
tenths  of  a  radian,  corresponding  to  a  radial  movement  of  a  few  tenths  of 

an  earth  radius.  We  estimate  D  by 

0  -  (9) 

B  -  2  Td 

where  Ar  is  a  fraction  of  an  earth  radius,  and  ?D  the  drift  time  (say 

m-io-q 


n  1 

one  or  two  hours).  An  upper  limit  for  D  night  be  1  R  day"  .  This 
sort  of  diffusion  cannot  occur  below  L  =  5,  where  the  magnetosphere  begins 
to  look  axially  symmetric. 

It  is  clear  from  the  crudity  of  the  estimates  of  this  section 
that  it  is  difficult  to  pin  down  a  good  value  for  the  radial  diffusion 
coefficient.  All  things  considered)  it  is  probably  fair  to  tako  !)  at 
L  ■  6  of  the  order  of  1  /day.  It  would  be  not  at  all  surprising  to 
learn  that  other  workers  have  made  different  (and  better)  estimates  of 
the  diffusion  coefficient,  but  the  author  knows  of  no  published  values  for 
processes  other  than  magnetic  impulses. 
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3.  CO:©IHED  DIFFUSION  AND  NON  LINE  AH  Ia'JSo  RATE 

scattering  of  Davis-Williamson  protons  by  Pc  1  can  be  described 
as  a  diffusion  process  acting  on  the  protons’  pitch  angles,  which  causes  the 
protons  to  leak  into  the  loss  cone.  3ecause  the  protons  themselves  produce 
Pc  1  waves,  the  pitch-angle  diffusion  process  is  nonlinear  [Cornwall,  1966; 
\ennei  and  Detschek],  For  a  sufficiently  strong  source  of  protons,  satura¬ 
tion  occurs  and  the  proton  flux  is  effectively  bounded.  It  would  be  virtu¬ 
ally  impossible  to  combine  the  nonlinear  partial  differential  equations  of 
pitcu— angle  diffusion  and  wave  growth  with  a  radial  diffusion  dif.erential 
equation.  There  is,  however,  a  simple  approximation  which  allows  us  to 
describe  pitch-angle  diffusion  by  an  effective  nonlinear  loss  term,  whose 
magnitude  has  recently  been  calculated  [Cocke  and  Cornwall,  1967  j. 

Let  f(r,v,t)  be  the  Boltzmann  distribution  function  for  the 
protons.  The  time  rate  of  change  of  f  by  pitch-angle  diffusion  can  be 
written: 


df 

dt 


x-W>1/2 


De 


(10) 


where  x  (not  to  be  confused  with  the  spatial  coordinate)  is  the  sine  of 
the  pitch  angle,  D  =  75  ( L/U )“3  sec-1,  and  t  is  essentially  b2/V 
(b  =  Pc  1  magnetic  field,  3  =  earth's  field),  (dee  Cocke  and  Cornwall 
for  details.)  Ne  propose  to  replace  equation  (10)  by  a  much  simpler  form: 


df 

at 


-  x(r)r 


(11) 
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TO  do  30,  we  neeu  to  Know  some  average  behavior  of  f  and  *  .  We  can 
make  an  empirical  fit  to  the  calculations  of  Cocke  and  Cornwall;  this  fit 
is  shown  in  Figure  1.  where  the  data  points  are  computer  calculations,  and 
1  is  inversely  proportional  to  the  omnidirectional  flux  J  of  protons. 

A  good  fit  is: 


«(.*) 


-1* 

1.7  x  10  exp 


-2  X  10 


(12) 


where  J  is  the  omnidirectional  flux  in  cm"2  sec"1.  Tne  number  of  protons 
entering  the  I033  cone  per  second  is 


j  dx  x(l-x2)  f 

at 


-1/2  Qf 

x(l-x  )  D* 


U3) 


x«x. 


where  xq  is  tne  sine  of  the  loss  angle.  Experimentally,  f~x2  (rouSnly), 

so  we  ca rry  out  the  differentiation  in  (13),  multiply  by  a  factor  repre¬ 


senting  the  fraction  of  time  nicropulsations  are  present  (say,  0.3),  «d 
assume  it  is  sufficiently  accurate  to  replace  f  by  J.  The  resulting 


equation  of  the  type  (11)  is: 


dJ 

5t 


a 


(HO 


;cte  t  :e  saturation  effect:  as  T  becomes  large,  tne  loss  rate  grows 
exponentially. 


Ill -10-12 


The  next  3tep  is  to  add  to  (lM  an  effective  source  term  from 
radial  diffusion.  On  very  general  grounds  [Hungey;  Falthammar,  1966]  the 
diffusion  term  looks  like 


±  /nr-2  dJ \ 

dr  \Dr  dr) 


(15) 


(again  we  replace  the  doltzmann  distribution  function  f  by  T  ) ,  where 
r  is  tiie  radius  variable:  r  =  LRe  .  In  equilibrium, 

n  3J  2  3/  -2  3J  \  _w  tv 

0  "  at  ’  r  8?(Dr  a?)  - 


This  is  a  nonlinear  equation  in  J,  which  a  computer  could  easily  handle. 
The  important  nonlinear  feature  means  that,  given  the  spatial  dependence 
of  J,  the  magnitude  of  J  can  be  determined  at  a  particular  point.  To 
illustrate,  suppose  J  ~  r  ,  D  ~  r  ;  equation  (lo)  becomes 

!i(N+3)DJr-2  ■  Jk(J)  (17) 


Take  r  =  6R  ,  D  =  1  R  2  day"^-  ,  N  -  -6  to  find 
e  e 

7  — ?  “1 

J(L=w)  =  n  cm”*'  sec  .  This  is  quite  comparable  to  observed  flaxes, 
dote  that  the  value  of  J  depends  only  logarithmically  on  t  ie  details 
of  the  diffusion  process,  which  is  fortunate  since  these  details  are  very 
poorly  hnovn. 

Of  coarse,  there  are  other  loss  processes  (e.g.,  charge  exchange. 
Coulomb  scattering)  not  considered  here,  but  apparently  in  L  =  U— 6 
pitch-angle  diffusion  is  the  main  1033  mechanism,  on  tue  basis  of  rough 
estimates. 
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k.  CONCLUSIONS 


A  wide  variety  of  physical  effects  contribute  to  radial 
diffusion,  and  these  are  difficult  to  estimate  severally  anu  collectively. 
However,  the  exponential  dependence  of  loss  rate  on  J  minimizes  the 
effect  of  this  uncertainty,  rfe  are  left  with  the  task,  of  calculating 
accurately  the  coefficient  of  in  the  exponent  in  loss  rate  (equar- 

tion  (1?)).  The  work  of  Cocke  and  Cornwall  may  well  contain  systematic 
errors  of  a  factor  of  2  or  3,  which  might  best  be  revealed  by  fitting  the 
observed  data  to  the  solution  of  equation  (lo),  or  a  similar,  more 
ucc irate  equation. 

Little  can  be  said  about  the  radial  dependence  of  J  until 

equation  (l6)  is  solved  numerically.  It  is  unlikely  tn&t  J  can  stray 

-  1* 

too  far  from  a  depcnuence  like  L  ,  since  J  -  L  means  that  the 
exponential  term  in  (lU)  depends  rather  insensitively  on  L.  This  sort 
of  radial  dependence  i3  reasonably  consistent  with  experiment. 

One  major  problem  we  have  left  untouched  is  that  most  of  the 
processes  listed  in  Table  1  give  quite  small  diffusion  coefficients  at 
L  =  3-*»,  with  at  least  the  exception  of  electrostatic  processes,  about 
which  we  know  little.  The  Davis-*filliamson  protons  are  certainly  seen 
in  this  L  range,  so  something  must  be  diffusing  them  inward.  One  is 
tempted  to  speculate  that  universal  instabilities  might  be  important  near 
the  density  knee,  near  L  *  3-1*,  where  plasma  density  gradients  are 
larger  than  elsewhere  in  the  magnetosphere. 
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MEASUREMENT  AND  INTERPRETATION  OF  ULF  AND  VLF  POWER  SPECTRA 


J.  F.  Kenney,  H.  B.  Knaflich  and  H.  B.  Liemohn 
Boeing  Scientific  Research  Laboratories 
Seattle,  Washington  98124 

Recently  a  method  was  proposed  for  estimating  the  pitch  angle  and 
energy  distributions  of  nonthermal  particles  in  the  magnetosphere  based 
on  the  cyclotron-resonance  amplification  of  VLF  and  ULF  whistlers.  The 
method  requires  measurements  of  the  power-spectral  density  of  two  or 
more  consecutive  whistler  elements  which  have  traversed  the  same  path. 
Multi-hop  VLF  whistler  events  recorded  by  a  ground  receiver  (courtesy 
D.  L.  Carpenter)  and  by  satellite  (courtesy  R.  E.  Barrington)  have  been 
analyzed  for  power  spectra  by  playing  a  tape  loop  repeatedly  through  a 
swept -frequency  filter.  Several  multi-hop  ULF  whistlers,  otherwise 
known  as  pearl  micropulsations,  have  been  similarly  analized  for  power 
spectra. 

The  relationship  between  the  observed  power  spectra  of  successive 
whistler  echoes  is  shown  schematically  in  Figure  1.  is  the  initial 
spectrum  incident  on  the  top  of  the  ionosphere,  Pg  and  P^  are  the 
observed  source  and  response  spectra  and  M,  T,  and  R  are  linear  power- 
transfer  coefficients  for  the  magnetosphere  and  ionosphere.  Observed 
ratios  give  the  product  MR,  which  requires  independent  measurements  of 
M  or  R  to  test  the  theory. 
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An  ampligram  of  one  of  the  multihop  VLF  events  (D.  L.  Carpenter) 
which  we  have  analyzed,  and  some  of  its  fine  structure  is  shown  in  Fig¬ 
ure  2.  The  power  spectra  of  elements  in  this  event  are  shown  in  Figure 

3  and  their  successive  ratios  PR/Pg  are  shown  in  Figure  4.  Using  pub¬ 

lished  ionospheric  absorption  curves  to  estimate  R  (Helliwell,  Whis¬ 
tlers  and  Related  Ionospheric  Phenomena)  the  magnetospheric  amplifica¬ 
tion  can  be  estimated  as  shorn  in  Figure  5.  The  cutoff  ratio  of  propa¬ 
gation  frequency  <d  to  minimum  electron/frequency  is  based  on  con¬ 

ventional  dispersion  analysis  (R-mode) . 

The  cyclotron  resonance  interaction  between  VLF  whistlers  and  non- 
thermal  electrons  (0.1-10  keV)  gives  amplification  or  absorption  of 
the  signal  in  the  magnetosphere  depending  on  the  phase  space  distribu¬ 
tion,  F.  The  theory  for  whistler  propagation  in  a  hot  plasma  gives  a 
complex  propagation  vector  k  ■  +  ik^  where  k^  is  given  by  the  cold 

plasma  approximation  and  k^  depends  strongly  on  F.  The  growth  or  decay 

along  the  propagation  path  is  given  by  A(db)  ■  const  /  k^ds  *  10  log^M. 

path 

The  logarithmic  magnetospheric  power  transfer  has  been  evaluated  for 
energy  (E)  and  pitch-angle  (a)  distributions  of  the  form  F  ■  E  nsinma 
(Liemohn,  J.  Geophys.  Res.  12 !,  39-55,  1967)  and  a  typical  result  is 
shown  in  Figure  6.  The  curve  m  ■  0.5,  n  ■  2.5  agrees  fairly  well  with 
the  data  in  Figure  5. 

An  ampligram  of  a  typical  multihop  ULF  whistler  event  (pearl)  is 
shown  in  Figure  7.  The  unsmoothed  power  spectra  of  the  elements  in 
this  event  and  their  successive  ratios  are  shown  in  figures  8  and  9, 
respectively.  There  appears  to  be  a  steady  downward  shift  in  the  peaks 
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of  maximum  power.  Since  reflection  can  only  reduce  the  power  in  suc¬ 
cessive  hops  it  is  evident  that  some  amplification  is  present  along  the 
magnet ospherlc  path.  From  dispersion  analysis  (L-mode)  the  propagation 
path  is  along  L  -  5.3  which  has  a  minimum  proton  gyrofrequency  3.1  Hz. 

The  ionospheric  reflection  coefficient  is  not  sufficiently  well 
known  to  correct  the  observed  ratios.  The  published  reflection  curve  = 
(Field  and  Greifinger,  J.  Geophys.  Res.  70,  4885-4899,  1965)  for  px*  .'f  le¬ 
gation  parallel  to  a  vertical  magnetic  field  (polar  case)  which  st 
shown  in  Figure  10  are  not  quantitatively  applicable  to  our  dal?> 
deep  resonance  spikes  in  |r|  depend  on  the  model  but  may  exp)'.,  r.  r.h* 
large  fluctuations  in  the  power  spectra  that  have  been  deto. d  in 
pearl  elements. 

Typical  theoretical  power-transfer  curves  for  th..*  yclotror.  reso¬ 
nance  interaction  between  ULF  whistlers  and  nonther  protons 
(0.1-10  keV)  in  the  magnetosphere  are  shown  in  FJf  11.  Although  the 
power  transfer  to  Individual  elements  cannot  be.  v.te  mined,  an  estimate 
of  m  =  2,  n« 2  hi  i  been  Inferred  from  certain  ^  -eral  properties  of 
pearls. 

has  been  demonstrated  that  the  product  Miv  can  be  determined  ex¬ 
perimentally  for  both  VLF  and  ULF  whistlers.  For  VLF,  R  is  sufficient¬ 
ly  well  known  to  allow  quantitative  estimates  of  M,  but  for  ULF.  mere 
realistic  estimates  of  R  are  needed  before  an  adequate  measure  cf  M  can 
be  attained. 
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EM  AND  HM  GUIDED  WAVES  IN  THE 


EARTH'S  ATMOSPHERE 

Neil  Brice 

CENTER  FOR  RADIOPHYSICS  AND  SPACE  RESEARCH 
CORNELL  UNIVERSITY 
ITHACA,  NEW  YORK 

In  thia  paper  7  shall  consider  guiding  of  waves  in  the  upper  atmos¬ 
phere  and  then  discuss  the  guided  waves  which  are  of  natural  origin. 
Guiding  of  waves  in  the  earth's  upper  atmosphere  may  be  brought  about 
by  either  inhomogeneity  or  anisotropy.  The  inhomogeneity  may  take  the 
form  of  a  discontinuity  in  refractive  index  (such  as  is  found  at  the 
earth's  surface)  from  which  partial  reflection  occurs  or  of  a  gradient 
in  refractive  index.  In  the  latter  case  in  an  isotropic  medium,  the 

h 

curvature  of  the  ray  is  given  by  the  ■  ■*■  where  n  is  the  refractive 

A 

index  and  k  is  a  unit  vector  in  the  wave  normal  direction. 

The  general  condition  for  trapping  or  ducting  of  waves  along  a 
surface  is  that  the  component  of  refractive  index  normal  to  the  surface 
be  sufficient  to  give  the  ray  a  curvature  at  least  equal  to  the  curva¬ 
ture  of  the  surface.  For  example,  the  ducting  of  high  frequency  vaves 
in  the  magnetosphere  requires  a  gradient  of  refractive  index  in  the 
meridian  plane  normal  to  the  magnetic  field  which  produces  a  curvature 
of  the  ray  at  least  equal  to  the  curvature  of  the  magnetic  field 

Of  greater  interest  to  conjugate  point  studies  is  the  guiding  of 
waves  due  to  anisotropy  in  the  propagation  which  arises  from  the  earth's 
magnetic  field.  In  a  homogeneous  magnetoplasma  such  as  the  upper 
atmosphere,  the  group  ray  (or  vector  group)  velocity  W  of  a  wave 
with  frequency  u>  and  wave  vector  k,  propagating  at  an  angle  e  to 
the  magnetic  field  is  given  by 

.  * 

a  £  dco  .  _©  du: 

V  ’  k  mi  I®  +  1ST  35  |k| 
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Hie  first  ten  gives  the  component  of  group  ray  velocity  in  the 
direction  of  propagation  while  the  second  term,  the  component  normal  to 
the  direction  of  propagation,  arises  from  thr  anisotropy  of  the  medium 
and  produces  the  guiding  to  be  discussed  below.  Hie  angle,  a ,  between 
the  ray  direction  and  wave  normal  direction  is  given  by 

n  *  --1  /-l  dnv  .  -1  /I  dWv 
a  .  tan  (-5  35)  -  tan  g§) 

where  W  is  the  phase  velocity  of  the  wave.  For  guiding  in  the  magnetic 
field  direction,  we  require  a  to  be  negative  and  hence  ~  to  be 

Uv 

positive.  The  magnitude  of  dn/dO  may  be  estimated  from  the  shape  of 
the  phase  velocity  surface  (a  polar  plot  of  phase  velocity,  W,  versus  9). 
If  the  phase  velocity  surface  is  elliptical,  jg  (and  hence  dn/d&)  is 
generally  small,  corresponding  to  a  weak  guiding  condition.  If  the  phase 
velocity  surface  has  a  "figure  8"  shape,  dW/dO  is  generally  large, 
corresponding  to  strong  guiding.  Thun  a  general  idea  of  the  conditions 
under  which  strong  guiding  will  occur  can  be  obtained  from  examination 
of  the  phase  velocity  surfaces  for  a  CMA  diagram.  Figure  1  shows  a  CMA 
constructed  for  conditions  appropriate  to  the  magnetosphere,  i.e.,  a 
plasma  consisting  of  electrons  and  protons. ^  Hie  vertical  scale  is 
the  ratio  of  electron  gyrofrequency  squared  to  wave  frequency  squared 

p  2  2 

r  ■  -At/ of)  and  the  horizontal  scale  is  the  ratio  of  electron  plasma 

frequency  squared  to  wave  frequency  squared  (3^  -  7 r^/cn).  For  each 
area  of  the  diagram,  the  shape  of  the  phase  velocity  surface  is  given  at 
the  right  of  the  figure,  the  vertical  direction  being  that  of  the 
agnetlc  field.  For  strong  guiding  in  the  direction  of  the  magnetic 
field  we  require  a  vertical  "figure  8"  phase  velocity  surface,  which  is 
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found  in  regions  7  and  8  for  the  right-hand  polarized  (whistler)  mode 
and  region  10  for  the  left-hand  polarized  (Alfven)  mode.  Thus  the  strong 
guiding  along  the  field  of  right-hand  polarized  waves  occurs  at  fre¬ 
quencies  below  the  plasma  frequency  (Xg  >l)  and  the  electron  gyro- 
frequency  (Ye>  1)  but  above  the  electron-ion  hybrid  (or  lower  hybrid) 
resonance  frequency,  while  for  the  left-hand  polarized  waves  it  occurs 
at  frequencies  below  the  ion  gyrofrequency. 

For  a  sufficiently  dense  plasma  the  phase  velocity  for  the  right- 
hand  circularly  polarized  mode  (W  )  at  low  frequencies  is  given  by 
(ignoring  ion  effects) 


e 


where 


~e  is  the  electron  plasma  frequency  and 
/"l  is  the  electron  gyrofrequency  and 
is  the  wave  frequency. 


For  the  left-hand,  (L)  ar.d  plasma  (?)  modes  the  phase  velocities 


are  approximately 


and 
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Die  As t roe  equation  for  the  ph.  se  velocity,  H,  as  a  function  of 


propagation  angle,  0, 


(w2-  wf  Kw2-  w?)cos2©  +  (w2-  w?)( w2-  w2)*1®2®  -  0 


then  reduces  to 


(W2-  W^KW2  +  W^)coa2©  +  W4sin2©  «  0  , 
W4(cos2©  +  sin2©)  -  W4  cos2©  ■  0  , 


so  that 


and 


W  -  WR(cos  ©)*  , 


n  -  Dp(cos  ©)"^ 


For  this  case  the  angle  between  the  wave  normal  and  the  ray  is  given  by 


tana  -  -  - 


1  dn 
n  d5 


-  (l/2) tan  © 


The  angle  between  the  ray  and  the  aagentic  field  is 


4 

tan  4 


a  +  0  so  that 

tan  © 

- jp 

2  +  tan  © 


for  aaall  values  of  ©,  4  ■  i  ©  and  for  large  ©,  4  tends  to  xero, 

i.e.,  perfect  guiding.  At  frequencies  well  below  the  ion  cyclotron  fre¬ 
quency, 


W?  -  W2  m  W2 
L  •  x  •  A 


where  V.  is  the  Alfven  velocity  and  the  Astroai  equation  reduces  to 
A 
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■•"'a  '  T  > 


-  /(*.  -  V  )bIvT9  m  0 


•o  that  for  the  alow  mode 


W2  -  W?  cos20 
A 


from  which 


tan  a  ■  i  ||  ■  -tan  0 


ao  that 


-0  and  the  ray  direction  ia  the  direction  of 


the  magnetic  field  (i.e.,  perfect  guiding.) 

The  approximations  uaed  here  break  down  for  the  right-hand  mode 
at  frequencies  below  but  comparable  with  the  electron  gyrofrequency  and 
between  the  ion  gyrofrequency  and  the  electron-ion  hybrid  (or  lower  hy¬ 
brid)  resonance  frequency,  and  for  the  left-hand  mode  at  frequencies 


below  but  comparable  with  the  ion  gyrofrequency. 

In  the  first  of  these  cases  i^),  we  may  use 


W2  - 


c\:<(\  cos© 


to  obtain 


tan  u 


•tan  0 


(1-  <®  sec 


so  that  3  is  negative  and  larger  in  magnitude  than  0  for  small  0, 
and  increases  in  magnitude  more  rapidly  than  0  . 

For  the  right  hand  mode  at  the  ion  gyrofrequency,  we  may  use 
■  0  »  in  the  Astrom  equation  and  obtain 


tan  3  ■ 


tan  0 


2  +  tan  0 


which  provides  reasonable  guiding  only  at  small  angles. 

For  the  left  hand  mode  near  the  ion  gyrofrequency,  we  may  use 


0  ?  V2  «  W2,  V?  «  W2 


to  obtain 


tan  a  ■  - 


tan  0(  14-tan  0) 
2  +  tan2  0 
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It  is  seen  that  in  the  case,  n  t»  -9/2  for  small  0,  and  that 
*>.  ■  -6  for  large  6,  so  that  this  situation  is  very  similar  to  the  low 
frequency  whistler  mode  case  discussed  above. 

We  turn  our  attention  now  from  the  guiding  of  waves  to  the  waves 
which  are  guided  in  the  direction  of  the  magnetic  field  in  the  upper 
atmosphere,  i.e.,  right-hand  polarized  (electromagnetic)  waves  at  fre¬ 
quencies  between  the  electron  and  ion  cyclotron  frequencies  and  left- 
hand  polarised  (hydrcuagnetic)  waves  at  frequencies  below  the  ion  cyclotron 
frequency.  In  the  upper  atmosphere,  the  naturally  occurring  waves  which 
fall  into  the  former  category  are  whistlers  and  so-called  "VU"  emissions 
while  in  the  latter  are  ion  whistlers  and  micropulsations  or  "hydro- 
magnetic"  emissions.  The  range  of  frequencies  encompassed  in  this  paper 
extends  over  7  decades  from  a  lower  limit  of  about  0.1  Hz  (for  which 
the  wavelength  in  the  magnetosphere  is  at  least  several  earth  radii)  to 
about  1.5  MHz  (the  electron  cyclotron  frequency  in  the  ionosphere). 

Since  the  w*  tiers  are  better  understood  then  the  emissions,  I  will  con¬ 
sider  them  first. 

Whistlers  have  their  origin  in  impulses  radiated  by  lightning  dis¬ 
charges,  these  impulses  being  dispersed  by  propagation  in  the  ionosphere 

and  magnetosphere  to  produce  the  whistler.  The  energy  propagates  from 

(2) 

hemisphere  to  hemisphere  along  the  earth's  magnetic  field.  '  Whistlers 
have  many  fascinating  properteis,  a  few  of  which  will  be  mentioned. 

(See  Helllwell^  for  a  more  complete  description).  Whistlers  which  make 
many  successive  traverses  of  the  magnetosphere  produce  echo  trains,  the 
ratio  of  dispersions  of  successive  whistlers  being  2:4:6,  etc.  (if 
initiating  lightning  discharge  was  in  the  same  hemisphere  as  the  ob¬ 
server),  or  1:3:5,  etc.  (if  the  lightning  was  in  the  opposite  hemisphere). 
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Frequently  whistler  energy  is  trapped  in  field-aligned  columns  of 

enhanced  ionization  or  "ducts",  ^  producing  whistlers  with  a  number  of 

discrete  traces,  each  corresponding  to  a  particular  field-aligned  outer- 

magnetoapheric  path.  An  example  of  such  a  multipath  whistler  is  shewn  in 

Figure  2.  Occasionally  the  whistler  follows  different  paths  on  successive 

hops  (mixed-pa  th  whistlers)  and  multiple-flash  lightning  produces  multi- 

flash  whistlers.  Whistlers  have  long  been  used  to  measure  electron 

(2  1)  5) 

densities  in  the  magnetosphere.  *  *  '  More  recently  radial  motions  of 

the  ducts  have  been  measured  using  whistlers,  and  from  these  measurements 

the  magnitude  of  the  electric  fields  in  the  magnetosphere  have  been 

deduced.  ^  An  example  of  movement  of  ducts  is  illustrated  in  Figure  3* 

(7) 

Audio  frequency  receivers  on  board  satellites  have  recorded  proton' 

/a  \ 

and  helium'  '  (ion)  whistlers  in  which  the  dispersion  arises  primarily 
from  the  proximity  of  the  wave  frequency  to  the  appropriate  ion  cyclo¬ 
tron  frequency. ^  An  example  of  a  proton  whistler  is  shown  in  Figure  4. 
The  electron  nr.d  ion  whistler  traces  meet  at  the  crossover  frequency^®^ 
from  which  the  fractional  abundance  of  hydrogen  ions  at  the  satellite 
may  be  determined. The  dispersion  of  the  proton  whistler  is  depen¬ 
dent  on  the  proximity  of  the  wave  frequency  to  the  proton  gyrofrequency 
at  the  satellite  and  is  proportional  to  the  proton  plasma  frequency  at 
the  satellite.  These  properties  of  proton  whistlers  have  been  used  to 
determine  the  earth's  magnetic  field  strength  at  the  satellite  with 
considerable  accuracy ^  and  also  the  proton  number  density. 

The  proton  whistler  is  cut  off  by  cyclotron  damping  at  frequencies  a 
few  cycles  below  the  proton  gyrofrequency.  The  difference  between  cut¬ 
off  frequency  and  gyrofrequency  arises  from  the  thermal  motion  of  the 
ions,  and  has  been  used  to  determine  proton  temperature .  ^ 
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Whistlers  have  been  recorded  in  the  frequency  range  from  a  fev  tena  of 
Ha  to  a  fev  tena  of  kHs. 

Ve  torn  now  to  wavea  which  are  generated  in  the  ionized  atmosphere, 
and  generally  referred  to  as  "emissions". 

(lM 

A  number  of  generalizations  which  may  be  made  about  these  emissions' 
are  illustrated  in  figure  5  •  The  scale  at  the  top  of  this  figure  shows 
the  range  of  frequencies  occupied  by  the  emissions,  from  a  (rather 
arbitrary)  lower  limit  of  about  0*1  Hz  to  an  upper  limit  of  about  10^  Hz. 
The  emissions  observed  on  the  ground  fall  naturally  into  two  calsses, 
those  below  about  10  Hz  which  we  call  hydromagnetic  (HM)  and  those  above 
100  Hs  which  we  call  electromagnetic  (EM). 

At  intermediate  frequencies,  ground-based  observations  of  waves  of 
upper  atmospheric  origin  are  hampered  by  Schumann  resonances  and  power¬ 
line  interference,  and  the  dearth  of  observed  emission  may  be  due  to  this 
interference,  to  the  opacity  of  the  ionosphere  at  these  frequencies  (due 
to  multiple-ion  propagation  effects)  or  to  an  absence  of  noise  genera¬ 
tion  in  this  frequency  range  in  the  magnetosphere. 

Each  of  these  classes  of  emissions  (EM  and  HM)  may  be  further  sub¬ 
divided  into  "discrete"  emissions  which  consist  of  a  series  of  gliding 
(rising  and/or  falling)  tones  and  "diffuse"  emissions,  or  hiss,  which 
most  resemble  band-limited  white  noise. 

As  illustrated  in  Figure  5,  for  both  EM  and  HM  emissions  there  are 

tiro  types  of  diffuse  emission  and  one  type  of  discrete  emission.  The 

diffuse  impulsive  HM  emissions  (P^  micfopulsations)  occur  primarily  at 

night (n)  when  the  energetic  particle  precipitation  contains  both 

electrons  (e)  and  protons  (p).  The  diffuse  continuous  HM  emissions 

(P  micropulsations)  occur  primarily  in  the  morning^)  (M)  when  the 
c 
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energetic  particle  precipitation  consists  primarily  of  electrons  (e). 

f  it) 

Discrete  HM  emissions  (pearls)  occur  in  the  afternoon'  '  (A),  pre- 

f  *  fi  i 

sumably  in  association  with  proton  precipitation.1 

At  night  one  also  finds  a  for  diffuse  EM  emissions  at 

frequencies  generally  above  4  kHz  (auroral  hiss)  While  in  the  morning 
one  finds  diffuse  EM  emissions  below  4  kHz  (polar  chorus  and  mid- 
latitude  hiss)  as  well  as  discrete  EM  emissions  ( chorus ).^19^  The  fre¬ 
quency  ranges  generally  occupied  by  these  emissions  are  also  shorn  in 
Figure  5*  The  night-time  diffuse  emissions  (EM  and  HM)  are  often 
associated  with  auroral  substonns^9'2^ . 

EM  and  HM  discrete  emissions  have  a  number  of  fascinating  character¬ 
istics  in  common.  For  both  types,  a  predominance  of  rising  tones  is 
observed.  (^2,23)  ^  occasion,  periodic  emissions  are  detected,  (^3,24) 
an  example  of  periodic  EM  emissions  being  shown  in  Figure  6.  The 

/nc  q£\ 

emissions  are  observed  alternately  in  opposite  hemispheres '  '*  '  and 

the  period  ox  the  emissions  is  the  two-hop  hemisphere-to-hemlsphere 

(27) 

propagation  delay, '  '  the  delay  being  that  of  left-hand  polarised  waves 

for  HM  emissions and  right-hand  polarized  waves  for  EM  emissions. 

This  property  has  allowed  the  use  of  periodic  emissions  to  determine 
magnetospheric  electron  densities  by  techniques  similar  to  those  used  far 
whistlers/20'29^ 

From  observations  of  periodic  emissions,  it  is  apparent  that  the 
generation  of  one  emission  my  be  stimulated  by  the  arrival  in  the 
generation  region  of  the  echo  of  the  previously  generated  emission, 
and  generation  of  discrete  VLF  emissions  may  also  be  stimulated  by 
whistlers ^  or  man-made  signals  from  VIF  transmitters,^1^  an  exuqple 
being  shown  in  Figure  7*  Multi-phase  periodic  EM  emissions  have  been 
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observed,  i.e.,  evenly  spaced  sets  for  which  the  apparent  period  is  a 
subharmonic  of  the  true  period.  An  example  of  multiphase  emissions  is 
shewn  in  Figure  8.  It  is  found  that  two  sets  or  phases  represents  an 
unstable  configuration  while  three  phases  is  much  more  stable. 

Discrete  EM  emissions  sometimes  show  "frequency  modulation",  i,e.,  a 
quasi-per iodic  change  in  the  center  frequency  of  the  emission,  ^  an 
example  being  shown  in  Figure  9*  The  discrete  emissions  are  generally 
observed  at  frequencies  less  than,  but  comparable  with,  the  minimum 
gyrofrequency  for  the  magnetic  field  line  path  along  which  they 
propagate.  The  generation  mechanism  for  discrete,  emission*:  is  not  well 
understood.  The  emission  certainly  involves  non-linear  interaction  be¬ 
tween  the  waves  and  energetic  charged  particles,  very  likely  in  a  cyclo¬ 
tron-resonance  interaction,  probably  at  several  earth  radii  in  the 
vicinity  of  the  equatorial  plane  and  possibly  involving  trapping  of 
energetic  particles  in  potential  wells  of  coherent  wave  fields.  It  is 
believed  that  discrete  EM  emissions  are  generated  by  energetic  electrons, 
discrete  HM  emissions  by  energetic  protons. 

We  come  now  to  the  generation  of  the  diffuse  emissions,  and  here  we 
find  a  substantial  divergence  between  EM  and  HM  emissions.  Certainly 
sene  and  possibly  all  of  the  diffuse  HM  emissions  discussed  here  are  not 
generated  through  cyclotron-resonance  plasma  instabilities  with  protons, 
whereas  we  believe  that  all  the  diffuse  emissions  discussed  are  gen¬ 
erated  through  cyclotron-resonance  instabilities  by  energetic  elec- 
tron..<33,34) 

The  simultaneity  of  onset  of  the  night-time  micropulsations,  of 
the  cosmic  noise  (riometer)  absorption,  of  the  negative  bay  on  the 
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magnetometer  and  the  overhead  passage  of  the  expanding  auroral  bulge 
strongly  suggests  that  at  least  some  and  possibly  all  of  these  micro- 
pulsations  arise  from  high  frequency  components  of  ionospheric  currents. 
Hie  close  association  between  changes  in  amplitude  of  auroral  hiss, 
auroral  luminosity,  cosmic  noise  absorption  and  the  rate  of  change  of 
vertical  magnetic  field  is  illustrated  in  Figure  10. 

There  is  less  certainty  about  the  morning  P  micropulsations. 

c 

The  range  of  periods  typical  for  these  emissions  is  also  typical  for 
microbursts  (modulation  of  precipitation  of  energetic  electrons)  in 
balloon  X-ray  observations  and  for  the  fading  of  VHF  forward  scatter 
(of  A^  and  S  types),  all  of  which  generally  occur  at  about  the  same  time 
of  day.^)  However,  whether  the  modulation  of  the  magnetic  field  (the 
micropulsations)  causes  the  electron  precipitation  (the  microbursts)  or 
vice  versa  is  not  known.  It  is  probably  significant  that  the  periods 
found  are  of  the  same  order  as  the  mirror  periods  for  trapped  energetic 
electrons . 

We  turn  now  to  consideration  of  diffuse  EM  emissions.  In  Figure  11, 
contours  of  constant  average  amplitude  (in  db)  have  been  plotted  to  show 
the  variation  of  average  amplitude  with  frequency  and  time  of  day  at 
Byrd  Station.  The  main  features  of  interest  are  two  maxi ma,  one  in  the 
late  morning  hours  for  the  low  frequencies  (about  700  Hz)  and  one  at 
night  for  the  high  frequencies  (about  20  kHz.  While  there  is  little 
information  available  between  about  20  kHz  and  150  kHz,  the  ionosonde 
receiver  on  the  Alouette  II  satellite  often  records  intense  noise  at  hi|fr 
latitudes  extending  from  the  lowest  frequency  of  observation  (about 
150  kHz)  up  to  a  frequency  casiparable  with,  but  sosmvhat  less  than  the 
electron  gyrofrequency  at  the  satellite. It  seems  likely  then  that 
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the  night-time  auroral  hiss  extends  up  to  frequencies  of  the  order  of 
1  MHz*  The  higher  frequencies  must  certainly  he  generated  in  or  near  the 
ionosphere.  By  contrast,  the  morning  polar  chorus  is  confined  to 
frequencies  less  than  about  1*3  kHz,  and  is  presumably  generated  in  the 
equatorial  plane  at  several  earth  radii. 

As  ms  mentioned  above,  the  EM  diffuse  emissions  are  probably 
generated  by  energetic  electrons  through  the  cyclotron-resonance  plasma 
instability.  The  growth  rate  Y  of  small-signal  whistler-mode  waves  for 
this  instability  is  given  by^^ 

/\u-%)2  ’iCV  pV  -  e® 

where  is  the  electron  cyclotron  frequency  (e  bAO,  is  the 

ratio  of  wave  frequency  to  JCI ,  rj(vD)  is  proportional  to  the  fraction  of 
electrons  whose  velocity  along  the  magnetic  field  equals  the  resonance 
velocity,  V  and  A(Vp)  is  the  anisotropy  factor  for  the  resonant 
energetic  electrons.  The  form  of  the  equation  is  particularly  interest¬ 
ing,  the  dominant  factors  for  the  growth  rate  being  number  of  resonant 
electrons, ‘n  ,  and  the  anisotropy  factor,  A. 

There  is  not  sufficient  time  to  discuss  this  Instability  in  detail, 
but  so  e  general  observations  may  be  made.  For  stably- trapped  electrons 
in  the  magnetosphere  the  anisotropy  factor,  A,  is  approximately  indepen¬ 
dent  of  the  position  along  a  given  magnetic  field  line.  As  a  result, 
for  a  given  ou  ,  larger  fluxes  of  electrons  (i.e.,  larger  t),  and 
hence  larger  growth  rates)  are  found  in  the  equatorial  plane.  Thus  the 
emissions  generated  by  stably- trapped  electrons  will  be  expected  to  have 
their  origin  in  or  near  the  magnetic  equatorial  plane.  Polar  chorus  and 
aid- latitude  hiss  are  believed  to  be  in  this  category. 
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For  energetic  electrons  newly  injected  into  the  magnetosphere 
(presumably  from  the  tell)  the  anisotropy  factor.  A,  is  larger  in  the 
Ionosphere  than  in  the  equatorial  plane.  Dlls  allows  generation  at  larger 
values  of  where  the  number  of  resonant  electrons,  *1  ,  is  often  very 
much  larger.  Thus  newly  Injected  electrons  are  expected  to  generate  EM 
emissions  primarily  in  the  ionosphere,  where,  it  is  believed,  all  or  most 
of  the  auroral  hiss  is  generated. 

Thus  diffuse  EM  emissions  at  frequencies  less  than  the  electron 
gyrofrequency  in  the  equatorial  plane,  i.e.,  the  polar  chorus  and  mid- 
latitude  hiss,  which  occ^ir  in  the  morning,  are  probably  generated  by 
stably  trapped  electrons.  The  diffuse  EM  emission  at  frequencies  much 
above  the  electron  gyrofrequency  in  the  equatorial  plane,  i.e.,  auroral 
hiss,  which  occurs  at  night,  is  probably  generated  by  newly  injected 
electrons. 

This  work  was  sponsored  in  part  by  the  National  Science  Foundation  , 
Division  of  ; tmospheric  Sciences  under  Grant  NSF  GA-870. 
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Figure  1. 


Figure  2. 


Figure  3* 


Figure  U. 


■  Captions 

A  CMA  DIAGRAM  FOP  ELECTRONS  AND  FROTONS.  The  horizontal  scale 
(Xe)  is  the  ratlo-of  electron  plasma  frequency  squared  to  wave 
frequency  squared  and  the  vertical  scale  (Y^)  the  ratio  of 
electron  cylotron  frequency  to  wave  frequency  squared.  Lines 
indicate  zeros  or  infinites  of  refractive  index  for  the  right- 
handed  (r)  or  left-handed  (L)  waves  propagating  along  the 
magnetic  field  or  the  extraordinary  (X)  or  ordinary  (P)  modes 
propagating  normal  to  the  field.  Sketches  of  the  phase  velocity 
surfaces  (the  magnetic  field  direction  is  vertically  upwards) 
appropriate  to  the  various  areas  of  the  diagram  are  given  to 
the  right  of  th*  figure,  (After  Brice^). 

A  NOSE  WHISTIER  RECORDED  AT  EIGHTS  STATION,  The  zero  on  the 
time  scale  indicates  the  approximate  location  of  the  causative 
atmospheric.  (After  Brice 

THE  CHANGE  CF  LOCATION  OF  WHISTLER  PROPAGATION  PATHS  IN  THE 
MAGNETOSPHERE  ASSOCIATED  WITH  A  POLAR  SUBSTORM.  rJie  whistler 
paths  begin  to  move  in  toward  the  earth  at  about  0620  U.T., 
while  the  Byrd  Station  magnetometer  and  riometer  and  Great 
Whale  River  magnetometer  show  a  sudden  onset  of  activity  at 
O6U5.  (After  Carpenter  and  Stone, ^). 

A  FREQUENCY-TIME  SPECTRUM  OF  A  PROTON  WHISTLER  (below)  AND 
ILLUSTRATIVE  DIAGRAM  (above).  The  proton  whistler  rises  from 
che  crossover  frequency  and  asymptotically  approaches  the  proton 
gyrofrequency  until  cyclotron  damping  cuts  off  the  whistler 
(After  Gurnett  and  Brice^-^). 
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Figure  5- 


Figure  6 


Figure  7 


Figure  8 


Figure  9 


THE  FREQUENCY  RANGE.-  OCCUf  iKD  BY  VARIOUS  TYPES  OF  EM  AND  HM 
EMISSIONS-  The  emissions  art-  also  classified  as  diffuse  (for 
L road-band  emission)  or  discrete  (for  narrow  band),  the  times 
of  maximum  occurrence,  (A  for  afternoon,  M  for  morning  and 
N  for  night)  are  given,  and  che  "ype  of  energetic  particle 
precipitation  (electron  (’e)  or  proton  (p))  with  which  the 
emissions  are  associated.  HM  emissions  sure  also  called 
"micropulsations'1  end  FM  emissions  "VLF  esiissions" . 

.  SIMULTANEOUS  RECOF DINGS  OF  SINGLE-? dASE  PERIODIC  EMISSIONS 

AT  TOO  PAIRS  OF  MAGNETICALLY  CONJUGATE  RECORDING  STATIONS. 

( -?L  1 

(After  Brice v J  ' ) . 

.  RECORDINGS  FROM  GREENBANK  SHOWING  THE  GROUND -WAVE  SIGNAL  FROM 
STATION  N i\h  AT  iU.7  xc,  AND  TRIGGERED  EMISSIONS  RECORDED  ON  THE 
RESEARCH  SHI?  EITANIN.  (After  Helliwell  et  el*31'). 

.  SINGLE- PHASE  AND  SYMMETRICAL  THREE- PHASE  EMISSIONS  RECORDED 
AT  GREAT  WHALE  RIVIK,  3 HI WING  A  PRECISE  THREE- TO-CNE  RELATION¬ 
SHIP  BETWEEN  Tit?  APPARENT  PERIODS  OF  THE  EMISSIONS.  (After 
Brice1 3**h. 

.  QUASI-PEF.IODIC  EMISSIONS  RECORDED  AT  EIGHTS  STATION,  SHOWING 
SLOWLY  RISING  FREQUENCY  BANDS,  EACH  CONSISTING  OF  A  NUMBER  OF 
SETS  OF  PERIODIC  EMISSIONS.  (Recorded  at  Eights  Station,  1930- 
3U  UT  and  19LO-4?  UT  cn  13  Dec,  1961.)  (After  Brice^3*^)- 
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Figure  10.  AVERAGE  AMPLITUDES  OK  VLF  AUR0RA1  HISS,  OJRORAL  5577A  LINE, 
COSMIC  NOISE  AI’.;  'C  \FTION  AND  ULF  Z  (i.e.  the  rate  of  change  of 
vertical  magnetic  field)  RECORDED  AT  THE  ONSET  OF  A  POLAR  SUB- 
STORM  DISTURBANCE  AT  BYRD  STATION,  ANTARCTICA.  (After 
Morozuml  and  Helliwell^1^). 

Figure  11.  AVERAGE  RELATIVE  AMPLITUDE  IN  d.b.  OF  EM  EMISSION  AT  BYRD 

STATION  AS  A  FUNCTION  OF  FREQUENCY  AND  TIME  FOR  THE  MONTH  OF 
JUNE,  1962.  (After  Morozumi  and  Hellivell^*^). 
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Selected  Problems  in  Wave-Particle  Interaction, 
from  Conjugate  Point  Observations 

R.  M.  Gallet 

Institute  for  Telecommunication  Sciences  and  Aeronomy 
Environmental  Science  Services  Administration 
Boulder,  Colorado  80302 

Organized  examples  of  observations  simultaneously  made  at  conjugate 
points,  concerning  both  whistlers  and  some  types  of  VLF  emissions,  will 
serve  as  a  review  of  the  observational  knowledge  of  natural  VLF  phenomena 
in  the  magnetosphere.  This  however  will  be  restricted  to  questions 
touching  the  subject  of  this  symposium,  namely  conjugate  point  studies. 

Of  particular  interest  are  types  of  interactions  between  a  whistler  and 
a  triggered  VLF  emissions  which  could  be  disentangled  only  by  means  of 
conjugate  point  observations. 


III-13-1 


BLANK  PAGE 


b 


POLAR  EMISSIONS 


Leif  Owren 

Radiophysics  Laboratory 
Dartmouth  College 

The  natural  ELF  (0.3  -  3  kHz)  and  VLF  (3  -  30  kHz)  emissions 
occurring  in  the  polar  regions  have  been  observed  at  Vostok, 
Antarctica  and  Thule  (Qanaq),  Greenland  since  1961j..  The  stations 
are  located  one  degree  from  the  South  and  the  North  Geomagnetic 
Poles  respectively.  Recordings  were  started  at  Vostok  (gm.  lat. 
89.1°  S)  in  February  1 961+  and  at  Thule  (gm.  lat.  88.9°  N)  in 
October  1 96I4..  The  Thule-Qanaq  station  lies  about  100  km  north 
of  Thule  Air  Force  Base.  The  Technical  University  of  Denmark 
maintains  ionospheric  stations  including  VLF  equipment  at 
Narssarssuaq  (gm.  lat.  71.2°  N)  and  at  Godhavn  (gm.  lat.  79.9°  N). 
The  chain  of  stations  Narssarssuaq,  Godhavn  and  Thule  spans 
the  polar  cap  region  from  the  auroral  zone  to  the  geomagnetic 
pole  at  about  9°  intervals  of  geomagnetic  latitude  within  a 
1*0°  sector  of  geomagnetic  longitude. 

Initially  the  ELF-VLF  recordings  at  the  geomagnetic 
pole  stations  were  made  with  Stanford-type  IGY  whistler 
recorders  using  the  standard  two-minute  per  hour  sampling 
schedule.  Later  hiss  recorders  were  installed  in  addition 
to  the  tape  recorders.  At  Thule  a  panoramic  receiver  covering 
the  0.5  -  20  kHz  frequency  range  was  put  in  operation 
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in  February  1965*  The  amplitude  vs.  frequency  oscilloscope 
display  is  photographed  every  six  minutes.  At  Vostok  continuous 
hiss  recordings  covering  the  0.2  -  8  kHz  range  were  started 
in  January  1966  using  a  dual-channel  strip-chart  system  which 
displays  the  integrated  rms  amplitude  and  the  spectral  intensity 
distribution  for  each  one-second  frequency  sweep. 

One  problem  arising  when  recording  phenomena  observed 
from  stations  near  the  geomagnetic  poles  is  the  definition  of 
local  magnetic  time.  In  view  of  the  fact  that  auroral  and 
related  ionospheric  phenomena  appear  to  be  centered  on  the 
eccentric  dipole  axis,  we  base  local  magnetic  time  on  the 
eccentric  rather  than  the  earth-centered  dipole.  With  this 
definition  magnetic  noon  at  Vostok  corresponds  to  about  1300  UT 
and  at  Thule  occurs  at  about  1500  UT,  so  while  the  geomagnetic 
longitude  difference  is  about  180  degrees  the  local  magnetic 
time  relative  to  the  eccentric  axis  poles  differs  by  only  two 
hours  for  the  two  stations. 

For  an  understanding  of  the  ELF-VLF  emissions  observed 
in  the  polar  regions  it  is  important  to  take  account  of  the 
propagation  characteristics  of  the  earth- ionosphere  waveguide. 
The  basic  properties  may  be  inferred  from  a  paper  by  Wait  (T957) 
in  the  1957  VLF-issue  of  Proc.  IRE.  According  to  the  results 
of  this  paper  VLF  propagation  over  long  distances  can  take 
place  with  little  attenuation  in  the  fundamental  mode  for 
frequencies  below  2  kHz  and  in  the  first  mode  for  frequenices 
above  8  kHz. 
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In  the  2-8  kHz  band  the  VLF  waves  are  subject  to  strong 
absorption.  Typically,  the  attenuation  is  of  the  order  of  3  db 
per  1000  km  at  1  kHz  and  2  db  per  1000  km  at  10  kHz  but  exceeds 
10  db  per  1000  km  at  1;  kHz.  Thus  the  propagation  conditions 
tend  to  divide  the  ELF-VLF  emissions  observed  in  the  highest 
lati tude' zone  of  the  polar  regions  into  three  categories 
characterized  by  the  three  frequency  bands  0.3-2  kHz,  2-8  kHz 
and  8-20  kHz.  Phenomenologically  there  is  a  natural  division 
between  the  emissions  observed  below  and  above  about  2  kHz  at 
any  polar  cap  station. 

Polar  Chorus 

The  natural  emissions  occuring  between  0.5  and  2  kHz  are 
generally  referred  to  as  polar  chorus  although  in  fact  they 
consist  of  a  mixture  of  chorus-type  emissions  and  hiss.  The 
term  polar  chorus  was  introduced  by  Ungstrup  and  Jackerott  (1963) 
to  distinguish  the  chorus  observed  below  1.5  kHz  at  Godhavn, 
Greenland  during  1957-1961  from  the  middle  latitude  chorus 
occurring  mainly  in  the  2-4  kHz  range.  The  proportion  of 
chorus  and  hiss  in  the  polar  chorus  varies  with  the  sunspot 
cycle,  the  chorus  component  being  more  prominent  near  sunspot 
maximum  and  receding  into  the  hiss  background  with  declining 
solar  activity. 

The  polar  chorus  observed  at  Vostok  during  1964  showed 
the  same  diurnal  and  seasonal  variation  as  that  found  for 
Godhavn  and  other  polar  stations  with  a  diurnal  maximum  at 
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10  -  11  local  magnetic  time  and  a  yearly  maximum  during  the 
local  summer.  The  polar  chorus  at  Vostok  In  1964  consisted 
mainly  of  the  hiss  component.  A  comparison  of  the  occurrence 
and  intensity  of  polar  chorus  at  the  three  Greenland  stations 
Narssarssuaq,  Godhavn  and  Thule  strongly  suggests  that  the 
polar  chorus  is  generated  in  the  auroral  zone  and  propagated 
to  higher  latitudes  in  the  earth-ionosphere  waveguide.  This 
conclusion  was  reached  by  Jackerott  et  al.  (1964)  from  a  detailed 
study  of  the  Narssarssuaq  and  Godhavn  recordings  and  is  further 
strengthened  when  the  Thule  data  are  taken  into  account. 

Auroral  Hiss 

Wide-band  hiss  occurring  in  the  4  “  16  kHz  frequency  range 
and  associated  with  auroral  activity  was  first  observed  at 
Byrd  Station,  Antarctica  by  Martin  et  al.  (I960).  The  hiss 
activity  typically  follows  in  the  wake  of  magnetic  storms  and 
has  a  diurnal  maximum  shortly  before  local  magnetic  midnight. 

The  hiss  observed  at  Vostok  during  1964  at  frequencies  above 
4  kHz  showed  a  similar  diurnal  variation  and  a  similar  relation 
to  magnetic  activity. 

Jorgensen  (1966)  has  studied  the  hiss  events  occurring  at 
Narssarssuaq,  Godhavn  and  Nord  (gm.  lat.  80.8°  N) ,  Greenland 
during  November  and  December  1964  using  continuous  recordings 
at  6  kHz.  He  concludes  that  the  auroral  hiss  is  a  rather  local 
phenomenon.  Drawing  on  observations  in  the  4  -  9  kHz  frequency 
range  from  a  number  of  high  and  middle  latitude  stations, 
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Jorgensen  derived  contour  maps  of  spectral  density  for  wide¬ 
band  hiss  bursts  showing  a  zone  of  maximum  occurrence.  This 
zone  has  the  shape  of  a  horseshoe  with  its  opening  toward  the 
magnetic  noon  meridian.  The  highest  occurrence  rate  and  spectral 
density  of  hiss  events  is  observed  in  the  auroral  zone  at 
about  70°  invariant  latitude  an  hour  before  magnetic  midnight. 

In  the  morning  and  afternoon  hours  there  is  more  hiss  activity 
at  80°  latitude  than  at  70°.  The  auroral  hiss  observed  in 
middle  latitudes  appears  to  be  propagated  from  the  auroral  zone 
in  the  earth-ionosphere  waveguide  as  indicated  by  a  decrease 
in  maximum  spectral  density  by  approximately  10  db  per  1000  km 
and  the  simultaneity  of  the  hiss  bursts  recorded  at  different 
stations. 


Local  Hiss  Events  at  the  Geomagnetic  Poles 

The  continuous  hiss  recordings  at  Vostok  during  1966 
provided  an  opportunity  for  a  systematic  study  of  hiss  events 
observed  near  one  of  the  geomagnetic  poles.  The  0.2  -  8.  kHz 
frequency  range  of  the  hiss  recorder  favors  the  detection  of 
local  hiss  events  over  hiss  propagated  from  the  auroral  zone. 

Two  main  types  of  hiss  bursts  are  observed:  (1)  band-limited 
events  with  a  bandwidth  of  about  3  kHz  and  (2)  wide-band  events 


having  an  essentially  uniform  intensity  over  the  8  kHz  range. 
The  band-limited  bursts  occur  predominantly  in  the  lower  half 
of  the  recorded  frequency  range  but  sometimes  in  the  upper 
half.  Both  types  may  be  represented  during  a  single  burst. 
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Storm  bursts  often  show  a  drift  of  the  maximum  intensity  from 
low  to  high  or  from  high  to  low  frequencies  in  the  course  of 
a  few  minutes.  Spot  checks  of  the  0.3  ■  20  kHz  Vostok  tape 
recordings  coinciding  in  time  with  hiss  events  indicate  that 
the  activity  in  the  2-8  kHz  range  is  the  dominant  feature. 

The  highest  level  of  hiss  activity  tends  to  follow 
moderate  or  moderately  severe  magnetic  disturbances.  Intense 
hiss  storms  lasting  from  6  to  9  hours  occurred  on  27  May, 

1  June  and  31  August  1966  in  the  wake  of  the  magnetic  storms 
of  26  May,  31  May,  and  30  August.  The  severe  magnetic  storm 
of  3  —  4-  September  was  not  followed  by  increased  hiss  activity 
at  Vostok.  This  simply  reflects  the  well-known  migration  of 

the  major  magnetic  disturbances  to  lower  latitudes.  On  the 
other  hand,  an  intense  hiss  storm  lasting  11  hours  occurred  on 
22  July  preceded  only  by  a  weak  magnetic  disturbance  which 
apparently  must  have  been  associated  with  a  polar  storm. 

Short  periods  of  high  hiss  activity  occassional ly  reaching 
hiss  storm  levels,  occur  frequently  in  the  evening  hours, 
part?cularly  during  the  local  winter  season.  These  evening 
"flashes”  between  18  and  24  hours  magnetic  time  typically  last 
1  or  2  hours  but  may  extend  over  periods  of  4.  -  5  hours.  It 
is  the  cumulative  effect  of  these  flashes  which  give  rise  to 
the  diurnal  evening  maximum  in  the  statistics  of  hiss  activity. 

Turning  then  to  the  Thule  six-minute  samplings  of  the 
VLF  hiss,  the  hiss  recorded  during  1965  appeared  mostly  in  the 
6  -  20  kHz  range  and  was  of  the  wide-band  variety.  This 
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situation  changed  in  1966.  The  hiss  storm  activity  now  became 
most  noticeable  in  the  2-8  kHz  range  with  a  significant 
occurrence  of  band-limited  emissions.  The  hiss  activity  at 
the  two  geomagnetic  pole  stations  Vostok  and  Thule  thus  showed 
similar  overall  trends  during  1966. 

Origin  of  Local  Hiss  Events 

It  is  generally  accepted  that  ELF-VLF  emissions  are  caused 
by  an  interaction  of  charged  particle  streams  flowing  along  the 
lines  of  force  of  the  earth's  magnetic  field  with  the  ionosphere, 
although  the  detailed  mechanism  for  the  generation  of  the 
omissions  is  not  understood  (Hel 1 iwel 1 ,  1965).  The  Injun  3 

observations  (Gurnett  and  O'Brien,  1961^)  demonstrated  the 
direct  relationship  between  electron  precipitation,  VLF  hiss 
and  aurora  in  the  magnetic  shell  near  the  bot ndary  of  the  outer 
radiation  zor-  ,  elated  with  the  auroral  zone.  The  occurrence 
of  locally  generated  hiss  in  the  polar  cap  regions  therefore 
suggests  an  occassional  occurrence  of  localized  electron  streams 
in  the  regions  of  the  magnetosphere  beyond  the  boundary  of  the 
outer  radiation  zone.  Evidence  of  such  electron  streams  has 
been  presented  by  Gringauz  and  Khokhlov  (1965)  using  data  from  the 
Russian  satellite  Electron-2  and  the  spacecraft  Mars-1.  These 
data  showed  the  occurrence  of  high  fluxes  of  100  ev  -  L|.0  kev 
electrons  in  high  latitudes.  Similarly,  McDlarmld  and  Burrows 
(1965)  have  reported  from  Alouette  observations  on  the  incidence 
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of  high  fluxes  of  electrons  with  energy  above  I4.0  kev  in  the 
form  of  narrow  spikes  occurring  in  high  latitudes  (70  -  85° 
invariant  latitude)  outside  the  outer  radiation  zone.  Jorgensen 
(loc.  cit.  )  found  a  good  time  correlation  between  the  spikes 
observed  from  Alouette  in  1963  over  the  north  polar  cap  region 
and  the  occurrence  of  strong  hiss  activity  at  Godhavn,  Greenland. 

Magnetic  Storm  Events  During  1965 

The  VLF  emission  activity  at  Thule  and  Vostok  associated 
with  the  sudden  commencement  magnetic  storms  of  1965  has  been 
studied  in  order  to  determine  the  "con jugacy"  aspects  of  the 
VLF  emissions  at  the  two  stations  during  incidences  of  low 
energy  solar  particles  on  the  magnetosphere.  For  each  storm 
and  post-storm  period  the  available  tape  and  .  m  recordings 
were  examined  for  enhancements  in  hiss  or  chorus  activity  over 
the  pre-storm  days.  The  results  are  shown  schematically  in  the 

figure  which  indicates  for  each  magnetic  storm:  (1)  the  time 
and  position  on  the  sun's  disk  of  the  solar  flare  responsible 
for  the  storm,  (2)  the  3-hour  Kp-flgures  for  the  storm  period, 
and  (3)  the  occurrence  and  Intensity  of  VLF  emission  enhance¬ 
ments  at  Thule  (labelled  Kanak)  and  Vostok. 

It  will  be  seen  that  in  ij.  of  the  5  cases  only  one  of  the 
two  stations  gave  evidence  of  increased  VLF  activity.  During 
the  two  June  storms,  designated  as  moderately  severe  and  compris¬ 
ing  the  strongest  magnetic  disturbance  period  of  the  year,  both 
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Thule  and  Vostok  recorded  VLF  emissions.  It  may  be  mentioned 
that  a  comparison  of  the  local  K-figures  for  the  two  stations 
during  each  storm  indicates  the  same  trend  as  the  VLF  emissions. 
The  station  showing  VLF  activity  also  had  relatively  the  higher 
local  K-figures.  The  observations  therefore  suggest  an 
asymmetrical  particle  distribution  and  magnetic  disturbance 
field  over  the  two  polar  regions. 

Consideration  of  the  position  of  the  solar  flare  provides 
no  obvious  clue  to  the  difference  in  response  to  the  solar 
particles  in  the  two  regions.  This  is  perhaps  not  surprising 
since  the  particle  ejection  from  the  sun  is  not  necessarily 
in  a  radial  direction.  Moreover,  the  inclination  of  the  geo¬ 
magnetic  axis  to  the  sun-earth  direction  might  be  expected  to 
be  a  more  important  cause  of  an  asymmetrical  response,  but  the 
observations  at  hand  give  nc  clear  indications  to  this  effect. 
More  detailed  and  comprehensive  information,  including  space¬ 
craft  data,  will  be  required  to  resolve  the  problem. 
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RECENT  RESEARCH  ON  THE  llAGNETOSPHERIC  PLASHAPAUSE 
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Stanford,Callfornla 

ABSTRACT 

The  plasnapause  is  a  three-dimensional  field  aligned  boundary  that 
divides  the  closed  field-line  portion  of  the  earth's  magnetosphere  into 
two  physically  distinct  regions.  The  boundary  is  asymmetric,  usually 
exhibiting  a  minimum  geocentric  range  near  dawn  and  a  maximum  near  dusk 
under  conditions  of  moderate  but  steady  geomagnetic  agitation  (K  >2-4). 

p 

The  mean  equatorial  radius  of  the  plasnapause  is  typically  about  4R_, 
but  it  may  vary  from  about  5.5Rg  during  periods  of  extreme  quiet  to  the 
range  2-3Rg  during  great  storms/  The  approximately  corotatlng  thermal 
plasma  within  the  boundary  exhibits  two  types  of  radial  drift  motions. 
These  may  be  visualised  as:  1)  slow  "breathing"  motions  that  follow  the 
radial  variations  in  a  fixed,  asymmetric  boundary,  2)  more  rapid,  tran¬ 
sient  (1-2  hour)  notions  that  occur  when  the  boundary  position  varies, 
as  is  the  case  during  a  polar  substoim. 

The  asymmetry  of  the  boundary  and  the  nighttime  drifts  of  tubes  of 
ionisation  probably  cause  the  mid-latitude  daytime  and  nighttime  iono¬ 
spheres  to  be  decoupled  from  one  another  insofar  as  the  balance  of  up¬ 
ward  and  downward  fluxes  of  thermal  protons  is  concerned. 
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The  plasmapause  involves  an  abrupt  change  in  electron  density,  in 
tube  electron  content  above  1000  Ion,  and  possibly  in  plasma  bulk  velocity 
and  mean  thermal  energy.  To  the  ionosphere,  the  protonosphere  inside 
the  plasmapause  appears  as  a  large  reservoir  of  thermal  protons,  while 
the  region  outside  appears  virtually  empty.  At  the  plasmapause,  equatorial 
values  of  electron  density  change  by  a  factor  of  10  -  100  within  less  than 
0.15  earth  radii.  Satellite  vlf  experiments  suggest  that  the  change  nay 
be  far  more  abrupt  than  this,  possibly  on  a  scale  of  a  few  kilometers. 
Studies  of  the  distribution  of  electron  density  along  the  field  lines  in 
the  plasmasphere  have  shown  that  earlier  empirical  models  of  the  type 

_3 

N  a  R  are  not  in  fact  compatible  with  recent  satellite  data  on  topside 
electron  concentrations.  Instead,  the  theoretically  palatable  diffusive 
equilibrium  model  has  been  found  to  be  an  appropriate  description  for 
aost  of  the  plasmasphere,  while  a  more  rapidly  varying  model  appears 
necessary  to  describe  the  tenuous  outer  region.  Details  of  the  latter 
distribution  may  vary  in  important  ways  as  a  function  of  local  time. 

Many  wave  propagation  phenomena  of  conjugate  Interest  are  strongly 
affected  by  the  presence  of  the  plasmapause.  For  example,  satellites 
moving  poleward  through  the  boundary  observe:  a  cutoff  in  whistlers 
propagating  from  the  conjugate  hemisphere,  a  decrease  in  the  intensity 
of  fixed-frequency  vlf  signals  propagating  upward,  and  dramatic  changes 
in  vlf  noise  such  as  the  LHR  resonance  noise.  In  ground  recordings  made 
at  Eights  (L  ~  4)  and  Byrd  (L  —  7)  in  the  austral  winter,  four  distinct 
aagnetospherlc  regions  of  propagation  may  be  identified:  (I)  the  outer 
part  of  the  plasmasphere;  (II)  the  outer  "surface"  of  the  plasmapause; 

(III)  s  belt-like  region  extending  1>2I|  outward  from  the  near  vicinity 
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of  the  plasaapause;  (IV)  a  region  beginning  ~  l.SR£  beyond  the  plasaa- 
pauae  and  extending  aeveral  earth  radii  outward.  Each  region  exhibits 
special  properties  with  respect  to  the  occurrence  and  spectral  behavior 
of  vlf  noise,  and  in  particular  noise  triggered  by  whistler  components. 

The  occurrence  of  one-hop  whistlers  propagating  in  the  outer  regions 
III  and  IV  is  relatively  low,  and  is  particularly  low  near  sidnight  and 
in  region  III  at  all  times  except  in  the  aid-afternoon.  Whistlers  propa- 
gatlng  in  regions  I,  III  and  IV  exhibit  an  abrupt  half-gyrofrequency 
upper  cutoff,  while  whistlers  propagating  in  region  II  do  not.  The 
ground  data  on  the  occurrence  of  whistlers  and  noise  are  broadly  consis¬ 
tent  with  recent  surveys  of  magnetospherlc  noise  carried  out  on  satellitea. 
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REVIEW  OF  MAJOR  FEATURES  OF  THE  PLASMAPAUSE 


The  plasmapause  la  a  three-dimensional  field-aligned  boundary  that 
divides  the  closed  field-line  portion  of  the  earth's  magnetosphere  into 
two  physically  distinct  regions.  The  purpose  of  this  paper  is  to  review 
some  of  the  features  of  the  plasmapause  and  to  describe  a  few  of  the  vlf 
wave  phenomena  whose  characteristics  appear  to  be  strongly  influenced  by 
the  boundary. 

The  boundary  during  moderately  disturbed  conditions.  Under  conditions 
of  moderate  but  steady  planetary  geomagnetic  agitation  (Kpe2-4),  the 
plasmapause  exhibits  an  equatorial  radius  similar  to  that  shown  by  the 
dark  curve  in  Figure  1.  The  curve  is  based  on  whistler  observations  near 
the  prime  geomagnetic  meridian  in  July-August,  1963  [Carpenter,  1966a]. 
Important  features  of  the  asymmetric  plasmapause  are  a  broad  minimum  in 
geocentric  radius  near  dawn  and,  near  the  dusk  meridian,  a  bulge  with 
equatorial  radius  1  to  2  earth  radii  larger  than  that  of  the  dawn  minimum. 

Variation  in  the  plasmapause  boundary  with  magnetic  activity.  The 
plasmapause  is  known  to  move  inward  with  increasing  magnetic  agitation, 
but  relatively  little  is  known  about  the  details  of  the  process  by  which 
the  accretion  or  erosion  of  ionization  takes  place.  The  statistics  are 
clear,  however,  in  support  of  a  boundary  that  is  typically  at  4  earth 
radii,  reaches  as  little  as  2  earth  radii  during  great  storms,  and 
extends  to  S  or  more  earth  radii  during  very  quiet  conditions.  Satellite 
measurements  using  a  mass  spectrometer  on  OGO-I  [Taylor,  et  al,  1965] 
and  a  Faraday  cup  on  IMP- 1 1  [Binsak,  1967]  have  provided  excellent 
documentation  of  this  effect. 

An  illustration  of  the  boundary  position  during  4  successive  days 
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is  shown  In  Figure  2,  where  the  plasmapause  radius  is  plotted  on  an 
hourly  basis  in  the  upper  part  of  the  figure  and  the  Kp  index,  with  values 
Increasing  downward,  is  shown  below  (from  [Carpenter,  1966a]).  Following 
a  period  of  quieting  there  is  a  weak  gradual-commencement  magnetic  storm, 
during  which  the  boundary  moves  inward  and  assuaes  the  type  of  diurnal 
behavior  illustrated  in  the  previous  figure.  Some  details  of  the  variation 
of  the  plasmapause  during  periods  of  Increasing  disturbance  and  of  quieting 
are  shown,  respectively,  by  the  dotted  and  dashed  curves  in  Figure  1. 

The  dashed  curve  in  Figure  1  is  a  sisoothed  representation  of  the  data  in 
the  interval  marked  'Q*  in  Figure  2. 

Statistics  on  the  geocentric  radius  of  the  plasmapause  as  a  function 
of  planetary  magnetic  activity  are  shown  in  Figure  3,  where  the  geocentric 
radius  near  dawn  on  48  days  in  July-August  1963  is  plotted  versus  maximum 
Kp  in  the  preceding  24  hours  (from  [Carpenter,  1967]).  For  radii  less 
than  4R  ,  there  appears  to  be  a  grouping  of  the  data  into  2  bands  separated 
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by  about  0.7  earth  radius.  This  appears  to  be  a  real  effect  associated 
with  the  as  yet  unknown  tine  sequence  of  events  during  the  development 
of  a  weak  storm.  During  the  first  2  days  of  a  storm,  points  tend  to  lie 
in  the  right  hand  band,  and  then  appear  to  concentrate  on  the  left  side 
during  the  latter  stages. 

Plasma  motions  and  the  role  of  the  substom.  There  appear  to  be  two 
basic  types  of  drift  motions  of  the  plasma  Inside  the  boundary.  The 
first  type  is  slow  breathing  motions,  during  which  the  boundary  may  be 
thought  of  as  fixed  in  space,  with  the  plasma  Inside  approximately  co- 
rot*  ting  with  the  earth  and  also  drifting  Inward  and  outward  in  conformity 
with  the  asymmetries  of  the  boundary  (with  the  possible  exception  of  the 
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1800-0000  LT  region)  [Carpenter,  1966a| .  The  other  type  of  motion  la  of 
abort  duration  (~  1  hour),  and  appeara  to  Involve  changea  In  the  boundary 
ltaelf.  At  auch  a  time,  the  plasmapauae  and  the  denae  plaama  within  It 
are  observed  to  move  radially  at  apeeds  of  the  order  of  O.AR^/hr,  The 
more  pronounced  Inward  displacements  appear  to  take  place  during  polar 
substorms  and  to  occur  on  the  nightslde  of  the  earth  near  the  longitudes 
where  the  substorm  bay  Is  observed  [Carpenter  and  Stone,  1967],  The 
extent  to  which  the  overall  plasmapauae  boundary  la  modified  la  not  yet 
known,  but  It  appears  that  some  subatorm  activity  occurs  on  nearly  all 
but  the  most  quiet  days.  Due  to  the  rotation  of  the  plasmaaphere,  the 
temporal  sequence  of  compressions  In  the  post-midnight  region  may  lead 
to  longitudinal  variation  in  the  boundary  radius,  and  possibly  to  spatial 
variations  in  ionospheric  parameters  that  depend  on  magnetoapherlc  drift 
motions. 

An  example  of  a  substorm  effect  ia  shown  in  Figure  4  (from  [Carpenter 
and  Stone,  1967]).  In  the  upper  part  of  the  record  are  shown  transcrip¬ 
tions  of  the  magnetic  bay  observed  at  conjugate  stations  Byrd,  Antarctica, 
and  Great  Whale  River,  Canada  (L  ~  7),  with  additional  information  on 
the  Byrd  riometer,  the  magnetometer  at  Eights  (L  ~  4,  one  hour  ahead  of 
Byrd  in  local  time),  and  the  elf  envelope  at  Eights  (0. 2-0.8  Hz).  The 
substorm  was  isolated,  in  that  the  magnetic  records  showed  very  little 
activity  in  the  previous  20  hours.  Planetary  magnetic  conditions  were 
relatively  quiet,  with  near  2.  The  inward  motion  of  the  plasma  is 
demonstrated  by  the  whistler  data  at  the  bottom  of  the  record.  At  the 
left  is  a  scale  for  the  geocentric  radius  of  a  whistler  path,  at  the 
right  a  corresponding  scale  for  measured  whistler  nose  frequency.  A 
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series  of  arrows  shows  various  Important  features  of  the  event,  Including 
an  inward  drift  of  plasma  that  began  prior  to  the  bay.  For  a  period  of 
1-2  hours,  the  plasm  drifted  Inward  at  a  velocity  on  the  order  of 
700  m/sec,  corresponding  to  a  v  X  B  electric  field  of  some  2-3  kV/earth 
radius  at  4  earth  radii  In  the  magnetosphere.  The  plasmapause  boundary 
was  not  Itself  observed  during  the  event,  but  was  subsequently  found  to 
have  shifted  Inward  to  about  4  earth  radii  (further  details  of  the  event 
are  presented  In  the  original  reference). 

Asymmetry  of  the  plasmapause  as  a  function  of  magnetic  activity. 
During  very  quiet  planetary  magnetic  conditions,  the  plasmapause  appears 
to  assume  a  more  nearly  circular  configuration  at  Its  larger  radius. 

This  Is  supported  strongly  by  the  recent  data  of  Blnsak  [1967]  from  a 
Faraday  cup  on  IMP-II.  During  great  storms  the  degree  of  asymmetry 
probably  becomes  more  pronounced,  although  details  are  not  yet  known. 
Corcuff  and  Delaroche  [1964]  reported  knee  effects  between  2  and  3  earth 
radii  during  great  storms,  and  the  present  author  has  observed  depressions 
in  the  magnetosphere  during  the  August  16-19,  1959  storm  which  are  con¬ 
sistent  with  a  knee  at  less  than  2  earth  radii  [Carpenter,  1962]. 

Further  description  of  the  plasmapause  is  divided  into  two  sections, 
one  on  the  changes  in  the  medium  across  the  boundary  and  another  on  vlf 
wave  phenomena  that  are  apparently  affected  by  the  presence  of  the 
plasmapause. 


CHANGES  IN  THE  MEDIUM  AT  THE  BOUNDARY 

At  the  plasmapause,  there  is  an  abrupt  change  in  plasma  density, 
and  possibly  corresponding  changes  in  plasma  bulk  velocity  and  mean 
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thermal  energy.  A  satisfactory  description  of  these  variations  la  not 

yet  available,  although  much  is  already  known  about  the  two  regimes  near 

the  boundary,  If  not  precisely  at  the  boundary  Itself. 

The  change  In  tube  electron  content  at  the  boundary.  At  nighttime, 

inside  the  plasmapause,  the  electron  content  in  a  tube  of  force  extending 

from  a  square  centimeter  at  1000  km  to  the  equatorial  plane  may  be  as 

2 

much  as  an  order  of  magnitude  larger  than  the  total  content  in  a  cm  - 

column  In  the  ionosphere  below.  From  low  latitudes  the  tube  content 

rises  steadily,  reaching  a  maximum  at  the  plasmapause  and  then  falling 

to  extremely  low  values  in  the  tenuous  region  outside.  This  Is  illustrated 

in  Figure  5  from  the  work  of  Angerami  [1966].  The  solid  curve  shows  the 

content  measured  on  a  very  quiet  day,  when  the  boundary  was  presumably 

beyond  the  observing  range  of  Eights  Station.  The  open  and  filled 

circles  show  individual  data  points  from  moderately  disturbed  days 

(Kp=2-4).  It  is  clear  that  the  nighttime  levels  of  2-4.10*3  el/cm^- 

tube  inside  the  plasmapause  provide  at  least  numerically  a  substantial 

12  2 

reservoir  for  an  ionosphere  of  say,  4,10  el/cm  -  column,  while  the 
12  2 

values  of  1-2.10  el/cm  -  tube  for  the  outer  region  are  low  by  comparison. 

The  possible  insufficiency  of  the  protonosphere  as  a  reservoir  for  the 

nighttime  ionosphere  at  middle  to  high  latitudes  may  be  relevant  to  the 

mid-latitude  trough  phenomena  described  by  a  number  of  investigators, 

e.g.  [Muldrew,  1965;  Sharp,  1966;  Stanley,  1966]. 

The  change  in  field-line  distribution  of  ionization  at  the  boundary. 

Until  the  recent  work  of  Angerami  [1966],  whistler  data  were  believed  to 
-3 

support  an  N  «  R  or  'gyrofrequency '  model  of  the  distribution  of 
electrons  along  the  field  lines  of  the  magnetosphere  [Smith,  1961a]. 


Meanwhile,  theoretical  considerations  favored  the  much  more  slowly  decreasing 

_3 

hydrostatic  model  (e.g.  [Johnson,  I960]).  The  convenient  R  model,  still 
useful  in  many  essentially  model-independent  calculations,  was  developed 
using  what  are  now  known  to  be  unrealistic  assumptions  about  the  latitudinal 
distribution  of  electrons  in  the  upper  F-reglon.  Using  extenswe  new 
sources  of  whistler  data  and  satellite  information  on  ionospheric  electron 
density  at  1000  km,  Angeraml  demonstrated  that  the  hydrostatic  model  is 
in  fact  the  best  choice  for  describing  the  region  inside  the  plasmapause. 

In  the  same  investigation,  he  found  that  the  tenuous  outer  region  or 

'trough'  cannot  usually  be  described  by  the  hydrostatic  model,  but  In 
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fact  appears  to  require  a  distribution  varying  roughly  as  R  .  Thus  the 
gyrof requency-type  model  may  again  be  Invoked,  with  the  supporting  physical 
argument  that  the  outer  region  is  descrlbable  by  a  'collisionless'  model 
in  which  particles  move  in  ballistic  orbits  (see  [Evlatar,  et  al,  1964]), 
and  in  which  forbidden,  collision-associated  orbits  become  populated  only 
gradually.  Much  remains  to  be  learned  about  electron  density  in  the 
outer  region,  in  particular  the  manner  in  which  the  field-line  distribution 
varies  with  time  from  after  midnight  to  late  afternoon.  At  night  the 
distribution  may  be  essentially  'collisionless,'  while  by  afternoon  it 
may  approach  a  collision-dominated  model. 

Change  in  the  equatorial  profile  of  electron  density  at  the  boundary. 
Fine  details  of  the  electron-density  profile  at  the  boundary  are  not  yet 
known,  and  perhaps  may  not  be  known  until  high-time-resolution  experiments 
are  flown  on  'satellites.  From  whistler  data,  it  was  found  that  the  electron 
density  decreases  by  between  one  and  two  orders  of  magnitude  within  a 
distance  of  less  than  .15  earth  radius  [Angeraml  and  Carpenter,  1966]. 

ni-15-9 


The  aatellite  data  of  Taylor,  et  al  [1965]  and  .Unsak  [1967]  Indicate 
that  plasna  denalty  at  the  plasmapauae  decreaaea  by  e  factor  greater 
than  10,  but  these  experiments  do  not  provide  fine  details,  say  on  a 
scale  of  10  km.  There  are  some  indications  from  satellite  vlf  experiments 
that  the  boundary  may  at  times  be  only  a  few  km  in  thickness.  Abrupt 
changes  in  whistler  rate  and  vlf  noise  phenomena  are  observed  over 
distances  of  the  order  ci  10  km  on  Alouette  I  at  1000  km  [Carpenter,  et 
al,  1967]  and  over  distances  of  20-100  km  on  OGO-III  at  20,000  km 
(K.  Howell,  private  communication). 

A  summary  of  the  variation  in  electron  density  across  the  'knee' 
is  shown  in  Figure  6,  which  is  based  on  whistler  data  for  June-July  1963 
[Angerami,  1966].  Outside  the  knee  nighttime  densities  (filled  points) 
are  extremely  low,  so  low  that  during  moderately  severe  storms  near 
sunspot  minimum, the  plasma  frequency  in  the  trough  may  fall  below  the 
local  gyrof requency  over  a  large  portion  of  the  path.  As  pointed  out 
by  Bell  [1966]  this  condition  will  permit  the  propagation  of  the  left- 
hand  mode,  and  may  give  rise  to  enhanced  Landau  damping  of  whistler  waves. 

The  density  outside  the  knee  is  sensibly  higher  on  the  afternoon 
side  of  the  earth  than  at  night,  and  appears  to  reflect  the  filling  from 
below  of  the  tubes  of  ionization.  As  noted  above,  the  manner  in  which 
the  tubes  are  filled  has  not  yet  been  studied,  and  could  provide  good 
research  opportunities.  An  example  of  a  problem  would  be  to  check  on 
models  of  high-latitude  convective  motions  in  which  the  daily  cycle  of 
filling  of  tubes  at  high  L  values  is  interrupted  in  specific  ways  (e.g. 
[Nishida,  1966,  Brice,  1967]). 

Attention  will  now  focus  upon  a  number  of  wave  phenomena  that  are 
apparently  affected  by  the  presence  of  the  plasmapause. 
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WAVE  PHENOMENA  AFFECTED  BV  THE  PRESENCE  OF  THE  PLASMAPAUSE 


Satellite  observations  of  a  whistler  cutoff  and  vlf  noise  effect* 
at  the  boundary.  It  has  been  known  for  some  time  that  the  probability 
of  observing  a  whistler  component  propagating  In  the  tenuous  region  out¬ 
side  the  plasmasphere  is  relatively  low  by  comparison  to  the  probability 
of  observation  in  the  dense  Inner  region  [Carpenter,  1963].  A  recent 
statistical  study  of  the  occurrence  of  whistlers  as  a  function  of  observing 
station  latitude  and  magnetic  activity  supports  this  conclusion  [Allcock, 

1966] .  It  has  recently  been  shown  that  the  whistler  rate  observed  on 
satellites  Alouette  I  and  II  drops  abruptly  as  the  vehicle  moves  to  in¬ 
variant  latitudes  beyond  the  plasmapause  [Carpenter,  et  al,  1967]. 
Frequently  associated  with  this  whistler  cutoff  are  abrupt  changes  in 
the  structure  and  intensity  of  vlf  noise,  one  such  change  being  an 
apparent  'breakup'  in  the  lower  hybrid  resonance  (LHR)  noise  band.  A 
summary  of  12  case  studies  comparing  ground  whistler  data  on  the  plasma¬ 
pause  position  with  the  invariant  latitude  of  the  whistler  cutofi  and 
noise  breakup  effect  is  present  in  Figure  7  (from  [Carpenter,  et  al, 

1967] ).  The  cases  are  ordered  according  to  maximum  K  at  Byrd  Station  in 
the  preceding  12  hours  (4th  column),  and  for  the  group  corresponding  to 
K  =  5,  according  to  plasmapause  position.  There  is  clearly  a  close 
spatial  relation  between  the  plasmapause  and  the  whistler  cutoff  and 
noise  breakup  (when  the  latter  is  observed).  The  agreement  persists 
over  a  wide  range  of  variation  in  magnetic  activity  and  a  corresponding 
Invariant  latitude  range  of  52°  -  66°.  The  apparent  tendency  for  the 
whistler  cutoff  and  noise  breakup  to  occur  slightly  equatorward  of  the 
plasmapause  is  not  yet  understood,  but  may  be  attributable  to  uncertainties 
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in  the  model  of  tne  earth's  field  used,  in  the  calculations  (the  plasma- 
pause  position  at  1000  km  was  calculated  in  dipole  coordinates).  In 
particular,  use  of  a  model  in  which  the  field  at  several  earth  radii  is 
dilated  by  plasma  currents  (see  [Mead  and  Cahill,  1966;  Heppner,  et  al, 
1967])  would  tend  to  reduce  the  offset  shown  in  the  figure. 

Reviewing  this  initial  study,  it  appears  that  broadband  vlf  measure¬ 
ments  on  satellites  may  provide  one  of  the  most  precise  methods  of 
determining  the  location  of  the  plasmapause  and  of  identifying  the  range 
over  which  the  transition  from  the  inner  to  the  outer  regime  takes  place 
in  a  physical  sense. 

Ground  observations  of  propagation  from  the  conjugate  hemisphere. 

The  probability  of  observing  a  whistler  trace  propagating  to  a  ground 
station  from  the  conjugate  hemisphere  is  a  complicated  function  of 
latitude,  local  time,  and  magnetic  activity  (e.g.  [Laaspere,  et  al,  1964; 
Yoshlda,  1965;  Allcock,  1966]).  To  partially  clarify  this  picture,  there 
is  offered  in  Figure  8  a  crude  summary  of  whistler  occurrence  in  four 
field-aligned  regions  inside  and  outside  the  plasmapause.  The  figure  is 
based  on  a  visual  survey  of  spectrographic  records  from  Eights,  Antarctica 
(L  ~  4)  for  15  days  in  July-August,  1963.  On  these  days  the  plasmapause 
was  typically  near  L  =  4,  and  planetary  magnetic  agitation  was  moderate 
and  steady  (Kp=2-4),  The  four  regions  of  propagation  are: 

(I)  in  the  plasmasphere,  the  field-aligned  region  within  about 
one  earth  radius  of  the  boundary  in  L  space; 

(II)  the  outer  'surface'  of  the  plasmapause; 

(III)  in  the  plasma  trough,  the  field-aligned  region  extending 
outward  roughly  1  earth  radius  in  L  space  from  the  plasmapause; 
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(IV)  In  the  plasma  trough,  the  field-aligned  region  beginning 
roughly  1  earth  radiua  beyond  the  boundary,  and  extending 
possibly  several  earth  radii  beyond. 

Approximate  local-time  periods  represented  are  2300-0100,  0500-0700,  and 
1400-1600.  For  each  period,  three  2-minute  synoptic  recordings  were 
available.  The  numbers  on  the  graph  represent  the  number  of  periods 
(of  a  possible  15)  during  which  at  least  one  whistler  component  was 
identified  at*  propagating  from  the  conjugate  hemisphere.  Major  features 
of  the  pattern  are:  1)  generally  low  values  beyond  the  plasmapause,  but 
high  ones  in  the  plasmasphere  and  on  the  outer  'surface*  of  the  plasma- 
pause;  2)  extremely  low  values  in  the  plasma  trough  near  local  midnight; 
3)  relatively  low  values  in  a  belt-like  region  outside  the  boundary 
(region  III). 

Propagation  in  region  I  is  frequently  characterized  by  periodic 
emissions,  low  decrement  whistler  echoing,  and  strong  vlf  noise  activity 
near  local  dawn.  Triggering  of  vlf  'chorus'  by  whistlers  is  often 
observed  in  region  II,  and  the  noise  may  persist  for  from  a  fraction  of 
a  second  to  tens  of  seconds  following  the  whistler.  Region  III  appears 
to  be  the  'quietest '  region,  and  its  described  by  whistler  traces  only 
rarely,  the  most  probable  time  of  occurrence  being  mid-afternoon.  Region 
IV,  which  may  possibly  be  associated  with  the  relatively  low-frequency 
'polar  chorus'  band,  is  a  region  of  active  noise  and  frequent  strong 
triggering  of  noise  by  whistlers. 

Certain  spectral  properties  of  whistlers  may  be  identified  with  the 
various  regions.  In  all  but  region  II,  whistlers  observed  on  the  ground 
tend  to  exhibit  a  well-defined  half-gyrofrequency  upper  cutoff  frequency, 
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at  which  the  wave  Intensity  drops  by  10  db  or  more  within  a  fraction  of 
a  kHz  (Carpenter,  1966b].  Such  a  cutoff  appears  to  be  related  to  the 
predictions  of  ray  theory  for  propagation  in  ducts  of  enhanced  ioniza¬ 
tion  (Smith,  1961b].  In  region  II,  the  traces  do  not  usually  exhibit  the 
abrupt  half-gyrofrequency  effect,  but  Instead  show  more  gradual  cutoff 
effects  and  a  wide  range  of  upper  limiting  frequencies,  sometimes  extend¬ 
ing  to  within  10  per  cent  of  the  minimum  gyrofrequency . 

The  occurrence  pattern  in  Figure  8  is  broadly  consistent  with  other 
types  of  observations,  including  the  Alouette  cutoff  effects  described 
above.  The  lack  of  whistler  activity  in  region  III  may  be  related  to 
some  observations  of  Heyborne  (1966],  who  made  a  satellite  study  of  the 
field  strength  of  fixed-frequency  signals  (~18  kHz)  from  vlf  transmitters. 
Digital  data  from  the  polar- orbiting  satellite  OGO-II  showed  that  as  the 
satellite  moved  poleward  near  L  ~  4,  there  occurred  field-strength 
decreases  of  about  25  db  in  less  than  5°  in  daytime  (1600  LT)  and 
about  20  db  in  less  than  10°  at  night  (0400  LT)  .  Heyborne  considered 
it  likely  that  increased  absorption  in  the  lower  ionosphere  accounted 
for  the  effect.  A  problem  in  interpreting  the  whistler  cutoff  in  terms 
of  ionospheric  effects  is  the  observation  on  satellites  of  continued 
upward  propagation  of  short  fractional  hop  (0+)  whistlers  at  latitudes 
poleward  of  the  cutoff  of  whistlers  from  the  conjugate  hemisphere 
(R.  E.  Barrington,  private  communication).  These  persisting  short-path 
whistlers  have  not  yet  been  studied  for  intensity  changes  at  the  plasma- 
pause.  Certainly  much  remains  to  be  done  before  the  propagation 
characteristics  of  the  outer  region  are  well  understood. 
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ILLUSTRATIONS 


Figure  1 


Figure  2 


Figure  3 


Figure  4 


Figure  5. 


.  Equatorial  radius  of  the  plasmapause  versus  local  time.  The 
solid  curve  represents  average  behavior  during  periods  of 
moderate,  steady  geomagnetic  agitation  (Kp=2-4).  The  observa¬ 
tions  were  made  In  July  and  August  1963  at  Eights,  Antarctica. 
The  dots  show  a  particular  example  involving  Increasing 
magnetic  agitation,  the  dashes  an  example  of  decreasing 
agitation  (see  Figure  2)  (from  [Carpenter,  1966a]). 

.  Comparison  of  the  equatorial  radius  of  the  plasmapause  and  the 
Kp  index  during  the  4-day  period  July  28-31,  1963.  The 
quieting,  interval  of  outward  displacement  in  the  plasmapause 
radius,  marked  'Q'  on  July  28-29,  is  shown  on  a  polar  plot  in 
Figure  1  (from  [Carpenter,  1966a]). 

The  geocentric  equatorial  distance  to  the  plqsmapause  near 

local  dawn  versus  maximum  3-hour  K  value  in  the  24  hours 

P 

preceding  the  dlsl  ice  measurement.  The  whistler  data  were 
recorded  in  the  Antarctic  at  Byrd  and  Eights,  and  represent 
a  roughly  30°  range  of  longitudes  near  the  prime  geomagnetic 
meridian.  Centered  dipole  coordinates  were  employed  in  the 
calculations  of  geocentric  distance  (from  [Carpenter,  1967]). 
Some  features  of  a  15  July  1965  polar  substorm  and  associated 
whistler  path  motions.  See  text  and  original  reference  for 
details  (from  [Carpenter  and  Stone,  1967]). 

Scatter  plot  of  day  and  nighttime  data  on  total  tube  electron 
content  above  1000  km  versus  dipole  latitude.  The  whistler 
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observations  were  made  at  Eights,  Antarctica  in  July-August, 
1963,  during  periods  of  moderate,  steady  geomagnetic  agitation 


(K  =2-4).  The  continuous  curve  corresponds  to  K  ~  0  (from 
P  P 

[Angerami,  1966]). 

Figure  6.  Scatter  plot  of  day  and  nighttime  data  on  equatorial  electron 
density  versus  geocentrid'  distance  in  earth  radii.  The  data 
points  correspond  to  those  in  Figure  5  (from  [Angerami,  1966]). 

Figure  7.  Summary  of  case  studies  in  which  the  whistler  cutoff  on 

Alouette  I  is  compared  to  simultaneous  ground  measurements 
of  the  plasmapause  position.  The  K-values  in  the  4th  column 
at  the  left  are  for  Byrd  Station.  See  original  reference  for 
details  (from  [Carpenter,  et  al,  1967]). 

Figure  8.  Crude  statistics  on  the  frequency  of  observation  of  whistler 
components  propagating  to  Eights  Station  from  the  conjugate 
region,  versus  local  tim«  and  path  location  in  the  magneto¬ 
sphere. 
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WHISTLER  PROPAGATION  IN  MAGNETOSPHERIC  DUCTS  * 

by 

R.  L.  Smith  and  J.  J.  Angerami 
Radioed  encc  Laboratory 
Stanford  University 

In  this  paper  we  will  review  some  of  the  theories  of  propagation  of 
whistlers  in  field  aligned  ducts  and  their  relation  to  recent  observations. 

A  number  of  the  features  of  propagation  deduced  from  the  ray  theory  approach 
have  been  verified  by  recent  ground  and  satellite  observations.  These 
observations  in  turn  have  led  to  further  ray  tracing  studies  which  explain 
some  interesting  high  frequency  leakage  effects  and  suggest  an  indirect 
experimental  determination  of  duct  dimensions  and  enhancement  factors. 

Whistlers  observed  by  ground  stations  have  been  explained  as  the  result 
of  the  propagation  of  energy  from  lightning  along  field-aligned  ducts  of 
enhanced  ionization.  The  theory  of  propagation  in  a  smooth  magnetosphere 
without  ducts  was  inadequate  to  explain  many  of  the  features  of  whistlers 
observed  on  the  ground,  such  as  the  multiplicity  of  discrete  components  and 
the  integral  relationships  between  whistler  echoes. 

The  theory  of  propagation  in  ducts  can  be  approached  in  many  ways.  The 
simplest  (and  the  most  productive  so  far)  is  through  the  use  of  ray  tracing 
techniques.  By  using  Poverlein*s  construction  and  by  making  ray  tracing 
calculations  in  a  medium  with  a  uniform  magnetic  field  and  a  variation  of 
electron  density  across  the  field  lines,  one  may  deduce  the  following: 

1)  Whistlers  are  more  readily  trapped  in  ducts  of  enhanced  (rather  than 
reduced)  ionization. 

2)  Frequencies  below  the  lower  hybrid  resonance  may  be  trapped  only 
in  enhancements. 

3)  The  upper  cutoff  frequency  for  trapping  in  enhancements  is  approxi¬ 
mately  one  half  of  the  minimum  gyrofrequency. 
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4)  If  the  electron  density  varies  slowly  across  the  duct,  the  energy 
below  one  half  the  minimum  gyrof requency  can  be  completely  trapped. 

5)  On  the  other  hand,  if  the  duct  has  sharp  sides,  so  that  the 
gradients  of  ionization  across  the  duct  are  quite  large,  then  the  duct 
will  usually  become  "leaky, "  and  some  of  the  energy  will  escape. 

6)  The  average  group  velocity  of  the  trapped  whistler  energy  is  nearly 
the  same  as  that  for  a  whistler  traveling  strictly  along  the  direction  of  the 
static  magnetic  field. 

Another  approach  has  been  used  by  several  workers  who  applied  to  the 
magnetosphere  the  mode  theory  of  propagation  as  developed  for  tropospheric 
ducting.  One  of  the  results  of  this  approach  is  a  trapping  condition  which 
is  a  function  of  the  curvature  and  gradient  of  the  duct.  Another  result 
gives  a  discrete  set  of  trapped  modes  and  a  minimum  size  ("track  width")  of 
the  duct.  We  note  two  points  on  which  some  care  must  be  taken:  1)  the 
anisotropy  of  the  medium  must  be  taken  into  account,  and  2)  one  must  choose 
the  proper  phase  of  the  reflected  wave  relative  to  the  incident  wave. 

Still  another  approach  is  a  full  wave  treatment  using  Maxwell's  equations 
directly.  However,  in  order  to  make  the  problem  tractable,  the  duct  is 
usually  considered  to  have  a  discontinuous  boundary.  Because  of  this,  authors 
following  this  approach  deduce  that  only  depressions  of  ionization  are  suit¬ 
able  for  trapping  the  whistler  energy,  or  that,  if  one  must  study  trapping  in 
the  enhancements,  that  the  ducts  are  very  leaky.  We  see  that  both  results 
depend  directly  on  the  sharp  boundary  assumption  and  hence  do  not  necessarily 
represent  the  physical  case. 

The  ray  theory  of  duct  propagation  has  been  applied  to  ground-based 
whistler  observations  with  considerable  success.  The  most  easily  tested 
prediction  of  the  ray  theory  approach  is  that  the  upper  cutoff  of  whistlers 
observed  on  the  ground  should  be  approximately  one-half  of  the  minimum 
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gyrof roquency  along  the  path.  Carpenter  has  measured  the  upper  cutoff  of 
hundreds  of  nose  whistlers,  and  his  results  show  a  very  strong  tendency  for 
upper  cutoff  frequency  to  equal  one-half  the  minimum  gyrofrequency . 

One  might  ask  what  happens  to  the  whistler  energy  which  is  not  trapped 

in  ducts  of  enhanced  ionization.  The  answer  to  this  question  has  been 

found  in  0G0  satellite  observations  of  a  new  type  of  whistler  which  is 

described  as  "magnetospherically  reflected."  The  wave  normal  of  this  non-ducted 

o 

whistler  rapidly  approaches  90  ,  and  consequently  it  cannot  propagate  through 
the  lower  ionosphere  to  be  observed  on  the  ground. 

Ducted  whistlers  are  also  seen  in  satellites,  mainly  at  L  values 
greater  than  about  3.  At  high  altitudes  (several  thousand  km),  they  exhibit 
frequency -time  spectra  which  are  characterized  by  discrete  changes  in  shape, 
corresponding  to  the  motion  of  the  satellite  through  successive  ducts.  It 
is  therefore  possible  to  identify  the  different  ducts  and  to  measure  their 
size  and  spacings. 

The  best  example  of  these  ducted  whistlers  is  found  in  records  from  an 
inbound,  nearly  equatorial  pass  of  OGO-III  on  15  June  1966.  The  whistlers  are 
classified  as  short  fractional  hop  (0+)  .  We  note  the  following  results: 

1)  In  the  region  between  4.2  and  4.8  earth  radii  (inside  the  plasma- 
pause),  the  satellite  traversed  five  ducts,  their  sizes  ranging  from  200  km 
to  400  km  measured  radially.  The  spacings  between  ducts  range  from  a  few  tens 
of  km  to  1000  km.  Such  density  of  occurrence  of  ducts  is  quite  comparable  to 
the  duct  occurrence  inferred  from  ground  observations  of  whistlers,  under 
similar  conditions  of  magnetic  activity  (very  quiet).  This  has  led  us  to 
suggest  that  the  ducts  are  E-W  elongated,  rather  than  circular. 

2)  The  measured  nose  frequencies  agree  (within  a  measurement  uncertainty 
of  3%)  with  the  predictions  based  on  purely  longitudinal  propagation, 


hydrostatic  distribution  of  ionization  along  the  magnetic  field,  and  the 
measured  intensity  of  magnetic  field  at  the  satellite  (about  6%  lower  than  the 
Jensen  and  Cain  field).  This  comparison  provides  further  support  for  the 
electron  density  model  and  the  longitudinal  approximation  discussed  above. 

3)  Whistlers  are  observed  to  have  an  upper  cutoff  as  high  as  0.75  times 
the  gyrof requency  at  the  satellite  even  inside  a  duct.  This  result,  which 
might  be  thought  to  be  contradictory  to  the  ray  theory  of  ducting,  is 
actually  very  consistent  with  it.  The  frequencies  above  half  the  gyrof requency 
become  untrapped  (as  seen  below)  .  As  the  propagation  changes  from  a  trapped 
mode  to  an  untrapped  mode,  the  energy  nevertheless  is  close  to  the  original 
enhancement  for  some  distance.  This  is  confirmed  by  ray  tracing.  Since 
frequencies  above  hal:’  the  gyrof  requency  are  observed  in  situ  but  not  on 

the  ground,  it  is  suggested  that  the  upper  cutoff  frequency  of  Whistlers  is 
not  due  to  Landau  damping  effects,  as  proposed  by  Scarf  and  Liemohn. 

4)  The  high  frequency  portions  of  whistlers  (lying  roughly  between 

0.5  and  0.75  of  the  local  gyrof  requency)  become  unducted  as  the  gyrof  requency 
decreases  along  the  propagation  path.  They  propagate  then  from  the  outer 
(high  L-value)  ducts  toward  the  inner  region.  We  can  predict  the  behavior 
of  these  unducted  components  from  previous  ray  tracing  experience:  the  wave 
normals  will  tend  to  rotate  away  from  the  magnetic  field  in  an  outward  direction, 
and  the  ray  will  tend  to  rotate  inward.  (Recall  that  the  medium  is  highly 
anisotropic.)  As  a  consequence,  as  the  satellite  moves  inward,  not  only  are 
signals  from  a  given  duct  observed  but  also  the  high-frequency  unducted 
components  escaping  from  ducts  at  higher  L  values. 

The  original  proposals  concerning  the  untrapping  of  whistler  energy 
were  based  on  a  simple  application  of  Poverlein’s  diagram.  A  more  convincing 
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demonstration  of  the  transition  from  the  trapped  to  the  untrapped  mode  can 


be  made  by  tracing  rays  in  a  more  realistic  magnetosphere  including  at  least 
e  dipole  magnetic  field  and  a  small  field-aligned  enhancement  of  ionization 
superimposed  on  a  reasonable  background  ionization  model.  Such  ray  tracings 
have  shown  at  least  qualitative  agreement  with  the  observations.  The 
frequencies  at  which  untrapped  signals  from  an  outer  duct  are  observed  in  an 
inner  duct  may  be  used  to  set  an  upper  limit  for  the  density  enhancement 
factors  in  the  intervening  ducts.  Preliminary  results  based  on  such 
procedure  indicate  an  upper  limit  of  60%. 

♦This  research  was  supported  in  part  by  the  National  Aeronautics  and  Space 
Administration  under  grant  NsG  174-SC/05-020-008  and  contract  NAS  5-2131 
and  in  part  by  the  National  Science  Foundation  through  the  Office  of  Computer 
Sciences  in  the  Mathematical  Division  under  grant  NSF  GP-948. 


FIGURE  CAPTIONS 


Figure  1.  Ray  tracing  of  waves  initially  trapped  in  a  duct  of  enhanced 
(107c)  ionization.  Frequencies  above  half  the  minimum  gyro- 
frequency  (7.55  kHz)  leave  the  duct  (become  untrapped  or  non- 
ducted).  The  points  where  the  local  gyrof requency  are  twice 
the  wave  frequencies  are  indicated  by  the  small  circles. 


Figure  2.  (a)  -  (d).  Ducted  whistlers  and  related  untrapped  components 

(leakages)  as  observed  on  OGO-III.  Figure  2a  shows  whistlers 
observed  while  the  satellite  is  in  duct  no.  1.  Figure  2b  shows 
the  whistlers  observed  in  duct  no.  2.  Figures  2c  and  2d  show 
the  leakage  signals  from  duct  no.  2  observed  in  the  region 
between  ducts  no.  2  and  3.  Note  that  as  the  satellite  moves 
inward  the  observed  low  frequency  cutoff  of  the  leakage  signals 
increases. 

Figure  2.  (e)  -  (h).  Figure  2e  shows  observations  in  duct  no.  3,  where 

the  related  ducted  whistler  is  seen  as  well  as  the  leakages 
from  duct  no.  2.  The  leakage  signals  from  ducts  no.  2  and  3 
are  observed  in  the  interval  between  ducts  no.  3  and  4  as  seen 
in  Figure  2f.  Similarly  wc  observe  in  Figures  2g  and  2h  the 
ducted  whistlers  in  ducts  no.  4  and  5,  as  well  as  the  leakages 
from  the  prior  ducts. 
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WAVEGUIDANCE  IN  THE  MAGNETOSPHERE 
ALONG  FIELD-ALIGNLu  IRREGULARITIES 

Nicholas  Gothard 

Texas  A&M  University,  College  Station,  Texas 


ABSTRACT 

A  study  is  made  of  waveguidance  along  an  inhomogeneous  field- 
aligned  plasma  stratum.  A  time  invariant  model  of  the  ionosphere 
imbedded  in  the  earth's  dipole  field  is  postulated.  Along  the  lines 
of  flux  of  the  earth's  magnetic  field  bell  shaped  plasma  strata  are 
aligned.  The  theory  presented  here  is  based  on  the  mode  of  theory 
of  propagation  in  an  inhomogeneous  medium  and  is  modified  to 
account  for  anisotropy.  The  two  dimensional  problem  is  reduced  to 
a  one  dimensional  integral  equation,  the  mode  condition,  by  using 
Booker's  q  parameter. 

The  results  of  the  solution  are  given  in  terms  of  mode  angles 
and  track  widths. 

The  numerical  results  are  calculated  for  the  locked  modes 
only,  although  the  theory  is  capable  of  handling  the  leaky  modes 
as  well. 

The  numerical  results  indicate  the  existence  of  total  guidance 
between  the  conjugate  points  for  some  frequencies,  for  others  only 
partial  guidance  exist. 
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J.  A.  Thomas 

Physics  Department  (KAAF  Academy) 

University  of  Melbourne,  Australia 

Many  attempts  have  been  made  to  observe  conjugate  ducting  of  HF 
radio  waves  from  the  ground.  Some  long  range  signals  of  more  or  less 
the  correct  radar  range  have  been  detected  by  Obayashi  (1959),  Gallet  and 
"Jtlaut  (1961)  and  Du  Castel  (1965),  but  they  were  not  definitely  established 
as  having  travelled  the  field  line  paths  involved.  Muldrew  (1963)  estab¬ 
lished  the  nature  of  certain  long  delay  echoes  observed  by  Alouette  I 
above  the  ionosphere  -  they  were  attributed  to  multi  hop  ducted  propagation 
between  the  conjugate  reflection  levels  for  the  frequencies  involved.  Such 
echoes  have  since  been  detected  by  the  fixed  frequency  topside  sounder 
Explorer  XX  and  by  Alouette  II.  I  wish  to  briefly  outline  certain  aspects 
of  investigation  into  such  propagation,  with  particular  emphasis  on 
Explorer  XX  observations. 

Loftus,  VanZandt  and  Calvert  (1966)  have  made  an  analysis  of  HF 
conjugate  ducted  echoes  occurring  during  two  3-month  periods  of  19oh  and 
1965,  and  established  the  latitudinal  variation,  diurnal  variation,  KS  duct 
thicknesses  and  separations  and,  (following  Booker  (1962)),  minimum  relative 
enhancements  of  electron  density  required  for  trapping  the  observed  signals. 
Subsequently  (in  conjunction  with  I.T.S.A.),  work  was  commenced  on  follow-up 
problems  related  to  the  E-W  thickness  and  lifetime  of  the  ducts. 

Zz  should  be  pointed  cue  that  observations  of  a  related  phenomenon 
have  given  seme  estimate  of  minimum  E-VJ  duct  widths  for  what  are  believed 
•cc  be  the  lower  ends  of  ducts  which  may  on  occasion  extend  to  the  conjugate 
point.  Duces  may  be  detected  non  only  when  the  satellite  is  directly  within 
the  duct  but  more  often  by  means  of  "combination  mode"  propagation  (Muldrew 
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(1963),  Calvert  and  Schmid  (1963)).  Dyson  (1967)  has  shown  that  field 
aliened  sheet-like  irregularities  with  Ail/N  of  about  2$  and  thickness 
oerper.dicuiar  to  the  field  of  1  km,  can  trap  downward  signals  from  a 
satellite  well  above  the  duct  and  return  a  small  portion  of  the  signal  to 
the  vicinity  of  the  satellite.  Assuming  meridional  propagation,  the 
sequential  observation  on  a  number  of  Alouette  ionogram  frames  of  traces 
due  to  a  duct  on  the  same  L  shell,  gives  an  estimate  of  the  E-W  extent  of 
the  duct  since  the  sub-satellite  path  is  usually  inclined  to  the  magnetic 
meridian.  Typical  maximum  widths  are  20-30  km,  but  these  should  not  be 
regarded  as  maximum  valu.s,  both  because  of  steepness  of  field-lines,  and 
because  the  possibility  of  non-meridional  entry  into  a  duct  should  be 
considered.  It  seems  likely  that  signals  entering  an  elliptical  duct  from 
I  or  1C  of  the  meridian  plane  should  also  be  trapped  under  appropriate  con¬ 
ditions  bun  detailed  ray  tracing  .still  has  to  be  carried  out  for  suitable 
models  and  estimates  made  of  the  power  returned  to  the  vicinity  of  the 
satellite.  This  problem  also  arises  when  such  a  duct  is  at  or  near  the 
satellite  height  -  i.e.,  if  ducts  are  elliptical  in  sections,  how  far  to 
ohe  E  or  \i  of  the  duct  can  it  be  detected  by  the  observation  of  radio 
propagation  to  the  near  or  far  end  reflection  levels? 

Two  experimental  observations  have  been  attempted: 

(a)  Attempts  have  been  made  to  identify  the  same  duct  at  both 
ends  of  a  line  of  force  when  the  satellite  passes  nearly 
along  the  magnetic  meridian  (or  at  least  close  to  the 
conjugate  points). 

and  (b)  Attempts  were  made  to  identify  the  same  duct  observed  by  both 
Explorer  XX  and  Alouette  I  at  times  when  they  passed  through 
the  same  area  close  in  time  or  longitude. 
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If  such  identification  can  he  made  then  valuable  information  is 
available  concerning  the  E-W  extent,  drifting  and/or  lifetime  of  conjugate 
and  "near-end"  ducts.  No  positive  identification  has  yet  been  achieved, 
partly  because  of  time  and  space  difficulties  but  mainly  because  of 
difficulties  associated  with  the  accuracy  of  both  the  satellite  location 
end  magnetic  field  models  together  with  the  methods  of  tracing  field  lines 
to  the  conjugate  point.  Fortunately,  we  are  concerned  with  L  values  of 
not  more  than  2  and  generally  rather  less  than  this,  so  Main  Field  models 
might  be  expected  to  give  a  sufficiently  accurate  picture. 

Some  indication  of  the  problem  can  be  gained  from  the  set  of  observa¬ 
tions  of  Table  1  which  show  the  departures  of  the  magnetic  field  determined 
by  gyrofrecueney  plasma  resonance  of  observations  on  Alouette  I  in  the 
Australian  Area  from  the  values  computed  from  the  GSFC  12/66  set  of 
coefficients  (Cain  (1967)).  The  significant  feature  is  that  for  south- 
going  passes  the  departures  are  negative  whereas  for  north-going  passes  the 
departures  are  positive.  It  is  to  be  hoped  that  the  satellite  locations 
are  not  in  error  by  amounts  sufficient  to  account  for  these  differences. 

In  order  to  establish  the  duct  identity  we  would  like  to  be  able  to  locate 
the  sateliite(s)  to  within  2  or  3  km  and  also  trace  a  field  line  from  the 
satellite  to  its  conjugate  point  (or  at  least  to  the  level  at  which  the 
satellite  passes)  with  the  same  accuracy.  It  is  probable  that  this  is 
beyond  the  presently  available  techniques.  It  may  be  that  the  problem 
should  be  turned  around  and  a  detailed  comparison  of  the  duct  structure 
used  us  establish  a  better  knowledge  of  conjugate  point  locations  in  mid 
and  lev  latitude  regions. 

Some  attempt  has  been  made  to  investigate  the  reason  for  the  dawn 
peak  cf  conjugate  duct  observations  reported  by  Loftus  et  al.  (1966). 
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If  one  examines  the  CRPL  prediction  maps  of  f^F^,  for  1964  and  19° 5  the 
pre-davn  period  is  seen  to  be  one  of  extremely  low  critical  frequencies 
(2.p  to  2.0  K c/s)  and  this  is  followed  by  a  situation  lasting  an  hour  or 
more  in  which  the  critical  frequencies  vary  very  little  with  latitude  but 
steadily  increase  with  local  time.  This  suggests  that  either  (a)  propagation 
along  ducts  at  1.5  -  4  Mc/s  is  more  easily  maintained  under  such  ionospheric 
conditions,  or  (b)  that  ducts  are  more  readily  formed  and  maintained  then. 
Preliminary  tests  with  model  ionospheres  varying  realistically  with  dip 
latitude  and  time  of  day  indicate  that  while  same  latitude  zones  tend  to 
have  more  observable  ducts  crowded  into  them  than  others,  and  these  zones 
change  with  time  of  day,  such  changes  are  not  sufficient  to  account  for  the 
observed  diurnal  variation  of  conjugate  ducting.  The  curves  of  figure  1 
show  the  numbers  of  possible  conjugate  echoes  (0  4  X)  detectable  in  the 
first  four  "windows"  of  Explorer  XX  records  at  2  MHz  assuming  ducts  are 
uniformly  present  at  all  satellite  locations.  The  vast  majority  of  ducted 
echoes  observed  are  propagated  in  the  X-mode  and  all  the  computed  ducts  for 
1-30  Ii-2  are  in  fact  0  mode  (as  the  electron  density  at  the  satellite  is  then 
too  high  to  support  X  mode  propagation).  This  result  is  then  in  agreement 
with  the  observed  midday  minimum  of  conjugate  ducted  propagation  but  the 
midnight  propagation  conditions  are  not  too  dissimilar  to  those  at  0700  US 
and  on  the  basis  of  propagation  changes  alone  one  would  not  expect  a  dawn 
peah  of  occurrence.  The  second  alternative  is  in  many  ways  more  attractive 
since  a  given  fractional  enhancement  of  the  electron  density  is  achieved 
around  dawn  with  a  smaller  absolute  value  of  electron  density  enhancement, 
s.g.,  a  2.5-p  relative  enhancement  at  a  plasma  frequency  of  2.5  MHz  requires 
an  absolute  enhancement  of  2,000  el/cc,  but  at  5  Miz  requires  8,000  el/cc. 
Singleton  (ljc2)  has  c accented  to  some  extent  on  this  feature  in  spread  F  data 
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Finally,  it  is  worth  bringing  to  notice  a  feature  which  is  occasionally 
observed  on  Explorer  XX  records.  At  times  when  the  electron  density  at  the 
satellite  is  Just  below  that  required  for  the  X-mode  propagation  cut  off, 
uhe  virtual  range  of  the  signal  (at  a  fixed  frequency)  is  a  very  sensitive 
indicator  of  the  electron  density  -  greater  ranges  corresponding  to  smaller 
electron  densities  close  to  the  satellite.  Occasionally,  when  such  a 
condition  holds  for  the  1.5  Mc/s  frequency  of  the  sounder,  near-end  ducts 
and  conjugate  ducts  are  observed  on  the  other  frequencies,  and  these  occur 
most  strongly  at  times  of  relative  minimum  electron  density  at  the  satellite. 
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TABLE  I  -  Alouette  I  Gyrofreauency  Measurements 

Comparison  of  GSFC  12/66  Field  Model  with  Field  Values 
Seduced  from  Gyrofrequency  Harmonics  -  Australian  Region 
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Predicted  conjugate  echoes  observable  in 
Explorer  XX  "Windows"  0,  1,  2,  and  3  for 
both  0-  and  X-mode  propagation  at  2  MHz. 
(Satellite  ht.  1000  km)  Ducts  are  assumed 
co  be  available  at  all  satellite  locations 
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MEDIUM-FREQUENCY  CPU  JUGATE  ECHOES  OESERVED  ON 


TOPSIDE-SOUNDER  DATA 


D.B.  Muldrew 
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Defence  Research  Telecommunications  Establishment 

Ottawa,  Canada 


Introduction 

Medium-frequency  conjugate  echoes  (or  ducted  echoes) 
were  first  observed  on  the  Alouette  I  sweep -frequency  topside- 
sounder  ionograms.  An  example  of  an  Alouette  I  ionogram 
containing  con jugate -echo  traces  and  a  sketch  illustrating 
the  propagation  paths  corresponding  to  the  different  traces 
are  given  in  Figure  1.  Combinations  of  the  paths  along  the 
magnetic  field  to  the  north  and  to  the  south  of  the  satellite 
for  the  extraordinary  and  ordinary  waves  are  labelled  using 
the  symbols  N  ,  NQ,  Sx,  SQ.  Conjugate  echoes  were  also  observed 
on  the  Explorer  XX  fixed-frequency  topside -sounder  ionograms; 
an  example,  from  the  paper  by  Loftus,  VanZandt  and  Calvert  (1966), 
showing  two  distinct  ducted  echoes  at  1009:13  and  1009:14  U.T. 
at  apparent  ranges  of  300  and  1200  km  is  shown  in  Figure  2, 

The  time  interval  over  which  each  of  these  echoes  was  received 
indicates  that  the  radio  energy  is  ducted  in  a  region  which  has 
a  thickness,  in  a  direction  normal  to  the  magnetic  field,  of 
about  one  kilometer. 
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An  example  of  an  Alouette  II  ionogram  containing 
conjugate  echoes  is  given  in  Figure  3.  The  Alouette  II  sounder 
is  better  suited  to  the  study  of  conjugate  echoes  than  that  of 
Alouette  I  because  of  the  lower  low-frequency  limit  of  the 
sounder.  The  results  of  a  statistical  study  of  about  60,000 
Alouette  II  equatorial  ionograms,  of  which  about  3,000  (5*) 
contain  conjugate  echoes,  will  be  presented  in  this  paper. 

All  the  available  ionograms  from  eight  equatorial 
stations  were  classified  into  five  groups  depending  on  the 
value  of  f xS,  that  is,  the  frequency  of  the  extraordinary  wave 
which  has  zero  refractive  index  at  the  satellite  height. 

The  five  classes  are:  0.2  (  fxS  (  0.9  MKz,  0.9  <  frS  <  1.25  MHz, 
1.25  <  fxS  <  1.6  MHz,  1.6  <  fxS  <  2.0  MHz  and  fxS  >  2.0  MHz. 

uere  were  no  conjugate  echoes  observed  in  the  last  class  and 
consequently  it  was  excluded  In  the  statistical  study.  The 
fxS  was  chosen  for  the  classifications  because  it  gives  a  rough 
indication  of  the  electron  density  at  the  satellite  height 
and  is  readily  obtained  from  the  ionograms.  Conjugate  echoes 
tend  to  occur  at  frequencies  slightly  above  fxS  so  this  parameter 
also  gives  p.n  indication,  on  the  average,  of  the  frequency  of 
the  conjugate  echoes.  For  each  class  the  recording  time  was 
noted  for  those  ionograms  containing  conjugate  echoes  as  well  as 
for  those  ionograms  not  containing  conjugate  echoes.  The 
percentage  occurrence  of  conjugate  echoes  on  ionograms  v.as  then 
calculated  as  a  function  of  Kp  index,  L  value,  local  time  and 
calendar  date. 
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Tnere  are  two  importrnt  factors  which  must  be 
considered  in  an  Interpretation  of  the  results: 

1.  The  orbital  plane  of  the  satellite  precesses  1.78 
degrees  per  day  with  respect  to  the  sun  so  that  for  north-to- 
south  (or  south-to -north)  passes  of  the  satellite  the  local 
time  at  a  given  latitude  changes  slowly  with  the  date  (day  of 
year).  Thus,  a  plot  of  percentage  occurrence  of  conjugate 
echoes  as  a  function  of  local  time  can  also  be  interpreted  as 
a  plot  of  occurrence  as  a  function  of  date. 

2.  The  argument  of  perigee  rotates  1.89  degrees  per 
day  which  is  almost  the  same  as  the  precession  of  the  orbital 
plane  with  respect  to  the  sun.  Hence,  the  local  time  at 
perigee  (or  apogee)  changes  less  than  3  hours  per  year.  About 
one  year  of  data  is  presented  here  and  for  these  data  perigee 
occurs  near  noon  and  apogee  near  midnight.  Thus  the  satellite 
intersects  the  lov'est  L  shells  near  noon. 

Conjugate  Echoes  as  a  Function  of  Kp  Index 
The  percentage  occurrence  of  conjugate  echoes  on 
ionograms  as  a  function  of  Kp  index  for  all  the  data  with 
0,2  \  fxS  <  2.0  MHz  is  presented  in  Figure  4.  The  error  bars 
are  a  measure  of  confidence  for  the  points;  they  give  the 
standard  error  of  the  mean  assuming,  although  it  is  not  strictly 
true,  that  all  the  data  are  Independent  and  random.  It  can 
be  seen  that  the  occurrence  is  independent  or  decreases  very 
slightly  as  the  Kp  index  increases  from  0  to  6.  Muldrew  (1963) 
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and  Loftus  et  al  (1966)  also  concluded  that  the  occurrence  of 
conjugate  echoes  was  not  correlated  with  Kp  index. 

At  geomagnetic  latitudes  less  than  30°  Calvert  and 
Schmid  (1964),  using  Alouette  J  data,  found  a  weak  negative 
correlation  between  Kp  and  topside  spread  F  (below  the  satellite 
height)  attributed  to  aspect-sensitive  scatter.  This  agrees 
with  the  low  latitude  results  of  Shimazaki  (1959),  Shimazaki 
found  that  the  correlation  between  occurrence  of  spread  F  and 
geomagnetic  activity  is  strongly  negative  for  geomagnetic  latitudes 
less  than  20°  and  strongly  positive  for  latitudes  between  20° 
and  60°.  In  the  present  study  94#  of  the  conjugate  echoes 
occur  at  L  values  greater  than  1.25  and  the  geomagnetic  latitude 
corresponding  to  L  =  1.25  at  300  km  height  is  about  21°.  Thus, 
assuming  L  value  to  be  a  more  appropriate  parameter  than  geomagnetic 
latitude,  the  relationship  between  the  occurrence  of  conjugate 
echoes  and  Kp  does  not  agree  with  the  positive  correlation 
between  spread  F  and  Kp  found  by  Shimazaki  from  ground-based 
stations  at  latitudes  greater  than  20°.  However,  the  characteristics 
of  spread  F  are  dependent  on  the  epoch  of  the  solar  cycle 
(Shimazaki,  1959)  and  Shimazaki's  correlation  results  were  for 
high  sunspot  activity  (IGY)  whereas  the  data  presented  here  is 
for  lower  sunspot  activity  (1966). 

Conjugate  Schoes  as  a  Function  of  L  Value 

The  percentage  occurrence  of  conjugate  echoes  between 
21  hours  and  04  hours  local  time  is  plotted  as  a  function  of  L 
in  Figure  5.  This  period  of  local  time  was  chosen  because  the 
percentage  occurrence  of  conjugate  echoes  remains  fairly  constant 
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throughout  the  period  (see  Figure  6)  and,  as  mentioned  in  the 
introduction,  the  satellite  height  remains  near  apogee. 

To  date,  the  highest  L  value  for  a  conjugate  echo 
observed  with  Alouette  II  is  3.8.  This  is  not  a  real  cut-off 
however  since  the  sweep -frequency  receiver  is  no  longer  tuned 
to  the  frequency  of  the  echoes  that  arrive  after  the  long 
delay  for  propagation  along  a  duct  with  an  L  value  of  about 
4  or  greater. 

The  occurrence  of  conjugate  echoes  increases  with 
decreasing  L  down  to  about  L  =  1.4.  There  is  an  apparent  rapid 
drop  in  occurrence  for  L  less  than  about  1.35.  It  is  not 
certain  at  present  whether  this  drop  is  real  since,  due  to 
the  characteristics  of  the  satellite  orbit,  the  lata  at  these 
lower  L  values  may  not  be  representative.  If  it  is  real  the 
distribution  is  not  unlike  the  distribution  of  energetic 
electrons  (energy  >  1  meV)  with  L  value  (Katz,  1966). 

In  Figure  5  plots  for  0.2  ^  fxS  \  0.9  MHz  and  for 
0.9  <  fxS  <  1.25  MHz  are  given;  plots  for  other  classifications 
are  not  given  because  of  insufficient  data.  For  L  less  than 
2.5  the  occurrence  of  conjugate  echoes  for  0.9  <  fxS  <  1.25  MKz 
is  about  twice  the  occurrence  for  0.2  \  fxS  <(  0.9  MHz.  This 
could  be  due  to  the  lower  efficiency  of  the  antennas  at  lower 
frequencies;  however,  because  of  the  high  echo  signal  strength 
of  conjugate  echoes  this  seems  unlikely.  For  ducting,  the 
transverse  thickness  of  field  aligned  ionization  irregularities 
would  likely  need  to  be  in  the  order  of  a  wavelength  or  more, 
thus,  the  two  distributions  in  Figure  5  could  be  explained  if 
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for  L  values  between  about  1.3  and  2.5,  thin  field-aligned 
ionization  irregularities  of  thickness  less  than  about  300  meters 
were  more  numerous  than  those  of  thickness  greater  than  about 
300  meters. 

Conjugate  Echoes  as  a  Function  of  Local  Time  (Two -Hour  Averages) 

The  percentage  occurrence  of  conjugate  echoes  for 
0.2  <  fxS  <  1.6  MHz  (solid  line)  and  for  1.6  <  fxS  <  2.0  MHz 
(dashed  line)  as  a  function  of  local  time  is  given  in  Figure  6. 

The  data  in  the  three  classifications  between  0.2  <(  fxS  <  1.6  MHz 
were  similar  and  for  simplicity  were  included  together.  The 
data  were  averaged  in  two  hour  intervals  of  local  time.  It 
can  be  seen  that  the  occurrence  of  conjugate  echoes  is  less  than 
about  1$  from  about  lo  to  16  hours  local  time,  starts  to  increase 
before  sunset  and  reaches  a  maximum  of  about  ll«g  just  before 
sunrise.  In  this  plot  the  noon  values  are  for  heights  near  perigee 
(500  km)  and  the  midnight  values  are  for  heights  near  apogee 
(3000  km).  The  Alouette  I  analysis  (Muldrew,  1963)  Indicated 
maxima  at  both  midnight  and  noon;  the  disagreement  at  noon  might 
be  due  to  the  lower  height  of  Alouette  II  for  noon  data  (500 
kilometers  vs  1000  kilometers). 

The  occurrence  distribution  of  conjugate  echoes  for 
the  case  1,6  (  fx2  \  2.0  MHz  is  quite  different  than  that  for 
0.2  <  fxS  <  1.6  KHz.  Although  there  are  not  many  data  during 
the  night  with  1.6  i  fx3  (  2.C  MHz,  it  appears  that  the  occurrence 
of  conjugate  echoes  begins  to  increase  sometime  after  21  hours 
local  tine,  reaches  a  peak  somewhere  near  07  hours  and  drops 
off  to  near  zero  at  13  hours.  This  type  of  distribution  is 
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similar  to  that  obtained  by  Loftus  et  al  (1966)  from  the  fixed- 
frequency  topside-sounder  data.  They  obtained  conjugate  echoes 
on  1.50,  2.00,  2.85  and  3.72  MHz.  There  thus  appear  to  be  two 
quite  different  distributions  for  the  occurrence  of  conjugate 
echoes  as  a  function  of  local  time;  one  for  frequencies  below, 
and  one  for  frequencies  above,  about  1.5  MHz. 

Conjugate  Echoes  as  a  Function  of  Local  Time 
or  Date  (Half-Hour  Averages) 

The  percentage  of  conjugate  echoes  for  0.2  <  fxS  <  1.6  MHz 
as  a  function  of  local  time  is  given  in  Figure  7  (north- to -south 
passes  of  the  satellite)  and  Figure  8  ( sou th-to -north  passes). 

The  class  for  1.6  <  fxS  <  2.0  MHz  was  excluded  because  there 
are  too  few  conjugate  echoes  in  this  class  to  make  a  significant 
contribution  to  the  graphs.  The  data  for  Figures  7  and  8  were 
averaged  in  half-hour  intervals  of  local  time.  It  was  mentioned 
in  the  introduction  that  a  plot  of  this  type  can  also  be 
interpreted  as  a  plot  of  occurrence  as  a  function  of  date;  the 
day  of  the  year  is  thus  given  along  the  top  of  the  figures, 
increasing  from  right  to  left. 

In  Figures  7  and  8  it  is  seen  that  the  occurrence  of 
conjugate  echoes  with  date  or  local  time  does  not  vary  smoothly 
but  tends  to  have  a  considerable  fluctuation.  An  attempt  was 
made,  using  all  of  the  data  of  Figures  7  and  8,  to  determine 
a  regular  period  for  the  fluctuations.  The  best  fit  gives  a 
period  of  1.75  hours  in  local  time  or  14.8  days  in  date.  This 
periodicity  is  illustrated  in  Figures  7  and  8  by  the  arrows. 
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It  can  be  seen  that,  in  general,  the  peaks  in  the  occurrence 
distribution  lie  near  the  arrows.  There  are  some  exceptions 
however,  notably  the  large  maximum  in  Figure  7  at  2130  hours 
local  time  where  a  minimum  might  have  been  expected. 

If  it  is  assumed  that  the  percentage  occurrence  of 
conjugate  echoes  oscillates  in  local  time,  then  this  can  be 
interpreted  as  shown  in  Figure  9.  The  shaded  regions  are  fixed 
relative  to  the  sun  as  the  earth  rotates  beneath  them  and  are 
separated  by  1.75  hours  of  local  time.  The  shaded  areas  represent 
regions  in  the  earth's  magnetosphere  which  are  most  favourable 
for  the  formation  of  conjugate  ducts. 

A  different  interpretation  results  if  it  is  assumed 
that  the  periodicity  of  the  percentage  occurrence  of  conjugate 
echoes  is  one  in  date.  Twice  the  period  of  14.8  days  is  very 
nearly  the  synodic  period  of  the  moon  (29.53  days).  For  the 
data  presented  in  Figures  7  and  8  it  appears  that  conjugate 
ducting  is  most  likely  to  occur  when  the  position  of  the  moon, 
relaLive  to  the  sun  and  earth,  is  in  either  of  the  two  positions 
shown  in  Figure  10.  This  is  not  a  simple  tidal  effect  since 
the  location  of  Alouette  II  is  not  of  primary  importance  in 
determining  when  the  occurrence  of  conjugate  echoes  is  a  maximum. 

Further  studies,  using  both  Alouette  I  and  Alouette  II 
data  should  determine  whether  or  not  the  occurrence  of  conjugate 
echoes  oscillates  with  a  regular  period  and  whether  the  oscillation 
is  with  respect  to  local  time  or  with  respect  to  date. 
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Summary 

The  percentage  occurrence  of  conjugate  echoes  on 
Alouette  II  data  from  November  1965  to  about  October  1966: 

1.  is  either  independent  of  Kp  or  decreases  slightly 
with  increasing  Kp  for  Kp  between  0  and  6, 

2.  increases  for  night-time  data,  from  about  1$  at 

L  =  3  to  18#  (0.2  <  fxS  <  0.9  MHz)  or  32#  (0.9  <  fxS  <  1.25  MHz) 
near  L  ®  1.35  and  then  may  decrease  rapidly  for  L  less  than 
this, 

3.  increases  from  a  minimum  of  less  than  1#  near  noon 
(possibly  affected  by  low  satellite  heights)  to  a  maximum 
of  about  11#  before  sunrise  for  0.2  <  fxS  <  1.6  MHz  and  for 
1.6  <  fxS  <  2.0  MHz  increases  from  very  low  values  between 
13  and  21  hours  local  time  to  a  maximum  near  07  hours, 

4.  appears  to  have  an  oscillation  in  either  local  time  or 
calendar  date  with  periods  of  either  1.75  hours  or  14.8  days. 
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Figure  5  distribution  of  conjugate  echoes  with  L  value  for 
21  to  04  hours  local  tine  Use  Figi  re  6). 


Figure  r  Distributions  of  conjugate  echoes  with  local  tine  when  averaging  Is  done 
over  relatively  long  Intervals  of  local  tlae  (two  hours).  All  data  are 
grouped  Into  a  single  24  hour  local -tine  period. 
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PROPAGATION  C  HA  RAC  TERM  TIC  8  OF  THE  EQUATORIAL  I0N08PHERE  A8 
MEASURED  BETWEEN  MAGNETICALLY  CONJUGATE  POINTS  AT  VHF 

by  E.  D.  Bo  woo 
Smyth  Research  Associates 
3655  Aero  Court 
San  Diego,  California  92123 

A  transeqoatorlal  propagation  experiment  was  conducted  In  the  Far  East  In 
whloh  an  anomalous  propagation  mode  appeared.  This  anomalous  mode  has  been 
Interpreted  as  a  field  guided  mode  as  suggested  by  the  experimental  evidence. 

The  program  for  studying  the  propagation  characteristics  of  the  equatorial 
Ionosphere  was  conducted  from  December  1956  to  June  1968.  In  the  experiment, 
transmissions  at  34,  46,  64  and  77  MHz  from  Okinawa  were  monitored  at  receiving 
sites  at  Darwin,  Australia  and  at  FtJL  The  Okinawa-Darwin  terminals  are  very 
nearly  magnetically  conjugate  and  separated  by  4300  km.  They  occupy  positions 
on  the  safe*  L  shell  (L  »  1. 10),  but  are  separated  by  about  two  degrees  In 
latitude.  The  Okinawa- FIJI  path  of  7200  km  is  skewed  to  about  40*  with  die 
magnetic  equator.  Horizontally  polarized  yagl  antennas  were  placed  three  wave¬ 
lengths  above  a  smooth  foreground  at  each  experimental  site  and  at  eaoh  frequency. 
Phase  locked  receivers  with  narrow  tracking  loops  received  the  emmlsslons  from 
the  stable  frequency,  300  Watt,  CW  transmitters.  These  transmissions  were  on 
occasion  pulsed  at  a  25  cps  rate  and  with  a  100  microsecond  pulse  width  for  travel 
time  measurements. 

Figure  1  shows  the  signal  levels  of  received  frequencies  at  Darwin  (upper 
half  of  figure)  and  at  FIJI  (lower  half)  for  a  characteristic  day  near  the  equinox. 
Thirty-four  MHz  and  occasionally  45  MHz  occur  throughout  much  of  the  day  at 
both  FIJI  and  Darwin,  but  a  propagation  mode  with  markedly  different  character¬ 
istics  occurs  at  Darwin  in  the  late  evening  hours  (2200).  The  Incidence  of  this 
mode  Is  displayed  by  a  significant  Increase  in  frequency  support,  with  the  lowest 
frequency  appearing  first  and  the  higher  frequencies  progressively. 

The  diurnal  occuranoe  of  the  guided  mode  Is  limited  to  the  nighttime 
hours  with  a  maximum  at  2300  local  path  midpoint  line.  The  seasonal  variation 
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maximizes  at  the  equinox,  but  the  biannual  peaks  In  activity  are  broadened  to 
include  the  magnetic  equinox  whloh  follows  the  vernal  equinox  by  fourteen  days. 
This  skewing  of  the  seasonal  dependence  carries  over  to  the  solstlclal  periods 
as  a  higher  activity  at  northern  summer  than  at  northern  winter. 

Nearly  free  space  signal  levels  (*=  -70  dBm)  are  produced  by  the  guided 
mode  on  occasion.  Generally  a  frequency  independence  for  path  loss  is  shown 
that  persists  up  to  some  cut  off  frequency.  This  cut  off  frequency  varies  with 
each  occurance  of  the  mode,  and  severe  path  loss  is  noted  beyond.  An  average 
of  all  the  guided  mode  occurances  yields  a  \  dependence  for  the  path  loss. 

Figure  2  shows  three  fading  distributions  of  the  received  signal  at 
Darwin.  The  lower  histogram  is  the  fading  distribution  obtained  from  a  sample 
of  the  normally  propagated  F  layer  modes.  The  fading  frequency  Is  0.  26  Hz. 

By  contrast  the  fading  distributions  shown  in  the  upper  two  histograms  from  the 
guided  modes  have  more  than  an  order  of  magnitude  higher  fading  frequency. 

The  fading  frequency  was  found  to  be  transmission  frequency  Independent  and 
near  five  Hertz. 

l.xperimental  evidence  for  the  field  guided  mode  was  obtained  from 
elevation  angle  and  relative  travel  time  measurements.  The  experimental 
travel  time  difference  measurements  agree  with  differences  between  the 
calculated  travel  time  for  the  numerous  modes  oocurlng  over  the  trans- 
equatorial  path  and  the  calculated  travel  time  for  the  field  guided  mode. 

Figure  3  shows  a  field  guided  path  and  two  trapezoidal  paths  and  several  dipole 
magnetic  field  lines  along  the  experimental  magnetic  meridian.  The  trapezoidal 
mode  occur  ed  regularly  over  the  tr  ana  equate  rial  path  but  did  not  display  the 
high  fading  frequencies  of  the  field  guided  mode. 

The  occurance  of  this  mode  at  2300  local  time  colncldea  with  equatorial 
topalte  spread  F  which  displays  magnetic  field  aligned  irregularity  structure  and 
occurs  at  a  time  when  h^F^  (dotted  line  in  the  figure)  is  decreasing  to  below  the 
height  of  the  field  guided  path. 


III-20-2 


SIGNAL  STRENGTH  .  dbm  SIGNAL  STRENGTH  ,  dbm 


LOCAL  TIME  AT  PATH  MIDPOINT 


60 

70 

80 


MAR. 17,1966 
-  FIJI 


MAR. 16, 1966 


LOCAL  TIME  AT  PATH  MIDPOINT 

Figure  1,  Diurnal  variation  of  signal  level  at  Darwin  and  Fiji  for  the  period 
0800  March  16  to  0800  March  17. 
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